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Requirements  for  the  Degree  of  Doctor  of  Philosophy 

A STUDY  OF  NOISE  PROCESSES  IN  TWO-DIMENSIONAL  AND 
THREE-DIMENSIONAL  ELECTRON  GAS  DEVICES 

By 

Arthur  Dirk  van  Rheenen 
May  1937 

Chairperson:  G.  Bosnian 

Major  Department:  Electrical  Engineering 

The  three  types  of  electrical  noise,  observed  in  most  electron 
devices,  are  investigated  in  different  semiconductor  device  structures. 
Emphasis  was  placed  on  hetero junction  devices  that  exhibit  charge  trans- 
port in  a quasi-two-dimensional  momentum  space. 

For  a pseudomorphic  (AlGaAs/InGaAs/GaAS)  modulation-doped  field- 
effect  transistor  (MODFET),  noise  measurements  were  performed  in  the 
frequency  range  from  1 Hz  to  100  kHz.  The  observed  spectral  density  of 
the  current  fluctuations  as  a function  of  the  gate  voltage  shows  a 
strong  1/f  component  for  a gate  voltage  of  -94  mV  and  the  occurrence  of 
generation-recombination  (g-r)  components  at  other  gate  voltages. 
Explanations  are  presented  in  terms  of  a resonance  in  the  1/f  noise  and 
in  terms  of  the  Fermi  level  moving  through  the  quantization  levels  in 
the  conduction-band  well. 

Low-frequency  g-r  noise  measurements,  1 Hz  - 1 MHz,  performed  on 
planar  silicon  resistors  were  used  to  demonstrate  a method  developed 
here  to  extract  the  parameters  that  describe  the  electrical  nature  of 
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carrier  trapping  centers.  The  trap  parameters  were  used  as  input 
parameters  of  an  exact  multitrap  model  and  yielded  zero— frequency  pla- 
teau levels  and  characteristic  times  in  excellent  agreement  with  the 
measured  data. 

Computer  simulations  of  the  trajectory  of  an  electron  in  wave- 
vector  space  under  the  influence  of  an  applied  electric  field,  and 
suffering  many  collisions  with  different  types  of  phonons,  resulted  in 
drift-velocity  field  curves  and  diffusion  coefficient  field  curves.  A 
high  mobility  was  found,  much  larger  than  observed  in  experiments. 
High-frequency  (0.5  - 12  GHz)  noise  measurements  on  gateless  AlGaAs/GaAs 
MODFETs  only  partially  supported  the  theoretical  results  for  the  diffu- 
sion coefficient.  It  is  argued  that  the  discrepancy  results  from  the 
high-mobility  value  found  in  the  simulations. 
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CHAPTER  I 
INTRODUCTION 


One  of  the  most  severe  limitations  of  the  operation  of  semicon- 
ductor devices  is  the  fluctuation  of  the  output  signal  due  to  internally 
generated  noise.  Even  in  the  case  of  a noise-free  input  signal,  the 
output  signal  will  consist  of  the  sum  of  the  processed  input  signal  and 
the  added  noise  signal.  When  the  ratio  of  the  processed  input  signal 
power  to  the  generated  noise  power  in  the  device  goes  to  one,  the  output 
signal  will  not  be  an  accurate  representation  of  the  input  signal  any 
more.  In  this  way  the  noise  generated  by  the  device  limits  the  sensi- 
tivity of  the  device. 

It  has  always  been  the  goal  in  noise  studies  to  identify  possible 
sources  of  the  noise.  Once  identified  and  understood,  methods  might  be 
developed  to  reduce  the  noise.  For  example,  better  processing  tech- 
niques might  be  contrived,  or  regions  of  operation  of  the  devices  can  be 
found  where  the  produced  noise  is  minimal. 

The  electrical  noise  of  a device  is  measured  in  terms  of  the  power 
spectral  density  of  voltage  or  current  fluctuations.  Usually,  three 
types  of  noise  are  found  in  the  resulting  spectra,  each  of  which  has  its 
own  specific  frequency  dependence.  In  the  frequency  range  from  1 Hz  to 
100  kHz  the  two  competing  types  are  1/f  and  generation-recombination 
(g-r)  noise.  The  power  spectral  density  of  the  former  is  proportional 
to  1/f  where  f is  the  frequency.  Although  this  type  of  noise  has  been 
identified  in  many  physical  systems,  its  origin  is  highly  controversial. 
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Three  different  mechanisms  have  been  proposed:  a summation  of  g-r 

spectra  [1],  mobility  fluctuations  [ 2 , 3 J , and  fluctuations  of  the  cross- 
sections  of  elect ron— phonon  interactions  [4].  The  g— r noise— producing 
physical  processes  are  well  understood.  The  trapping  and  detrapping  of 
free  carriers  in  semiconductors  by  discrete  energy  levels  in  the  for- 
bidden band  give  rise  to  fluctuations  in  the  number  of  free  carriers  [5] 
and  hence  in  the  current  conducted  by  the  carriers.  The  noise  power 
spectrum  has  a typical  frequency  dependence:  S (0)/(l  + oTx2).  S1(0) 

is  the  zero— frequency  plateau  value,  w = 2mf,  and  x is  the  character- 
istic time  which  is  equal  to 


T " <2’f-3dB> 


(1.1) 


At  higher  frequencies  the  main  contribution  to  the  noise  power 
spectrum  comes  from  fluctuations  of  the  velocity  of  the  carriers  parti- 
cipating in  the  conduction  process.  The  fluctuations  in  the  carrier 
velocity  are  caused  by  scattering  of  that  carrier  with  lattice  vibra- 
tions, other  carriers,  etc.  The  resulting-  spectrum  is  independent  of 
frequency,  up  to  frequencies  comparable  to  the  reciprocal  average  inter- 
collision  time. 

As  is  indicated  above,  the  measured  noise  can  be  related  to  physi- 
cal microscopic  processes  that  pertain  to  the  charge  transport  in  the 
devices  under  consideration.  Therefore,  also  from  a theoretical  point 
of  view,  noise  measurements  are  of  particular  interest  since  a better 
understanding  of  these  processes  will  result. 

In  this  dissertation  the  three  types  of  noise  are  studied  in  great 
detail  as  it  is  measured  in  different  semiconductor  materials  and  device 
structures.  In  Chapter  II  1/f  noise  measurements  are  presented  that 
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were  performed  on  an  AlGaAs/InGaAs  pseudomorphic  modulation-doped  field- 
effect  transistor  (MODFET).  The  g-r  noise  measurements,  performed  on 
planar  silicon  resistors,  are  discussed  in  Chapter  III.  In  Chapter  V 
the  velocity-fluctuation  noise  measurements,  performed  on  a gateless 
AlGaAs/GaAs  MODFET,  are  presented. 

!• 1*  Band  Structure  of  AlGaAs/GaAs  and  AlGaAs/InGaAs/GaAs 

Hetero junctions 

The  difference  between  a hetero junction  and  a homojunction  is  that 
in  the  former  case  the  junction  is  formed  by  two  materials  with  dif- 
ferent energy  bandgaps,  while  in  the  latter  case  the  junction  is  formed 
by  identical  semiconductor  material  but  with  a different  doping  on 
either  side  of  the  junction.  In  this  dissertation  many  experiments  are 
described  which  are  performed  on  hetero junction  devices.  Therefore,  the 
band  structure  of  these  devices  will  be  reviewed  here. 

When  two  semiconductor  materials  with  different  energy  bandgaps 
are  brought  together,  discontinuities  in  the  conduction  band  and  valence 
band  at  the  interface  will  be  the  result.  The  way  in  which  the  total 
bandgap  difference  is  distributed  over  the  conduction-  and  valence-band 
discontinuities  is  still  under  discussion.  Initially,  Dingle  [6]  sug- 
gested that  for  the  AlGaAs-GaAs  system  85%  of  the  bandgap  difference 
would  be  accommodated  in  the  conduction— band  discontinuity.  However, 
lately  Nussbaum  [7]  suggested  that  the  conduction-band  discontinuity  is 
equal  to  65%  of  the  bandgap  difference.  In  the  case  of  a homo junction, 
the  difference  in  the  conduction-band  energies,  far  away  from  the  inter- 
faces, is  just  the  difference  between  the  positions  of  the  Fermi  levels 
in  the  two  materials  before  they  were  brought  together.  In  thermal 
equilibrium  the  Fermi  levels  have  to  line  up  (dEF/dx  = 0),  and  the  con- 


ductlon  and  valence  bands  are  shifted  with  respect  to  each  other  until 
the  Fermi  levels  do  line  up. 


For  heterojunct ions  the  shape  of  the  conduction  and  valence  bands 
is  governed  by  the  difference  between  the  electron  affinities  and  elec- 
tron work  functions  of  the  two  materials.  In  addition,  the  Fermi  level 
has  to  line  up.  The  difference  between  the  electron  affinities  is 
supposed  to  cause  the  conduction-band  discontinuity  at  the  interface. 
This  affinity  rule  AE^  = Ax  may  not  be  a valid  assumption.  However, 
choosing  AE^  as  an  adjustable  parameter,  the  Anderson  [3]  model  remains 
unaltered  [9],  The  difference  between  the  work  functions  determines  how 
much  the  bands  have  to  bend  to  allow  the  Fermi  levels  to  be  lined  up. 
Figure  1.1  illustrates  the  energy  band  diagram  of  a hetero junction 
between  n-type  AlGaAs  and  intrinsic  GaAs.  The  bending  of  the  bands  is 
caused  by  an  electron  transfer  from  the  larger  work-function  material  to 
the  smaller  work-function  material.  The  separation  of  charges  gives 
rise  to  an  electric  field  normal  to  the  interface.  In  the  positive 
space-charge  region  (uncovered  donors)  the  conduction  band  bends  con- 
cavely , and  in  the  negative  space-charge  region  (electrons)  the  band 
bends  convexly.  This  follows  directly  from  solving  Poisson's  equation, 


q dx2  dx 


£ 

e * 


(1.2) 


where  q is  the  charge  of  an  electron,  Ec  the  conduction-band  edge 
energy , F the  electric  field  strength,  p the  space  charge,  e the 
dielectric  constant,  and  x is  taken  normal  to  the  interface. 

When  the  electric  field  is  strong  enough,  the  conduction  band  is 
bent  so  much  that  the  carriers  will  be  confined  in  a region  close  to  the 
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Fig.  1.1.  Heterojunction  band  lineup  according  to  the 
Anderson  [8]  model. 
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interface.  In  Fig.  1.2  one  interface  is  drawn  and  the  carrier  confine- 
ment is  suggested.  This  confinement  leads  to  high  sheet-carrier  densi- 
ties desirable  for  large  current  conduction. 

Classically,  this  confinement  can  be  understood  by  realizing  that 
the  built-in  field  is  forcing  the  electrons  against  the  interface  where 
the  conduction-band  discontinuity  keeps  them  from  transfering  into  the 
larger  bandgap  material.  A slightly  different  type  of  confinement  is 
realized  by  making  the  small  bandgap  material  layer  very  thin  and  sand- 
wiched between  two  larger  bandgap  materials,  e.g. , AlGaA/GaAs/AlGaAs  or 
AlGaAs /InGaAs/Ga As . In  this  case  the  confinement  is  physically 
suggested  by  the  well  in  the  conduction  band.  For  very  thin  layers 
(~  200  A),  quantum  mechanical  effects  will  determine  the  charge  trans- 
port in  what  is  called  a two-dimensional  electron  gas  (2DEG).  The  size 
limitation  in  the  direction  perpendicular  to  the  interface  will  result 
in  a quantization  of  the  wave  vector  in  this  direction.  Therefore,  a 
quantum  mechanical  analysis  is  needed  to  find  the  electron  distribution; 
the  time-independent  Schrodinger  equation  has  to  be  solved. 

= Et|»  , (1.3) 


where  H is  the  Hamiltonian,  iji  the  wave  function,  and  E the  energy.  How- 
ever, this  equation  needs  to  be  solved  in  conjunction  with  the  Poisson 
equation.  To  see  the  relation  more  clearly,  eq.  (1.2)  is  rewritten  as 


. d2E 
_1  c 

q dx2 


f "S0 


44 

00 

r * 

J 44  dx 


(1.4) 


where  nsg  is  the  total  sheet-carrier  concentration.  The  simultaneous 
solution  of  eqs.  (1.3)  and  (1.4)  is  rather  complicated.  Therefore, 
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Fig.  1.2.  Carrier  confinement  at  the  interface  as  a 
result  of  the  band  bending. 
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approximations  are  made  to  gain  physical  insight.  The  conduction  band 
is  either  approximated  by  a triangular  (Fig.  1.2)  well  or  a rectangular 
well  in  the  case  of  a thin,  small  bandgap  material  layer.  In  both  cases 
the  energy  levels  and  wave  functions  are  well  known  and  the  electron 
distribution  easily  evaluated.  In  Chapter  IV  electron  transport  is 
described  in  square  quantum  wells.  The  barrier  heights  are  assumed  to 
be  infinite. 

In  real  devices,  however,  the  barrier  height  is  of  the  order  of 
200  meV  and  thus  far  from  infinite.  Even  for  moderate,  applied  electric 
field  strengths  parallel  to  the  interface,  the  average  kinetic  energy  of 
the  electrons  will  become  comparable  to  the  barrier  height.  Hence, 
carriers  can  transfer  back  Into  the  large  bandgap  semiconductor.  This 
real-space  charge  transfer  gives  rise  to  the  negative  differential 
mobility.  Increasing  the  field  strength  will  not  lead  to  a larger 
current,  but  instead,  due  to  the  transfer  of  high-mobility  electrons 
Into  regions  of  low  mobility,  to  a smaller  current.  This  is  called  the 
real-space-charge-transfer  effect  [10]. 

An  analogous  effect  is  the  wave-vector-space-charge  transfer. 
This  effect  is  explained  by  considering  the  intrinsic  conduction  band 
profile  of,  for  instance,  GaAs.  The  direct  semiconductor,  GaAs , has  a 
central  conduction-band  valley  at  the  T point.  In  addition,  local 
minima  in  the  conduction  band  are  found  located  along  the  main  crystal- 
lographic axes  [(100),  (010),  and  (001)]  and  the  diagonals  [(111), 

(111),  (111),  and  (Til)].  The  first  ones  are  called  X valleys  and  the 
last  ones  are  called  L valleys.  The  curvature  of  the  satellite  valleys 
in  k space  is  different  from  the  central  valley;  In  other  words,  the 
effective  mass  Is  different  in  the  different  valleys.  The  effective 


9 


mass  is  the  smallest  in  the  central  valley.  Hence  the  transfer  of  an 
electron  from  the  central  to  the  satellite  valley  will  reduce  its  mobil- 
ity, and  a negative  differential  mobility  is  observed.  As  the  electric 
field  is  increased,  the  electrons  gain  kinetic  energy  to  overcome  the 
energy  offset  between  the  central  and  satellite  valleys  and  they  trans- 
fer to  these  valleys,  thus  reducing  the  conductivity.  Figure  1.3  shows 
a schematic  band  diagram  for  GaAs . This  effect  is  called  the  RWH 
[11-14]  or  Gunn  [15,16]  effect. 

1.2.  Noise  and  Charge  Transport 

Measurements  of  the  noise  in  electrical  currents  passed  through  a 
semiconductor  can  be  used  to  obtain  parameters  that  are  of  importance  in 
the  charge  transport.  G-r  noise  measurements  are  used  to  find  the 

pertinent  parameters  that  describe  trapping  centers  (see  Chapter  III). 
The  lifetimes  of  the  carriers  are  inversely  proportional  to  the  concen- 
tration of  recombination  centers,  and  in  most  device  applications  the 
carr^er  lifetime  plays  an  important  role.  Therefore,  it  is  very 
important  to  characterize  the  impurities,  lattice  defects,  etc.,  which 
are  introduced  into  the  lattice  during  fabrication  of  the  device  and 
affect  the  lifetime.  The  only  widely  accepted  method  of  characteriza- 
tion of  deep  traps  is  deep-level  transient  spectroscopy  (DLTS) . In  this 
method  the  capacitance  of  the  space-charge  region  of  a reverse-biased 
p-n  junction  is  measured,  and  hence  large  electric  fields  are  present  in 
the  device  which  might  affect  the  dynamics  of  the  trap  centers.  The 
noise  measurements,  however,  are  performed  in  the  ohmic  device  operating 
regime  and  do  not  require  a p-n  junction.  The  noise  can  be  measured 
in  any  part  of  a device  without  special  preparation  of  the  specimen. 
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Fig*  1.3.  Schematic  band  diagram  of  GaAs  to  illustrate 
the  wave-vector-space-charge  transfer. 
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Hence,  g-r  noise  measurements  are  a valuable  tool  to  characterize  traps 
in  all  kinds  of  semiconductor  devices. 

The  velocity-fluctuation  noise  is  also  related  to  the  transport  of 
charge.  An  electron  traveling  in  the  channel  of  a MOSFET  or  a MODFET 
scatters  several  times  with  phonons,  impurities,  and  other  carriers 
before  it  reaches  the  drain  contact.  In  the  event  of  such  a scattering 
the  wave  vector  of  the  electron  changes  instantaneously,  and  therefore 
its  momentum  is  changed.  Consequently,  a scattering  event  gives  rise  to 
a fluctuation  in  the  velocity  of  the  carrier.  Since  the  current  is  pro- 
portional to  the  velocity  of  the  carriers,  it  is  obvious  that  a velocity 

fluctuation  will  give  rise  to  a current  fluctuation.  A measure  for  the 

spectral  density  of  the  velocity  fluctuations  is  the  diffusion  coeffi- 

cient. To  illustrate  this,  suppose  a sharp  pulse  of  carriers  is  in- 
jected into  the  source  region  of  the  channel.  If  all  of  the  carriers 
travel  at  the  same  speed,  without  being  scattered,  then  at  the  drain 
side  the  same  pulse  will  be  collected.  However,  if  the  carriers  suffer 
collisions  during  their  travel,  the  velocity  distribution  function  will 
widen  and  the  collected  pulse  at  the  drain  will  be  smeared  out.  This 
widening  of  the  velocity  distribution  function  is  a direct  measure  for 
the  diffusion  coefficient.  In  formula: 

00 

Sv(w)  = 4D(u>)  = / <Av(t)Av(t+x)>e  ^“Tdx  , (1.5) 

0 

where  D is  the  diffusion  coefficient  and  $(x)  = <Av(t )Av(t+x)>  the  auto- 
correlation function  of  the  velocity  fluctuations: 

1 T 

*(t)  = <Av(t)Av(t+r)>t=  - / Av(t)Av(t+x)dt  . (1.6) 
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By  definition: 

CO 

tq*$(0)  = / $(x)dT  , (1.7) 

0 

which  defines  the  correlation  time  x . From  eqs.  (1.5-1. 7)  it  is  clear 
that  the  low-frequency  limit  of  the  diffusion  coefficient  is  determined 
by  the  product  of  the  correlation  time  and  the  variance  of  the  velocity 
fluctuations.  For  randomizing  scattering  processes  (optical  phonon 
interactions)  the  memory  effect  will  be  small;  in  other  words,  the 
correlation  time  will  be  short.  However,  the  randomizing  scatterings 
involve  relatively  large  wave-vector  changes  and  hence  velocity  changes. 
Therefore,  the  variance  $(0)  of  the  velocity  fluctuations  will  be  rather 
large.  For  nonrandomizing  collisions,  such  as  polar  optical  phonon 
interactions,  which  favor  scattering  over  small  angles,  the  memory 
effect  will  be  long  or  the  correlation  time  will  be  large.  However,  the 
variance  of  the  velocity  fluctuations  will  be  small  since  small  angle 
scattering  involves  only  small  wave-vector  (velocity)  changes.  From 
these  observations  it  is  clear  that  the  two  factors  that  determine  the 
value  of  the  diffusion  coefficient  compete,  and  the  ultimate  test  to 
determine  the  influence  of  the  scattering  mechanisms  on  the  diffusion 
coefficient  will  be  the  experiment  or  an  extensive  Monte  Carlo  simula- 
tion. 

It  is  indicated  that  the  diffusion  coefficient  is  sensitive  to  the 
type  of  scattering  processes  the  charge  carriers  undergo  while  traveling 
from  source  to  drain.  A measurement  of  the  diffusion  coefficient  pro- 
vides information  essential  to  the  transport  of  charge.  The  tlme-of- 
f light  [17,18]  experiment  will  yield  the  dc  diffusion  coefficient. 
Noise  measurements  can  be  used  as  a tool  to  obtain  the  diffusion  coeffi- 


1.3 


cient  as  a function  of  frequency.  The  current  can  be  written  as 

i = qnvA,  (1.7) 

where  n is  the  density  of  carriers,  v the  velocity,  and  A the  cross- 
sectional  area  perpendicular  to  the  current  flow.  Then  it  follows,  when 
we  neglect  number  fluctuations,  that  the  spectral  density  of  the  current 
fluctuations  equals 

S1(u>  = Sv(a))  = D(u)  * (1.8) 

Here,  N is  the  total  number  of  carriers  and  L the  length  of  the  device. 
It  is  shown  here,  eq.  (1.8),  that  a measurement  of  the  spectral  density 
of  the  current  fluctuations  contains  direct  information  pertaining  to 
the  diffusion  coefficient. 

The  relation  between  the  1/f  noise  and  the  charge-transport  param- 
eters is  still  obscure.  It  has  been  suggested  by  Handel  [4]  that  the 
magnitude  of  the  1/f  noise  is  related  to  the  cross  section  of  the  scat- 
tering mechanisms.  Again,  noise  measurements  would  yield  information 
related  to  the  charge  transport.  Mobility  fluctuations  are  also  named 
[2]  as  a possible  origin  for  1/f  noise. 


CHAPTER  II 

1/F  NOISE  IN  PSEUDO-MORPHIC  (ALGAAS/INGAAS) 
MODULATION-DOPED  FETS 

2.1.  Introduction 

Modulation-doped  heterojunct ion  semiconductor  structures  have  been 

shown  [19]  to  be  excellent  candidates  for  field— effect  transistors 

(MODFETs).  The  vast  majority  of  the  reported  MODFETs  utilize  A1  5,  As/ 

x 1-x  ' 

GaAs  hetero junctions  in  which  electrons  are  transferred  from  the  larger 
bandgap  material  (A1xG^_xAs)  t0  smaller  bandgap  semiconductor 

(GaAs).  A very  important  reason  for  investigating  this  structure  made 
of  (Al,Ga,As)  is  that  it  is  relatively  easy  to  grow  perfect  lattices 
with  excellent  lattice  matching.  The  lattice  constant  of  GaAs  is  very 
close  to  the  lattice  constant  of  AlAs , as  it  is  for  every  mixture  of 
(Al,Ga)As,.  and  no  strain  is  introduced  in  the  lattice.  However,  there 
are  two  serious  disadvantages  associated  with  the  AlGaAs /GaAs  system: 
the  persistent  photoconductivity  effect  [20]  and  the  collapse  of  the 
drain  current  [21]  at  lower  temperatures.  These  two  effects  are  attri- 
buted to  the  introduction  of  deep  trap  levels  (DX  centers).  The  concen- 
tration of  these  trap  levels  increases  as  the  A1  mole  fraction  increases 
in  A1xGai_xAs.  For  a proper  operation  of  the  MODFET,  a high  sheet 
carrier  concentration  is  required  at  the  AlGaAs /GaAs  interface,  or,  in 
other  words,  the  electrons  that  spill  over  into  the  GaAs  from  the  AlGaAs 
must  be  confined  in  the  GaAs.  For  this  purpose  the  conduction  band  dis- 
continuity must  be  as  large  as  possible,  to  keep  the  carrier  from  trans- 
ferring back  into  the  AlGaAs.  The  conduction  band  discontinuity  is  pro- 
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portional  to  the  A1  mole  fraction  x up  to  x = 0.3  when  AE  = 200  meV. 
- c 

There  is  a tradeoff  between  having  a well-confined  electron  gas  and 
hence  a large  sheet  carrier  concentration  which  is  capable  of  carrying 
large  currents  on  the  one  hand,  and  on  the  other  hand  the  two  degrading 
effects  mentioned  before. 

It  has  been  predicted  [22]  and  demonstrated  [23]  that  a lattice- 
mismatched  layer  can  be  grown  sufficiently  thin  for  the  mismatch  to  be 
accommodated  entirely  as  elastic  strain.  In  this  situation  the  inter- 
face between  the  materials  is  essentially  free  from  misfit  dislocations, 
and  the  thin  layer  is  pseudomorphic . This  technology  opens  the  door  for 
heterojunction  systems  other  than  AlGaAs/GaAs  provided  only  a thin  layer 
of  one  of  the  materials  is  needed.  In  reference  [24]  current -voltage 
characteristics  and  values  for  the  transconductance  are  presented  for 
HEMT  s which  utilize  an  InQ . 1 5Ga0. 85As/GaAs  heterojunction  in  which 
electrons  are  transferred  from  the  heavily  doped  GaAs  to  the  narrower 
bandgap  InQ> 15Ga085As  (which  is  grown  pseudomorphically  on  a GaAs  sub- 
strate). The  A^Ga^As  is  eliminated  from  the  device  and  no  persistent 
photoconductivity  or  drain-current  collapse  is  observed." 

Ketterson  et  al.  [25]  replaced  the  GaAs  by  the  larger  bandgap 
material  AlGaAs , hence  increasing  the  sheet  carrier  concentration  at  the 
AlGaAs/InGaAs  interface.  They  report  that  "device  results  obtained 

using  this  InGaAs/AlGaAs  system  rival  the  best  reported  GaAs/AlGaAs 
MODFET  results.” 

In  this  chapter  low-frequency  spectral  density  measurements  of  the 
drain-current  fluctuations,  which  were  performed  on  the  InGaAs/AlGaAs 
structures,  are  presented.  The  noise  spectra  were  measured  as  a func- 
tion of  the  applied  (dc)  gate  voltage  and  the  (dc)  drain  current.  The 
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results  are  presented  in  section  2.5,  where  they  will  also  be  discussed. 
In  the  next  section,  2.2,  the  device  used  in  the  experiments  is 
described  and  a sketch  of  the  band  structure  is  shown.  A detailed 
account  of  the  dc  characterization  is  given  in  section  2.3,  and  the 
noise  measurement  setup  is  discussed  in  section  2.4. 

2.2.  Device  Description 

The  pseudomorphic  single  quantum  well  InCaAs/AlGaAs  structure 
studied  was  supplied  by  H.  Morkoc  of  the  University  of  Illinois.  The 
structure  of  the  MODFET  is  sketched  in  Fig.  2.1.  Using  molecular  beam 
epitaxy  (MBE) , a 200  A-thick  layer  of  intrinsic  InQ  < 15GaQ  ^As  is  grown 
on  top  of  a semi -insulating  GaAs  substrate.  Then  a 30  A intrinsic 
AlGaAs  spacer  layer  is  grown  followed  by  a 350  A A1Q  15  GaQ  354s  layer 
doped  with  3 x lO1^  Si  atoms/cm^.  This  layer  supplies  the  electrons 
that  form  the  two-dimensional  electron  gas  ( 2DEG) , the  hatched  region  in 
?ig.  2.1.  The  top  layer  is  an  n+-GaAs  layer  to  facilitate  the  formation 
of  an  ohmic  contact.  The  source  to  drain  spacing  is  4 ym  and  the  gate 
length  is  1 yn.  The  gate  width,  (V),  is  145  ym. 

A typical  band  diagram  for  this  structure  is  sketched  in  Fig.  2.2. 
Due  to  the  separation  of  the  ionized  donors  in  the  n-AlGaAs  and  the 
mobile  electrons  in  the  InGaAs , an  electric  field  builds  up  across  the 
AlGaAs /InGaAs  interface.  This  electric  field  gives  rise  to  a bending  of 
the  energy  bands.  Of  interest  here  is  the  bending  of  the  conduction 
band  (electron  transport  only).  By  applying  a positive  voltage  to  the 
gate,  the  built-in  field  is  aided,  and  there  will  be  more  free  carriers 
in  the  InGaAs.  Likewise,  a more  negative  gate  voltage  will  result  In  a 
reduction  of  the  number  of  mobile  carriers  in  the  InGaAs.  Consequently, 
the  gate  voltage  modulates  the  number  of  free  carriers  In  the  InGaAs, 
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Fig.  2.1.  Device  structure  (S  = source,  D = drain, 
G = gate). 
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Fig.  2.2. 


Conduction  band  diagram  for  the  structure 
in  Fig.  2.1. 
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and  therefore  the  conductivity  of  the  channel.  The  doping  of  the 
n-AlGaAs  Is  so  high  that  the  Fermi  level  (Ep)  lies  within  the  conduction 
band  of  the  InGaAs,  resulting  in  a degenerate  electron  gas. 

The  small  thickness  of  the  InGaAs  layer  (200  A)  results  in  a 
quantization  of  the  wave  vector  of  the  electrons  in  the  direction  per- 
pendicular to  the  interface.  This  makes  the  conduction  band  of  the 
InGaAs  split  up  into  subbands.  The  bottom  of  each  of  these  subbands  Is 

the  lowest  energy  state  of  a continuum;  it  is  not  a discrete  energy 
level. 

Due  to  the  built-in  electric  field  the  electrons  tend  to  cluster 
in  the  vicinity  of  the  heterointerface.  This  is  also  suggested  in 
r'ig.  2.2  by  the  shape  of  the  conduction  band,  which  reaches  its  lowest 
point,  with  respect  to  the  Fermi  level,  at  the  interface.  However, 
quantum  mechanical  effects  will  smear  out  the  electron  distribution  over 
the  whole  width  of  the  well.  The  solution  of  the  Schrodinger  equation 
will  require  the  wave  functions  to  go  to  zero  outside  the  well.  For  the 
finite-height  well  described  here,  the  wave  functions  will  extend  some- 
what into  the  barriers.  Still,  effectively  there  will  be  a vanishing 
probability  of  finding  the  electrons  at  the  interface.  A second  effect 
that  keeps  the  electrons  from  clustering  at  the  AlGaAs/InGaAs  interface 
is  the  Pauli  exclusion  principle.  A clustering  of  electrons  at  the 
interface  would  require  the  occupation  of  higher  energy  levels  which 
would  be  energetically  disadvantageous,  resulting  in  a spatial  redis- 
tribution of  the  carriers.  In  summary,  both  effects  will  yield  a finite 
width  of  the  spatial  electron  distribution  function.  The  hatched  region 
in  Fig.  2.1  suggests  this. 
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2.3.  Pc  Characterization 

The  main  objective  of  the  dc  characterization  is  to  describe  the 
properties  of  the  active  region  of  the  pseudomorphic  MODFET.  The  active 
region  is  the  region  under  the  gate.  Once  its  characteristics  are 
established,  the  source-to-drain  channel  can  be  modeled  as  three  lumped 
resistors  in  series,  see  Fig.  2.3.  The  resistors  Rss  and  describe 
the  resistance  between  the  source  and  the  active  region,  and  the  resis- 
tance between  the  active  region  and  the  drain,  respectively.  Unlike  RQ<, 
and  Rdd,  the  active  region  resistance  is  a voltage  controlled  one.  As 
is  explained  in  the  previous  section,  the  gate  voltage  controls  the 
electron  density  in  the  InGaAs  channel,  but  only  under  the  gate  area. 
The  first  step  in  the  dc  characterization  is  the  measurement  of  the 
current -voltage  (I-V)  characteristics  as  a function  of  the  gate  voltage 
(VG).  The  basic  setup  is  described  in  Fig.  2.4.  With  resistor  Rj  the 
gate  voxtage  is  controlled  and  with  R2  the  drain  current  is  controlled. 
In  a semi -automated  setup  using  the  Tektronix  digital  scope  driven  by  a 
Tektronix  desktop  computer,  the  voltage  VQS  and  the  voltage  across  the 
resistor  R2  are  measured,  and  these  data  are  converted  in  a point  of  the 
I-V  curve.  The  noise  measurements  were  to  be  performed  in  the  ohmic 
regime  (linear  part  of  I-V  characteristics)  to  avoid  excessive  Joule 
heating  and  simplify  the  analysis.  For  this  reason  only  the  linear  part 
of  the  I-V  characteristic  was  measured.  In  Fig.  2.5  these  results  are 
presented  for  the  gate  voltage  range  of  -0.5  V < VQ  < +0.2  V in  steps  of 
0.1  V. 

The  next  step  in  the  characterization  is  the  establishment  of  the 
relation  between  the  sheet  carrier  concentration  (ns)  in  the  InGaAs  well 
and  the  gate  voltage.  The  relationship  is  given  in  eq.  (2.3)  and 
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Fig.  2.3.  Lumped  resistor  model  for  the  source-to-dratn 

channel  in  a MODFET.  The  resistance  of  the  active 
part  is  controlled  by  the  gate  voltage. 
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2.4.  Setup  to  measure  the  channel  I-V  characteristic 
with  the  gate  voltage  as  a parameter. 
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DRAIN-SOURCE  VOLTAGE  (mV) 


Fig.  2.5.  Current  voltage  characteristics  in  ohmic  region. 


The  charge  Q on  the 
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is  derived  from  a simple  capacitance  argument, 
piste  of  a capacitor  G is  given  by 

Q = G • V 


(2.1) 


where  V Is  the  voltage  across  the  plates.  The  gate  control  of  the 

channel  conductivity  is  modeled  as  a capacitive  control.  The  charge  per 
unit  area  is 


Q'  = C' 


(2.2) 


where  the  prime  indicates  the  per  unit  area  quantity.  The  charge  per 
unit  area  in  the  channel  is  qn3,  while  the  capacitance  per  unit  area 

iS  £0er//d  ’ where  eQ  and  £r  are  the  vacuum  and  relative  dielectric  con- 
stant, respectively,  and  d is  the  distance  between  the  gate  and  the 
2DEG,  which  is  about  500  A,  see  Fig.  2.1.  Thus  eq.  (2.2)  reads 


qns  ■ -p-  <vv 


(2.3) 


where  VT  is  the  threshold  voltage  for  the  onset  of  introducing  free 
carriers  in  the  InGaAs  channel. 

To  find  the  contact  or  gate  access  resistances  R^  and  RDQ,  a plot 
of  the  total  resistance  of  the  drain-to-source  channel  (R  ) versus 

L/O 

<VV  1 has  t0  be  raade,  as  will  be  shown.  As  can  be  seen  from 
Fig.  2.3, 


RDS  = RSS  + Ract  + RDD 


Here,  Ract  is  modeled  as 


R 


act 


_ _l i 

A qunA  qungW  ’ 


(2.4) 


P 


(2.5) 
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where  p is  the  resistivity,  £ the  length,  and  A the  cross-sectional  area 
perpendicular  to  the  current  flow.  The  electron  charge  is  q,  the 
mobility  is  u,  and  n is  the  number  of  carriers  per  unit  volume. 

It  follows  from  eqs.  (2.3),  (2.4),  and  (2.5)  that 

“ RSS  + W * Vr<VVT)  + *0°  * (2'6) 

Hence,  a plot  of  versus  (V^V^y1  will  yield  Rss  + RQD  as  the  inter- 
cept of  the  Rds  axis.  The  slope  of  that  linear  graph  will  give  d, 
assuming  that  u is  known,  providing  a check  of  the  analysis,  since  d is 
also  found  from  the  device  dimensions  (Fig.  2.1). 

The  only  parameter  that  has  to  be  determined  yet  is  the  threshold 
voltage,  V,p.  For  FETs  the  saturation  current,  the  constant  part  in  the 
I V curve,  is  related  to  the  gate  voltage.  The  dependence,  however, 
distinguishes  two  cases:  the  so-called  long  channel  approximation  and 
the  so-called  short  channel  approximation  [26].  In  the  former  approxi- 
mation the  current  saturation  is  caused  by  the  depletion  of  the  con- 
ducting channel  of  charge  carriers;  i.e.,  the  channel  is  pinched  off. 
^D(sat)  proportional  to  (Vq-V^,)“.  In  the  latter  approximation 
the  current  saturation  is  caused  by  the  saturation  of  the  drift  velocity 
of  the  carriers.  In  this  case  ID(sat)  is  proportional  to  (Vq-Vt) . 

In  both  cases  the  threshold  voltage  can  be  found  from  plotting 

^(sat)  °r  ID(sat)  3S  3 function  of  Che  gate  voltage.  The  complete  I-V 
characteristics  were  measured  by  Ketterson  et  al.  [27]  and  reproduced 
here  in  Fig.  2.6.  The  next  two  figures,  2.7  and  2.8,  show  the  two 
curves  for  and  1^^,  respectively,  as  a function  of  the  gate 
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Fig.  2.6.  Current  voltage  characteristics 
measured  by  Ketterson  [27]. 
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Fig.  2.7.  Long  channel  approximation.  The  gate  voltage 
intercept  is  the  threshold  voltage  VT. 
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Fig.  2.8.  Short  channel  approximation.  The  gate  voltage 
intercept  is  the  threshold  voltage  VT> 
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voltage.  The  intercepts,  with  the  VG  axis,  are  shown  in  both  cases; 
they  are  the  threshold  voltages  V^. 

From  Fig.  2.6  one  finds  that  the  current  saturation  occurs  for 
VDS  > 1 V‘  For  a 2ate  length  of  1 pm,  assuming  that  all  of  the  applied 
voltage  drops  across  the  active  region  under  the  gate,  the  average 
electric  field  strength  is  equal  to  10  kV/cm.  For  these  field  strengths 
the  drift  velocity  is  saturated.  It  is  clear  from  this  consideration 
that  the  short  channel  approximation  (Fig.  2.3)  applies  to  the  device  we 
are  concerned  with.  The  resulting  threshold  voltage  was  found  to  be 

0.37  V (see  Fig.  2.3).  As  pointed  out  before,  with  the  knowledge  of 
the  value  of  VT  the  access  resistances  (R$s  + RDD)  can  be  found  from  a 
plot  of  Rds  versus  (VG-VT)_l,  see  eq . (2.6).  The  intercept  of  the 

resistance  axis  in  Fig.  2.9  is  found  to  be  17  ft,  which  is  equal  to  the 
sum  of  Rgg  and  RDQ. 

Utilizing  the  preliminary  analyses  presented  thus  far,  the  road  is 
open  to  establish  the  quantitative  relationship  between  the  sheet- 
carrier  concentration  and  the  gate  voltage.  As  will  be  shown  in 
section  2.5,  this  relationship  will  play  a major  role  in  the  inter- 
pretation of  the  measured  noise  spectra  presented  there.  For  each  gate 
voltage  the  total  resistance  of  the  channel  is  found  from  the  slopes  of 

the  curves  in  Fig.  2.5.  Subtracting  (Rss+Rdd)  will  then  give  R 

act 

and  from  eq.  (2.5)  the  sheet -carrier  concentration  is  found,  assuming 

2 

U = 'Vile  GaAs  = 0,66  m /Vs  • These  results  are  plotted  in  Fig.  2.10. 
The  check  on  the  analysis  is  the  slope  of  the  straight  line  found  in 

Fig.  2.10.  It  was  mentioned  before  that  the  slope  is  a measure  for  the 
distance  between  the  physical  gate  and  the  2DEG,  see  eq.  (2.3). 
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Channel  resistance  RDS  versus  (VG-VT)_1.  The  inter 
cept  of  the  resistance  axis  is  equal  to  RDD  + Rss. 


Fig.  2.9. 
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Fig.  2.10.  Sheet-carrier  concentration  as  a function 
of  the  gate  voltage. 
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Apart  from  that  the  voltage  axis  intersect  should  be  VT.  From  the  slope 
it  is  found  that  d = 500  A.  Compare  this  with  the  physical  device 
dimensions;  then  d is  expected  to  be  the  sum  of  the  thickness  of  the 
intrinsic  AlGaAs  layer  (30  A),  the  doped  AlGaAs  layer  (350  A),  and  the 
n -GaAs  layer  (200  A).  These  three  thicknesses  add  up  to  580  A.  The 
thickness  of  the  cap  layer  might  have  been  reduced  by  etching  when  the 
gate  was  laid  down.  Also,  the  mean  distance  of  the  2DEG  from  the 
interface  might  not  be  zero.  Consequently,  it  can  be  concluded  that  the 
distance  d that  is  found  from  Fig.  2.10  is  in  the  right  range. 

From  the  intercept  of  the  voltage  axis,  VT  is  found  to  be  -0.42  V, 
which  is  not  too  far  off  from  the  value  of  -0.37  V we  assumed  earlier. 
In  conclusion,  it  seems  that  the  dc  analysis  is  consistent,  more  so 
because  it  can  be  shown  that  the  slope  of  Fig.  2.10  is  very  critical  on 
the  choice  of  V^. 

Experimental  Setup  for  Low-Frequency  Noise  Measurements 

The  frequency  range  of  the  noise  measurements  discussed  here  is 
from  1 Hz  to  25  kHz  in  general,  and  uP  to  200  kHz  for  certain  spectra. 
The  generic  setup  is  shown  in  Fig.  2.11.  The  two  power  supplies  are 
ordinary  car  batteries,  to  avoid  any  60  Hz  pickup  caused  by  line  power 
supplies.  By  changing  Rg  different  voltages  can  be  supplied  to  the  gate 
(g)  of  the  MODFET.  The  resistor  RD  limits  the  current  flowing  through 
the  channel,  provided  that  RD  » RDg . The  FET  is  hooked  up  in  a common 
source  (s)  configuration,  and  the  noise  is  measured  at  the  drain  (d). 
The  signal  is  fed  into  a low-noise  amplifier  (CFW  101)  and  Fourier 
transformed  by  the  HP  3582A  spectrum  analyzer.  The  device  is  mounted  in 
a CTI  CRYOGENICS  Model  21  refrigerator. 
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Fig.  2.11.  Setup  for  low-frequency  noise  measurements . 
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In  this  setup  there  are  four  unknowns:  the  noise  of  the  device, 

the  noise  of  the  amplifier,  the  gain  bandwidth  product  of  the  system, 
and  the  signal  magnitude  of  the  calibration  noise  generator.  To  solve 
for  this  system  of  unknowns,  four  measurements  are  performed,  three  of 
which  are  just  the  three  different  settings  of  the  switch.  The  fourth 
is  the  calibration  of  the  noise  generator.  From  these  four  measurements 
the  four  named  parameters  can  be  determined,  most  importantly  the  spec- 
tral density  of  the  current  fluctuations  in  the  source-drain  channel. 
In  this  manner  noise  spectra  were  measured  for  different  gate  voltages 
and  drain  currents.  It  was  made  sure,  always,  that  the  FET  operated  in 
the  linear  regime  of  the  I-V  characteristic.  The  continuous  bias  points 
during  the  noise  measurements  were  compared  with  the  earlier  measured, 
pulsed  I-V  curves  to  avoid  nonlinear  operation  and  Joule  heating. 

2.5.  Results  and  Discussion 

The  range  of  gate  voltages  for  which  the  spectral  density  of  the 
current  fluctuations  was  measured  is  -0.3  V < V < + 0.2  V.  For  each 
gate  voltage  the  spectra  were  measured  for  two  drain  currents,  0.56  mA 
and  0.77  mA.  A representative  sample  of  these  measured  spectra  is  shown 
in  Figs.  2.12  to  2.14.  The  most  notable  feature  of  these  graphs  is  the 
shape  of  the  spectra  as  a function  of  the  gate  voltage.  For  V,  = -94  mV 
(Fig.  2.13),  pure  1/f  noise  spectra  are  observed  for  over  four  decades 
of  frequency.  However,  when  the  gate  voltage  is  either  increased  or 
decreased,  spectral  components  related  to  a different  noise  mechanism, 
most  likely  generation-recombination  (g-r)  processes  (see  Chapter  III), 
perturb  the  1/f  dependence  of  the  spectral  density.  Two  examples,  one 
at  VG  = -0.3  V and  one  at  Vg  = 0.2  V,  are  shown  in  Figs.  2.12  and  2.14, 
respectively  to  illustrate  this  phenomenon. 


POWER  SPECTRAL  DENSITY  (A2s) 
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FREQUENCY(Hz) 


Fig.  2.12.  Noise  spectra  for  VQ  = -293  mV  and  two  dif- 
ferent drain  currents.  The  lines  show  the 
f proportionality. 
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POWER  SPECTRAL  DENSITY  (A*s) 
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Fig.  2.13.  Noise  spectrum  for  VG  = -94  mV  and  two 

different  drain  currents.  The  lines  show 
the  f 1 proportionality. 


POWER  SPECTRAL  DENSITY  (A2s) 
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Fig.  2.14.  Noise  spectrum  for  VQ  = +214  mV  and  two 

different  drain  currents.  The  line  shows 
the  f~  proportionality. 
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Although  the  origin  of  1/f  noise  is  still  unknown,  it  has  been 
attributed  to  mobility  fluctuations  [2,3]  and  cross-sectional  fluctua- 
tions of  electron  scattering  processes  [4].  The  generation- 
recombination  processes  are  well  understood.  In  Chapter  III  the  result- 
ing noise  spectra  are  calculated.  Generation-recombination  noise  is 
recognized  by  the  specific  frequency  dependence  S ( 0) / ( 1 + u2t2)  of  its 
power  spectral  density.  S(0)  is  the  zero-frequency  plateau  level  and 
t - (27rf_3d3)  where  f_3dB  is  the  frequency  for  which  the  power  is 
reduced  to  1/2  S(0).  T is  often  called  the  characteristic  time  and  can 
be  associated  with  the  microscopic  trapping  process  of  a carrier  by  an 
impurity,  lattice  imperfection,  etc. 

It  can  be  shown  that  the  g-r  noise  reaches  its  maximum  when  the 
Fermi  level  coincides  with  the  energy  level  associated  with  the  trap 
under  consideration.  A simple  argument  to  confirm  this  statement  goes 
as  follows.  When  the  Fermi  level  is  separated  by  more  than  4 k T from 
the  energy  level  considered,  either  above  or  below,  then  the  levels  are 
either  all  occupied  or  all  empty.  In  either  case  there  will  be  almost 
no  transitions  from  the  conduction  band  (n-type  semiconductor)  to  the 
trap  level  or  vice  versa.  Hence,  the  electrical  noise  which  is  directly 
proportional  to  the  spectral  density  of  the  number  fluctuations  asso- 
ciated with  the  trapping  and  detrapping  will  be  zero.  For  positions  of 
the  Fermi  level  In  close  vicinity  of  the  trap  level,  however,  the  trap 
levels  will  be  partly  occupied,  and  transitions  will  occur.  In  fact, 
the  noise  will  be  maximal  when  the  relative  occupancy  of  the  traps  is 
50%. 

In  the  system  considered  here,  the  MODFET,  the  sheet  carrier  con- 
centration in  the  conducting  InGaAs  channel  is  controlled  by  the  gate 
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voltage.  The  sheet  carrier  concentration  is  a direct  measure  for  the 
position  of  the  Fermi  level  with  respect  to  the  conduction  band.  3y 
varying  the  gate  voltage,  the  position  of  the  Fermi  level  with  respect 
to  the  conduction  band  is  altered. 

A second  feature  of  the  system  studied  here  is  the  quantization  of 
the  wave  vector  and  splitting  up  of  the  subbands  associated  with  that. 
Although  the  energy  levels  calculated  from  the  quantization  of  the  wave 
vector  parallel  to  the  AlGaAs/InGaAs  interface  do  not  identify  discrete 
levels  as  in  the  case  of  the  levels  introduced  in  the  bandgap  by  carrier 
traps,  it  is  the  lower  boundary  for  a new  set  of  available  electron 
states.  This  set  is  a continuum  of  states.  As  was  suggested  in 
Fig.  2.2,  the  electron  gas  in  the  InGaAs  channel  is  degenerate  and  the 
Fermi  level  lies  in  the  conduction  band. 

From  these  observations  it  follows  that  the  necessary  ingredients 
to  observe  the  increase  and  decrease  of  the  magnitude  of  the  g-r  noise 
components  are  present.  By  changing  the  gate  voltage,  the  Fermi  level 
moves  through  the  conduction  band.  As  it  approaches  one  of  the 
quantization  levels  a new  set  of  energy  states  is  available  for  the 
electrons,  and  transitions  will  occur  between  the  so-called  subbands. 

As  in  the  case  of  discrete  levels,  the  spectral  density  of  the  fluctua- 
tions in  the  occupancy  of  the  subbands  will  maximize  when  the  Fermi 
level  coincides  with  one  of  the  quantization  levels.  On  the  other  hand, 
the  magnitude  of  the  g-r  noise  components  is  expected  to  be  minimal  when 
the  Fermi  level  is  right  in  between  two  of  those  levels. 

The  relation  between  the  sheet  carrier  concentration  and  the  Fermi 
level  for  a two-dimensional  electron  gas  follows  from 
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ns  = / g2D(E)f(E)dE 


E 

c 


(2.7) 


where  g2D(E)  is  the  density  of  states  per  unit  energy  for  a two- 

dimensional  gas  (see  Chapter  IV)  and  f(E)  is  the  Fermi-Dirac  function. 

For  quantization  levels  E (n  = 1,  2.  n ) 

n ' max7 


divided  by  2tt  . E^,  is  the  Fermi  energy. 

The  energy  levels  En  ara  found  from  solving  the  Schrodinger  equa- 
tion for  the  potential  well  in  the  conduction  band.  The  shape  of  the 
well  is,  in  essence,  a square  one.  However,  the  barriers  on  either  side 
of  the  well  have  different  heights,  due  to  the  difference  between  the 
conduction  band  discontinuities  at  the  AlGaAs /InGaAs  and  InGaAs/GaAs 
interfaces,  see  Fig.  2.15. 

The  total  wave  function  is  split  up  into  three  parts  as  indicated 
in  Fig.  2.15: 

AlGaAs : ^(x)  = A exp(<1(x+L))  x < -L  (2.9a) 
InGaAs:  ^(x)  = B cos(kx+a)  -L  < x < L (2.9b) 
GaAs:  ^XII(x)  = C exp( -^(x-L) ) L < x (2.9c) 


n 


s 


(2.3) 


where  m 0.04  mQ  is  the  effective  mass  of  an  electron  in  InGaAs,  '<3  is 
Boltzmann's  constant,  T is  the  temperature,  and  ft  is  Planck's  constant 


where  ^ and  <3  are  the  decay  constants  in,  respectively,  the  AlGaAs  and 
the  GaAs,  and  k is  the  wave  vector  of  the  electron  in  the  well: 


ENERGY 
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Fig.  2.15.  Schematic  conduction  band  diagram  for  the 
system  discussed  here.  The  energy  levels 
En  are  the  eigenvalues  of  the  Schrodinger 
equation  for  this  system. 
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'l  = O (AE  -E ) / h Z ) 


.2,1/: 


k = (2m*E/h2)1/2 


(2.10) 

(2.11) 


k3  - (2„*(4EC3-E)/h2)1/2  (2. 12) 

where  E is  Che  energy  of  the  electron.  Note  that  due  to  the  difference 
in  the  barrier  heights  on  either  side  of  the  well,  the  wave  function 
will  not  be  symmetrical  around  * = 0.  This  is  expressed  in  the  choice 
of  the  wave  function  in  the  well  (^(x))  by  adding  a phase  angle  to  the 
argument  of  the  cosine.  The  continuity  requirements  of  the  wave 

function  and  its  derivative  with  respect  to  x at  the  interfaces,  i.e., 
at  x - ±L,  yield  four  equations  in  the  variables  A,  B,  C,  E,  and  a.  The 
condition  that  the  probability  of  finding  the  electron  anywhere  equals 

1 , i • 6 • y 


i M’O'AxM*  - 1 , (2.13 

supplies  the  fifth  equation,  and  the  system  can  be  solved.  After  re 
uriting  the  equations,  the  eigenvalue  equation  that  has  to  be  solved  is 
k3 

E-  ’ tanUkL  - 3)  (2.14; 

where  tan(S)  - Kj/k  . This  is  essentially  an  equation  in  one  variable, 
the  energy  E.  The  eigenvalues  are: 


Ej  = 17  meV 

= 67  meV 

E-j  = 147  meV 


where  the  following  parameter  values  have  been  used:  AE  = 400  raeV 

AEC3  = 250  meV*  ra  = 0.04  mQ,  and  2L  = 200  A . 
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Once  the  energy  eigenvalues  are  found,  eq.  (2.8)  can  be  plotted  as 
is  done  in  Fig.  2.16.  With  the  help  of  the  previously  established  rela- 
tionship between  the  sheet  carrier  concentration  and  the  gate  voltage 

(Fig.  2.10),  a gate  voltage  axis  is  plotted  in  Fig.  2.16.  Hence,  an 
empirical  relationship  is  found  between  the  gate  voltage  and  the  posi- 
tion of  the  Fermi  level.  The  results  of  this  analysis  support  the 

interpretation  of  the  noise  spectra  suggested  before.  For  = -0.1  V 
the  position  of  the  Fermi  level  is  in  between  the  first  and  second  quan- 
tization levels  (see  Fig.  2.16  where  this  is  indicated).  Either 

increasing  or  reducing  the  gate  voltage  will  position  the  Fermi  level 
closer  to  the  quantization  levels.  As  was  argued  before,  this  will  give 
rise  to  an  increase  in  the  g-r  noise,  hence  obscuring  the  1/f  noise. 

To  find  conclusive  evidence  of  the  interpretation  presented  here, 
more  noise  measurements  need  to  be  done.  Lowering  the  lattice  tempera- 
ture will  move  the  Fermi  level  higher  up  in  the  band,  and  higher  quanti- 
zation levels  can  be  probed. 

The  interpretation  presented  thus  far  considers  transitions  from 
one  subband  to  another.  Essential  is  that  the  mobility  of  the  electron 
changes  in  the  process,  i.e.,  the  mobility  in  the  different  subbands, 
should  be  different.  If  the  mobility  of  the  electron  does  not  change  in 
the  transition  process,  no  low-frequency  current  fluctuations  will 
beobserved.  When  the  Fermi  level  is  low  in  the  conduction,  in  the 

vicinity  of  the  first  quantization  level,  it  must  be  concluded  that  the 
only  sources  that  can  supply  electrons  that  transfer  to  the  first 
subband  are  discrete  levels  in  the  bandgap.  The  higher  subband  is  not 

occupied  yet  by  electrons;  hence  the  electrons  must  come  from  energy 
levels  below  the  first  subband. 
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Fig.  2.16.  The  relationship  between  the  gate  voltage  and 
the  position  of  the  Fermi  level  as  found  from 
the  analysis  in  the  text. 


GATEVOLTAGE  (V) 
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A completely  different  interpretation  of  the  measured  spectra  is 
in  terms  of  a resonance  in  the  1/f  noise.  The  magnitude  of  the  1/f 
noise  is  often  presented  in  terms  of  Hooge’s  parameter,  which  is  defined 
by  the  following  empirical  relation, 

si(f)  “h 

2 " In  (2.15) 

where  S^f)  is  the  spectral  density  of  the  current  fluctuations,  IQ  the 
dc  drain  current,  f the  frequency,  and  N the  number  of  carriers  involved 
in  the  noise  process.  The  number  of  carriers  is  estimated  from  R 

act  ’ 

the  resistance  of  the  active  region  under  the  gate  of  the  channel.  It 
is  assumed  that  all  of  the  excess  noise  is  produced  in  the  active  region 
and  that  the  access  resistances  only  produce  thermal  noise. 

In  Fig.  2.17  Si*f  is  plotted  versus  the  gate  voltage  for  two  dif- 
ferent drain  currents.  It  can  be  seen  that  the  noise  scales  quad- 
ratically  with  the  current,  which  is  expected  for  the  operation  of  the 
MODFET  in  the  ohmic  regime.  Accordingly,  the  parameter  a is  indepen- 
dent of  the  drain  current  within  the  experimental  error,  as  can  be  seen 
in  Fig.  2.17  where  aR  is  also  plotted  as  a function  of  the  gate  voltage. 
With  increasing  gate  voltage,  aR  seems  to  decrease  slightly,  increases 
slightly  around  VG  = -0.1  V,  and  decreases  again.  The  slight  increase 
might  indicate  a resonance  in  the  1/f  noise,  which  could  overshadow  the 
g-r  noise  produced  in  the  channel.  At  present  there  is  no  explanation 
for  this  phenomenon. 

In  this  section  two  possible  interpretations  were  given  of  the 
measured  noise  spectra.  From  the  measured  data  so  far  it  is  impossible 
to  decide  in  favor  of  either  one  of  the  explanations. 
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Fig.  2.17.  The  magnitude  of  the  1/C  noise  (S1  - f ) versus 
the  gate  voltage.  Also  Hooge's  parameter  a 
(for  definition  see  text)  is  shown  as  a H 
function  of  the  gate  voltage. 


CHAPTER  III 

LOW-FREQUENCY  NOISE  MEASUREMENTS  AS  A TOOL 
TO  ANALYZE  DEEP-LEVEL  IMPURITIES  IN  SEMICONDUCTOR  DEVICES 

3.1.  Introduction 

One  of  the  limiting  factors  of  the  performance  of  semiconductor 
devices  is  the  electrical  noise  produced  in  these  structures.  From 
resistors  to  transistors  and  from  silicon  to  modern  ternary  and  quater- 
nary compounds  [23,29],  much  attention  has  been  given  to  the  character- 
ization of  the  noise  sources.  Once  the  noise  sources  are  identified, 
strategies  can  be  devised  to  reduce  the  influence  of  the  electrical 
noise  on  the  device  performance.  To  give  a few  examples,  better  pro- 
cessing techniques  might  be  developed  or  ranges  of  operation  of  devices 
might  be  found  for  which  the  noise  is  minimal. 

In  this  chapter,  measurements  of  the  spectral  density  of  the 
voltage  fluctuations  of  planar  n-type  silicon  resistors  as  a function  of 
temperature  are  reported.  The  advantage  of  studying  noise  in  a planar 
device  is  that  a one-dimensional  treatment  is  warranted.  The  advantage 
of  using  a simple  resistor  is  that  the  relation  between  the  distributed 
noise  sources  and  the  power  spectral  density  of  the  voltage  fluctuations 
developed  across  the  terminals  is  relatively  simple. 

Impurities  in  the  semiconductor  materials  will  give  rise  to 
discrete  energy  levels  in  the  bandgap.  These  levels  are  capable  of 
trapping  free  charge  carriers.  This  is  called  capture  or  recombination. 
The  inverse  process  is  called  emission  or  generation.  Due  to  recombina- 
tion and  generation,  i.e.,  the  transition  of  a free  carrier  from  the 
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band  to  the  trap  level  and  vice  versa,  the  nunber  of  free  charge 

carriers  will  fluctuate.  These  number  fluctuations  give  rise  to  current 

or  voltage  fluctuations  whose  power  spectral  density  can  be  measured 

when  a dc  electric  field  is  applied  and  result  in  Lorentzian-shaped 
2 2 

( 1/ ( 1+0)  T ))  components  in  the  spectrum. 

A Lorentzian  is  characterized  by  its  zero-frequency  plateau  level 
and  its  characteristic  time.  This  time  constant  is  equal  to  the  inverse 
of  2tt  times  the  corner  frequency.  These  two  quantities  can  be  expressed 
in  the  pertinent  trap  parameters  that  describe  the  electrical  nature  of 
the  trap.  In  the  case  of  the  interaction  of  a single  trap  with  the  con- 
duction band,  analytical  expressions  of  the  noise  parameters  as  a func- 
tion of  the  trap  parameters  can  be  found.  Also,  the  inverse  process, 
the  extraction  of  the  trap  parameters  from  the  noise  measurements,  is 
straig’ntf orward  [30].  Although  for  two  traps  interacting  with  the  con- 
duction band  the  noise  can  still  be  described  analytically  in  terms  of 
the  trap  parameters,  the  inverse  process  is  virtually  impossible. 

Our  measurements  indicated  the  presence  of  at  least  two  traps.  In 
order  to  justify  the  extraction  of  the  trap  parameters  by  using  a model 
that  describes  the  interaction  of  a single  trap  with  the  conduction 
band,  we  studied  the  conditions  under  which  a decomposition  of  a system 
of  two  traps  interacting  with  the  conduction  band  into  two  systems  of 
one  trap  interacting  with  the  conduction  band  is  allowed.  A detailed 
discussion  of  the  analysis  and  validity  of  the  decomposition  is  pre- 
sented in  the  next  section.  A condition  is  found  for  which  the  decom- 
position is  allowed. 

In  section  3.3  the  extraction  method  of  the  trap  parameters  is 
presented.  A description  of  the  sample  that  was  used  and  the 
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application  of  the  analysis  technique  to  the  results  of  our  measurements 
are  presented  in  section  3.4.  Also,  with  the  help  of  an  exact  multi- 
level computer  program  based  on  the  theory  of  Van  Vlist  [5],  extracted 
trap  parameters  are  used  as  the  input,  zero-frequency  plateau  levels, 
and  characteristic  times  are  generated  as  a function  of  temperature. 
Since  the  presence  of  some  generation-recombination  (g-r)  noise  compo- 
nents could  not  satisfactorily  be  explained  with  two  or  three  single- 
donor systems,  a double-donor  model  was  developed,  and  expressions  for 
the  noise  were  derived.  This  model  is  presented  in  section  3.5. 

3.2.  Theory 

This  section  is  organized  into  four  subsections.  In  subsection 
3.2.1  we  develop  an  exact  one-trap-level  model  that  describes  the  zero- 
frequency  plateau  levels  and  the  characteristic  times.  The  following 
subsection,  3.2.2,  is  a generalization  for  the  case  that  two  trap  levels 
interact  independently  with  the  conduction  band.  In  subsection  3.2.3 
the  zero-frequency  plateau  levels  and  characteristic  times  for  the  exact 
two-trap-level  model  are  derived.  In  the  final  subsection  we  compare 

the  results  obtained  in  subsections  3.2.2  and  3.2.3  and  derive  a con- 
dition pertaining  to  the  validity  of  the  decomposition. 

3.2.1.  Exact  one— trap— level  model 

Generation-recombination  noise  is  recognized  in  noise  spectra  by 
the  shape  of  its  Lorentzian  frequency  dependence,  namely  1/(1  +U^T2) 

The  two  quantities  that  determine  a Lorentzian  are  its  zero-frequency 
plateau  value  S&n(f  = 0)  and  its  corner  frequency  f = J-  , where  x is 

C ZTTX 

a characteristic  time  related  to  the  trapping  process.  In  this  section 
we  will  calculate  these  two  characteristic  quantities  for  the  case  that 
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a single  trap  in  the  bandgap  interacts  with  the  conduction  band  (see 
Fig.  3.1). 

The  number  of  electrons  in  the  conduction  band  and  in  trap  level  i 
are  n^  and  n^  respectively.  The  generation  rate,  g^,  is  equal  to  the 
number  of  electrons  going  from  trap  level  i to  the  conduction  band  per 
unit  time.  The  reverse  transition  is  called  the  recombination  rate  r^ . 
In  the  bimolacular  notation  these  two  rates  are 


S 

r 


i 

i 


dioVo  ■ dioVao  - "o' 


dOtI'oi,i  d0l:10<”l  !'l) 


(3.1a) 


(3.1b) 


where  Pq  is  the  number  of  vacancies  in  the  conduction  band  and  is  the 
number  of  vacancies  in  the  trap  level.  denotes  the  density  of  states 
times  volume  and  Nj  is  the  total  number  of  trapping  centers.  i^  and 

di0  are  caPture  coefficients.  Assuming  space  charge  neutrality  and 
neglecting  the  interaction  with  the  valence  band,  we  have 


n0  Ni  “ ni  * (3.2) 

The  rate  equation  for  trap  level  i reads 

*0=h~  rL  ' 

We  develop  this  equation  around  thermal  equilibrium  concentration 
values,  indicated  by  a bar, 


n 

n 


0 

i 


+ An 
+ An 


0 

i 


51 


valence 

band 


Fig.  3.1.  Energy  band  diagram  of  a one-trap-level 

system.  Interaction  of  the  trap  with  the 
valence  band  is  neglected. 
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d 


dt  (no  + Ano}  “ h(a o + An0’nL  + Ani}  " ri(no  + An0’ni  + Ani} 


’i(aO’ni)  " ri(nO’ni)  +df^  " ri(n0’ni)]An0- 


(3.3) 


Because  of  equation  (3.2)  we  only  have  one  independent  variable:  n0- 

The  first  two  terms  of  (3.3)  cancel  because  of  thermal  equilibrium. 


"rz"  n_  = 0,  equation  (3.3)  becomes 
at  0 


dT  An0  “ ['di0p0  ' diOHi  ' -Vi  - d0iKn]Ari3  (3.4a) 


Ti  An0 


(3.4b) 


Here  we  used  equations  (3.1)  and  (3.2).  The  condition  of  detailed 
balance  relates  g^  to  r^>  We  obtain 

d , „n  p , = d np 
iO  r 0 0i  0 l 


d , Pi 

di0  d0i  - - 

P0  ni 


(3.5) 


Now  we  insert  eq.  (3.5)  in  eq.  (3.4a)  and  we  find  for  An„ 


--  nn  P 


0 pi  - n0  Pi 


A"0=-d0i[p0^  ^ + ^ +Pi+n0]An0 


P0  ni 


po  ni 


= +-r-1  + p^1  +r^)]% 


(3.6) 
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with 


Ti1  = +-rL)  + p,(i  + ^)]< 


oi  * 


(3.7) 


To  calculate  the  noise  spectrum  we  add  the  Langevin  terms  to  equatLon 
(3.4b).  Hence 


An, 


dt  Ano  = 


+ F - F 

T-  *1  rl 


i 


(3.3) 


We  now  perform  a Fourier  transformation  on  this  equation  (3.8)  and 
obtain 


(J»  + - Fgi  - 5ri 


(3.9) 


where  the  tilde  denotes  the  Fourier  transform.  3y  definition 


SAn0(u) 


lim  f <AnQAn*> 

1 -fco 


- (SAg  (u)  + SAr.  (u)) 


i 


Ar, 


,.22 

1 + 0)  T 


where 


V = \ (f> 

8i 


and 


(3.10) 


V = SF  (f)  • 

1 ri 

The  angular  brackets  <C  > denote  the  ensemble  average, 
cussions  we  rewrite  (3.10)  as 

t2 

S.  (a.)  = <F  F*> 

An0  i i‘  1 + u2t? 


For  later  dis- 


(3.11) 
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with 

<FiV  =  1 * * *  SASl(a))  + SAr  <u)  ' (3-12) 

3-2.2.  Generalization  to  two  independent  trap  levels 
Applying  the  nodel  of  section  3.2.1  to  two  traps  interacting 
independently  with  the  conduction  band  results  in 


and 


(3.13a) 


(3.13b) 


as  expressions  for  the  corner  frequencies,  if  we  number  the  traps  1 and 
2 (see  eq.  (3.7)  and  Fig.  3.2).  The  noise  spectrum  will  look  like: 


<¥t> 


1 + U T 


TT  + <^2  V 


1 


. , 2 2 
1 + u t2 


which  follows  from  eq.  (3.11). 

3-2.3.  Exact  two-trap  level  model 
In  this  section  we  calculate  the  characteristic 
values  of  an  exact  two-level  trap  system.  Space  charg< 


(3.14) 


times  and  plateau 
neutrality  gives 


n0  ~ (N1  " ni>  + (N2  " n2)  * (3-15) 

In  this  model  we  have  two  independent  variables  tig.nj  because  n2  follows 
from  eq.  (3.15).  The  rate  equations  for  the  two  independent  variables 


are 
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VALENCE^ 
BAND 


Fig.  3.2.  Energy  band  diagram  of  a two-t rap-level  system.  Inter- 
action of  the  traps  with  the  valence  band  is  neglected. 
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n0  = S2 (nQ  »n2^  ~ r2^n9>n2^  + Si^no»ni)  ” rl^nO’nl^  (3.16a) 

nl  = “®i^no,n1^  + rl(aO’tll)  • (3.16b) 

■Ve  neglect  the  probability  of  a direct  translation  from  trap  1 to  trap  2 
or  vice  versa.  Again  we  develop  around  thermal  equilibrium 


no  " no  + An0’  ni  = ni  + Ani 
and  find  from  eq.  (3.16a) 


An0  S2(n0’n2)  r2(n0’n2)  +d ~ -Mno’n?)  ” r2*an»n?^  Anr 

0 r 


+ dir  [g2(n0’n2)  ■ r2(Vn2)]  Anl  + - Fr 

1 r 6 9 ^ 


+ g1(n0»n]L)  r1(nQ,n1)  +—  [gj^.n^  - r1(nQ,n1)]  Ann 

0 t 

+ dT[gl(n0-nl)'rl(n0’nl)lAnl+?J  * F • <3. 17a) 

1 t °1  rl 


Similarly,  we  find  from  eq.  (3.16b) 


Atll  ~Zl(nC),nl)  + r1(nQ,n1)  - [g^n^n^  - r^n^n^]  AnQ 

" 1^7  [Si(VV  " rl(n0’nl)]  Atli  - F + ?r  (3.17b) 

1 t S1  rl 

where  the  subscript  t stands  for  taking  the  derivatives  In  thermal 
equilibrium.  Because  of  detailed  balance  the  underlined 


terms  in  eqs. 


(3.17a)  and  (3.17b)  cancel.  We  write  eqs.  (3.17a)  and  (3.17b)  as 
follows : 


AnQ  = -Mn  AnQ  - A^  + Fj  + F£  (3.18a) 

aAx  = -M21  An0  - M22  Anx  - Fj  (3.18b) 


with 


fi  , 

M11  = ~d^7  ig2(n0’n2)  " r2(n0,tl9)  +Sl(n0^1)  - r (n  .n^] 

0 t 

(3.19a) 

M12  dn.  "g2^n0,n2^  “ r2^no,n2^  + S1 ’n i ^ “ ri^no’ni^i 

1 t 

(3.19b) 


M21  = dn7  ^lK’V  " rl(n0’tll)J  (3-19c> 

0 t 

M22  = d5“  ^gl (n0 ,nl ) “ rl(n0’nl)]  (3.1 9d ) 

1 t 

F = F - F 
1 gl  rl 


We  now  use  the  phenomenological  rates  and  eq.  (3.15)  to  eliminate  n2  and 
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ri  ' d<nVi  * doino(!,i  ‘ V 


(3.20b) 


g2  = d20n2P0  = d20(Nl  + N2  " n0  ' nl)(N0  ' V (3’20c) 


r2  d02n0P2  d02n0(n0  + nl  " • 


(3.20d) 


We  substitute  eqs.  (3.20a)  - (3.20d)  into  eqs.  (3.19a)  - (3.19d)  and 
obtain  for  the  phenomenological  relaxation  matrix  M: 

M11  = d20p0  + d02n0  + d20n2  + d02P2  + d10nl  + dQlpl  (3.21a) 


M12  d20p0  + d02n0  " d10p0  " d01n0 


(3.21b) 


M21  d10nl  d01pl 


(3.21c) 


M22  d10P0  + d01n0  • 


(3. 2 Id) 


We  use  detailed  balance  to  express  d1Q  in  terms  of  dQ1  and  rewrite 
expressions  (3.21a)  - (3.21d): 


11 


" d02n0^  +~)  + d02p2^ 1 + ^)  + d01pl^  (3.22a) 


d02n0*'1  + - ) d01n0^1  + 

n. 


12 


(3.22b) 


21 


= -doi?^1  + zr) 


(3.22c) 
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The  eigenvalue  equation  of  matrix  M reads 

X2  - <*11  + *22>X  + "u"22  - "12*21  ’ <X  - V<X  - 12)  ' 0 ' 

(3.23) 

The  reciprocal  eigenvalues  X^  are  the  characteristic  times  of  the 
noise  spectra  as  we  will  show  later  (see  eq.  (3.31)). 

For  Xt  we  find  from  eq.  (3.23) 

Al,2  = 1/2  Tr(M)  * l'2  (Tr2Cl)  - 4 Det(M))l72  . (3.24) 

We  now  calculate  the  noise  spectrum  from  equations  (3.18a)  and  (3.13b). 

rourier  transforming  these  equations  and  using  vector  notation,  we 
obtain 


+ Mu 


12 


An„ 


M 


21 


jm  + M 


22 


An, 


F + F 
1 2 


(3.25) 


where  the  tilde  denotes  the  Fourier  transform  of  the  quantities. 

We  will  show  that  this  matrix  is  nonsingular,  later.  We  can  now 
invert  the  matrix  and  solve  for  An^  and  An^ 


Ano 

1 

J“  + *22  - *12 

?1  +f2 

A“l 

Det 

*21  J“+"ll 

-?1 

(3.26) 

where  Det  w + M^^22  “ ^12^21  + + ^22^  = t*le  determinant  of 

the  matrix  in  equation  (3.26).  From  eq.  (3.26)  we  calculate  AnQ  : 
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i*0  ■ 0^  t(J“  - !,22K?1  + F2>  * ' 


(3.27) 


and  for  S,  (w)  : 

4no 


hn0M  ■ <15F  [<3“  + n22)(Fl  + V + MI2fl) 

5SS  [<-*•  + r,22K?l  + F2>  + 


= < 


| Det | ' 


r , z . o . ~ ~ ~w  ~ ~ ~-,v  9 ~ ^,v 

[(“  + * 22  F1F1  + F1F2  + F2F1  + F2F2)  + MI2  F1F1 


+ (jo.  + m22)m12(f1f1  + f2f^)  + (-jo.  + m22)m12(f1f^  + f1f*)]> 


The  mixed  Langevin  terms  cancel  if  we  calculate  the  average  value  and 
assume  that  there  is  no  correlation  between  the  Langevin  noise  sources. 
Hence 


SAn  ^ = 

A 0 I Det 


1 


[(m2  + 


(M 


12 


+H22>2) 


<FiFi>  + 


, 2 
(co  + 


M2  ) 

22 1 


<F2F2>] 


(3.28) 

with 


lDet 1 2 - I’"2  + Mlia22  - "U'V2  + “2<hl  + M22>2 


[oj2  - Det(M)]2  + 0)2Tr2(M)  . 


(3.29) 


The  right-hand  side  of  equation  (3.29)  can  be  written  as  the  product  of 
two  factors,  see  ref.  [4].  Comparing  eqs.  (3.23)  and  (3.29),  we  have 
the  r.h.s.  of  eq.  (3.29) 


2 2 

- Det  (M)  + joiTr(M)  ] [to  - Det(M)  - jioTr(M)  ] 
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' '°2  - "ll^J  +M12M21  + J“(HUf,22)1 

» [«2  - «n«22  + H12M2I  - j»(in  * „22 

= (j«  - ^jKju  - X2)(ju  + X1)(ju  + X ) 

= (u2  + A^)(u2  + X“) 


)] 


(3.30) 


where  are  the  eigenvalues  of  matrix  M.  We  can  now  spilt  up  the 

right-hand  side  of  equation  (3.28)  Into  two  terms,  two  Lorentzlans. 


2 

10 


+ 


2 

co 


where 


(3.31) 


a = 


[(hj2  + :i22) 


h2]  <hh>  + I«22 


2 ~ -w-A* 

V <F2F2> 


X2  - 


(3.32a) 


[X2  ' (*12  + M22)2^  <FlV  + [X2  ‘ fl22]  <?2F2> 


b = 


2 2 
X2  - A* 


(3.32b) 


It  is  clear  from  eq.  (3.31)  that  the  eigenvalues  A^^  can  be  Identi- 
fied with  the  measured  characteristic  times. 

3*2.4.  Comparison  of  the  models  described  In  sections  3.2.2 
and  3.2.3 

In  this  section  we  first  derive  a condition  under  which  the  char- 
acteristic times  of  the  exact  two-trap  level  noise  model  (section  3.2.3) 
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are  adequately  described  by  the  parameters  of  the  generalized  two-trap 
level  model  (section  3.2.2),  i.e.,  a”1  = x 

From  comparing  expressions  (3.13a)  and  (3.13b)  with  equations 
(3.22a)  and  (3.22b),  we  find  that  we  can  identify: 


T1  M22  M21 


T2  M11  + M21  • 


(3.33a) 

(3.33b) 


-1 


If  we  assume  A^,  - to  be  valid,  then  we  can  approximate  the 

eigenvalue  equation  for  the  matrix  M (eq.  3.23)  by 


<T"‘  - M22  +r,21><T_1  * "n  * V ‘ 0 


(3.34) 


or 


<t'V  - <MU  - M22)T->  + MuM22  - MnM21  + M21H22  - M2|  . 0 . 


(3.35) 


But  the  matrix  M has  the  form 


M11  M12 


M21  M22 


which  gives  the  following  eigenvalue  equation: 


-1, 


<T  P*  - (Mu  + M22)t_1  + MjjMjj  - M12M21  - 0 


(3.23) 


We  write  this  equation  (3.23)  in  the  form  of  equation  (3.35)  with  an 
added  term: 


<T-  >*  - <MU  + HJ2),-1  + MuM22  - «n«21  + M21M22  - M21 


• "l2M21  + M11M21  - M21M22  + M2l'  0 


(3.36) 
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In  order  to  approximate  the  exact  eigenvalues  by  t"1  and  t 
elude  that  we  must  have 


we  coq- 


I Mi2lM21 


MllM2l“  M21M22  + M2ll 


« |m 


M 

11  22 


- muj21  + m21m22  - m21  | 

(3.37) 


or 


+ M. 


11 


- m22  + m91)|  « 


21' 


| M M 
1 11  22 


M21(-Mll 


+ M_ 


22 


~ ^21 } I 


(3.33) 


Using  equations  (3.22a)  - (3.22d),  the  left-hand  side  of  eq.  (3.33) 
becomes 


01  0 


- r Pi.  _ Pi 

“nt1  + — ) + dn2nnf1  +-) 


L1 


doi"oC'  + ^ 

“l 


doiPiC1  +^)ll  - <11dffiS1i!(1  +^)2  . 


(3.39) 


The  right-hand  side  becomes 


I [d02n0^ 1 + - ) + d02p2^ 1 + ^ + doiPl^  + T~)  ] doinO^ 1 + 3“ 

n, 


n0w  . - . P2,  . - . V 


“ ^iPi^1  + _ ) [-  d02n0^  +~)  - dQ2p2^ 1 +~)  - doiPiU  +r^) 
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+ d0 1 ^0  ^ 1 +^) 

nl 


n0,  . p2' 


+ nQ p x ( 1 + — ) (l  + ^) 


(3.40) 


Inequality  (3.38)  becomes 


"0,2  2 


plp20  + - ) ^ n0O  + C1  + + nnP7(l 

n. 


0 2^ 


+ ViC1  C1  +~)  + P^U  +Z^)2  • (3.41) 

P0  n2  Po 

For  a nondegenerate  (n-type)  semiconductor  the  Fermi  level  is  well  below 

the  conduction  band.  The  number  of  unoccupied  states  in  the  conduction 

band  is  far  greater  _than  the  number  of  electrons  in  the  conduction  band 

n0 

and  the  factor  1 + — = 1 . Hence  eq.  (3.41)  becomes 

P0 


6 5 


plp2  « nQ(l  + ^)  (1  + — ) + nQP2(l  + —)  + n p ( 1 + ~)  + p p . 
n2  nl  1 ni  ' n2 

(3.42) 

Ue  will  consider  this  condition  for  the  following  three  cases:  1)  the 


Ferrai  level  Is  between  the  conduction  band  and  the  upper  one  of  the  two 
trap  levels;  2)  the  Ferrai  level  Is  below  the  lowest  of  the  two  trap 
levels;  3)  the  Fermi  level  Is  between  the  two  trap  levels  (see 
Fig.  3.3). 

1)  With  the  Ferrai  level  well  above  (more  than  4 k^T)  both  traps, 
the  traps  are  completely  filled. 


nl  = N1  ?1  <<n! 

n2  = N2  P2  < n2 

The  right-hand  side  of  inequality  (3.39)  is  then 

_2 

= nQ  + nQp2  + nQp L , 

which  is  certainly  much  larger  than  . 

2)  For  the  case  in  which  the  Ferrai  level  is  below  the  lowest  trap 
level,  both  traps  are  empty,  that  is, 


P1  >nl  P1  = Nl 

p2  * n2  p2  ~ N2  * 


The  right-hand  side  of  inequality  (3.42)  becomes 


'1 


1 


-2  P2  . 

n0  ~ ~ + V2  ~ + n0Pl  ~ + plp2 

n„  n,  n n 

2 1 1 2 
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CONDUCTIO  N 
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trap  1 


E 


trap  2 


Fig.  3.3.  Energy  band  diagram  of  a two-trap-level  system, 

showing  the  three  Fermi-level  positions  discussed 
in  the  text. 
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Divide  by  p^p^  » and  we  get 
2 

.4  ♦ i i ♦ 1 , 

"ln2  °1  "2 

which  is  certainly  larger  than  1. 

3)  When  the  Fermi  level  lies  between  the  two  trap  levels,  then 
trap  2 is  completely  filled  and  trap  1 is  empty. 

N2  = n2  > p2  ’ nl  € pi  = Ni  • 

Substitution  into  the  right-side  of  (3.42)  results  in 


-2  P1  , - - P1 


‘0  - + nop2  ~ + plp2  > p ip 9 

nl  nl 


If  we  divide  by  px,  then  we  find  that 

2 


0 + n0p2 
i n. 


+ a0  + P2  > P2  • 


So  for  these  three  cases  the  condition  expressed  by  equation  (3.42)  is 
always  fulfilled  if  the  Fermi  level  is  more  than  4 k^T  away  from  the 
trap  level. 

We  will  now  consider  the  noise  plateau  values.  Compare  eq.  (3.31) 
with  eq.  (3.14),  which  describe  the  noise  of  the  exact  two-level  model 
and  the  noise  of  two  independent  trap  levels,  respectively.  To  identify 
both  expressions  with  ench  other  we  expect • 


X1  H22 


*2  ' M22  + M12 


In  order  that 


(3.43) 


be  valid,  the  two  following  expressions  have  to 
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hold: 


'V  «l"22l 


(3.44a) 


and 


M11  + M21  M22  + M12  ' 


(3.44b) 


From  equation  (3.44a)  we  find,  using  equations  (3.22a)  - (3.22d) 


+r^)  <d0in0(1  +^)  . 


(3.45a) 


Using  again  equations  (3.22a)  - (3.22d),  equation  (3.44b)  leads  to 


d02n0^  + - ^ + d02P2^  + " d02n0^  + ^ 
n2  po  n: 


from  which  we  conclude: 


P2(!  + ~ ) < n0(l  + . 


(3.45b) 


Rewriting  equations  (3.45a)  and  (3.45b)  gives 


*r  + + i- 

no  po  ni  pi 


i_ 

n. 


♦ i-  < + 


(3.46) 


For  a nondegenerate  n-type  semiconductor  it  holds  that  nQ  < pQ  . If  in 
addition  nQ  > Nj  and  nQ  > N2  , then  Inequalities  (3.46)  are  fulfilled. 
In  conclusion,  it  is  shown  here  that  under  the  condition 
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Pl?2  < nQ(l  +—)(!  + 


P1  _ Pi  _ p 

— ) + nQP2(  1 + — ) + nnPi(l  + + PiP 

nl  nl 


n„ 


the  exact  two-trap  model  can  be  approximated  by  a model  In  which  each 
trap  interacts  Independently  with  the  conduction  band.  The  condition  Is 
satisfied  if  i)  the  Fermi  level  is  few  k^T  away  from  the  levels  con- 
sidered; ii)  nQ  > and  nQ  > N2  . 

3.3.  Extraction  of  Trap  Parameters 

In  analogy  with  Copeland's  approach  we  consider  the  interaction  of 
one  impurity  trap  level  with  the  conduction  band.  From  expressions 
desclbing  the  measured  g-r  noise  spectra  the  four  trap  parameters,  i.e. , 
energy  level,  trap  density,  spin  degeneracy,  and  capture  coefficient, 
are  extracted.  These  trap  parameters  completely  describe  the  electrical 
properties  of  a trap.  Although  our  treatment  is  analogous  to  the  one 
followed  in  ref.  [3],  we  use  a notation  more  commonly  used  in  g-r  noise 
description.  In  thermal  equilibrium  the  generation  rate  equals  the 
recombination  rate,  hence 


g = r . 

Copeland  formulates  expression  (3.47a)  as 


(3.47a) 


VTe  ‘ (t,T  ‘ nT)/Tc 


(3.47b) 


where  rip  is  the  number  of  electron-occupied  impurity  levels,  tip  the 
total  number  of  impurities,  a time  constant  associated  with  the 
generation  process  (emission),  and  ic  a time  constant  associated  with 
the  recombination  process  (capture).  (See  Figure  3.4.)  A small  per- 
turbation, AnT,  of  the  equilibrium  value  of  nT  relaxes  back  through  the 
mechanisms  of  generation  and  recombination: 
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valence 

band 


Fij.  3.4.  Interaction  of  a trap  with  the  conduction  band. 

g describes  the  generation  process  of  an  electron 
from  the  impurity  level  to  the  conduction  band  and 
r describes  the  inverse  process. 
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dt  Anr  = "Ag 


+ Ar 


AaT  A(NT-nT) 


T T 

e c 


T 


(3.43) 


where 


l 

T 


(3.49) 


Copeland  claims  that  the  measured  corner  frequency  of  a Lorentzian- 
shaped  spectrum  equals  1/(2ttt)  where  t is  given  by  eq.  (3.49).  It  is 
shown  here  that  this  is  an  approximate  expression  whose  validity  depends 
on  the  ratio  of  the  number  of  deep  traps  (N„)  and  the  number  of  shallow 
intended  dopants  (N^) . In  the  g-r  noise  formalism 


an,n„„ 
T SR 


[instead  of  n^/t  ] 


r = 


an(NT~nT) 


N - n 

[instead  of  ] 

c 


(3.50) 

(3.51) 


and 


1/T  = " h (s_r)t  • <3*52> 

In  these  expressions  n is  the  number  of  electrons  in  the  conduction 
band,  a is  a proportionality  factor  associated  with  the  recombination 
process  which  can  be  written  as  ( 5 ’v^/volume)  where  5 is  the  capture 
cross  section  and  vth  is  the  thermal  velocity.  nSR  is  the  Shockley-Read 
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quantity  which  is  aqual  to  the  number  of  electrons  in  the  conduction 
band  if  the  Fermi  level  were  to  coincide  with  the  trap  level  con- 
sidered. The  subscript  t denotes  that  the  derivative  is  calculated  in 
the  thermal  equilibrium  condition.  3y  using  the  condition  of  space- 
charge  neutrality,  one  of  the  variables,  n^,,  can  be  eliminated  from  eqs. 
(3.50)  and  (3.51),  and  eqn.  (3.52)  yields 


This  formalism  gives  an  extra  term,  1/x  , in  the  expression  for  the 

measured  time  constant,  t . It  can  be  neglected,  to  give  Copeland's 
result,  if 


(3.53) 


where 


(3.54) 


(3.55) 


Equating  eqs.  (3.50)  and  (3.51)  in  thermal  equilibrium  gives 


N, 


T 


n, 


T 


n 


SR 


(3.56) 


Expressing  Oj.  tn  terms  of  Fermi-Dirac  statistics: 


-1 


(3.57) 
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where  3 is  the  spin  degeneracy,  EL^,  is  the  energy  level  of  the  impurity, 
kg  is  Boltzmann's  constant,  and  T is  the  temperature.  Under  the 
assumption  that  the  number  of  carriers  in  the  conduction  band  is  given 
by  the  number  of  shallow  dopants  (all  donors  are  ionized)  and  using  eqs. 
(3.56)  and  (3.57),  eq.  (3.55)  will  read 

4 cosh  (y/2)  » Nt/Nd  . (3.58) 

Equation  (3.58)  holds  for  all  y if  Nt/Nd  « 4.  Under  this  condition  it 
is  appropriate  to  approximate  1/x  by  a(n  + n).  With  the  help  of  eqs. 
(3.56)  and  (3.57)  we  write 

1/t  = ct(n  + nSR) 

et  " ef 

= an[B  exp ( + l]  = an [ e xp ( y ) + 1]  . (3.59) 

B 

From  this  equation  the  general  temperature  dependence  of  t can  be 
found.  In  the  low-temperature  limit,  when  the  Fermi  level  lies  above 
the  trap  level,  the  exponential  in  eq.(3.59)  can  be  neglected  with 
respect  to  1 and  the  characteristic  time  is  virtually  a constant.  As 
the  temperature  increases,  specifically  when  the  position  of  the  Fermi 
level  is  below  that  of  the  trap,  the  exponential  will  dominate  the 

r.h.s.  of  eq.  (3.59)  and  a strong  temperature  dependence  occurs.  This 
is  depicted  in  Fig.  3.5. 

The  procedure  to  extract  the  trap  parameters,  according  to 
Copeland  [4],  goes  as  follows.  The  expression  for  the  spectral  density 
of  the  number  fluctuations  of  the  number  of  electrons  in  the  conduction 
band  is 

sn(“)  = 4gr2(l  + “2t2)  1 = Sn(0)(l  + u2t2)_1  . (3.60) 


7 4 


Fig.  3.5.  Qualitative  temperature  dependence  of  the 
characteristic  time. 
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The  zero-frequency  plateau  value  Is,  using  eqs.  (3.50)  and  (3.59) 


S(0)  = 4gx2  and  S (0)/x  = 4 
n n n + n 


(3.61) 


SR 


From  detailed  balance  n^  can  be  expressed  in  n,  and  n^.  Hence  eq. 

(3.61)  can  be  written  as 


V°> 


= 4 N„ 


SR 


(n  + 


nSR^ 


(3.62) 


If  we  now  use  eqs.  (3.56)  and  (3.57),  this  equation  transforms  into 


Sn(0) 

— - — = 4 NT[exp(y/2)  + exp(-  y/2)]--  . (3.63) 


Clearly,  if  the  r.h.s.  of  eq.  (3.63)  is  plotted  semilogarithmically  as  a 

o 

function  of  y,  an  inverted  cosh  dependence  is  observed.  It  is  from  the 
slopes  of  this  curve,  away  from  the  center,  that  the  energy  level  of  the 
trap  is  calculated.  But  the  parameter  y still  contains  the  temperature 
dependence  of  the  Fermi  level.  To  eliminate  this  we  make  use  of 
resistance-temperature  measurements  we  performed.  The  resistance  is 
given  by 


R 


= L/Ao  = L~/qpn 


(3.64) 


where  a is  the  conductivity,  L the  length  and  A the  area  of  the  device, 
q the  magnitude  of  one  electron  charge,  and  y the  electron  mobility. 
The  temperature  dependence  of  the  electron  mobility  in  low-doped  Si  is 
accurately  known  [31],  From  the  Boltzmann  expression  for  n and  eq. 
(3.64),  Ec  - Ep  can  be  expressed  in  terms  of  the  measured  resistance  and 
the  density  of  states.  In  the  low-temperature  limit  (y  < 0) , this 
procedure  results  in  the  following  expression 
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G 


S (0) 
n 

T 


L 

qpRN^A 


4N^B  exp 


(3.65a) 


Now  it  is  very  clear  that  from  a semilogarithmic  plot  of  G versus  1/T 
the  slope  will  yield  the  energy  level  of  the  trap.  Alternatively,  in 
the  high-temperature  limit  (y  > 0) 

S (0)  quRN„A  E - E 

G = — = 4NT3  exP( iTy^-)  • (3.65b) 

3 

Depending  on  which  side  of  the  curves  the  slope  can  be  determined  best, 
either  one  of  eqs.  (3.65a)  and  (3.65b)  can  be  chosen.  Another  feature 
of  expression  (3.63)  is  that  for  the  peak  value  it  holds  that  the 
parameter  Y = 0.  From  eq.  (3.59)  this  condition  gives  the  following 
relationship  between  the  energy  of  the  trap  level  and  the  spin  degen- 
eracy: 

n(V  Ec  - Et 

B ‘ HC(T0)AL  > 0.66) 

where  Tg  is  the  temperature  for  which  this  maximum  occurs. 

The  last  feature  of  eq.  (3.63)  that  needs  to  be  discussed  is  the 
value  of  the  maximum.  This  value  equals  the  number  of  impurities. 
Thus,  from  plotting  eq.  (3.63)  and  either  one  of  eqs.  (3.65a)  and 
(3.65b),  three  trap  parameters  can  be  resolved.  The  last  parameter  to 
be  extracted  is  the  capture  coefficient  a.  From  eq.  (3.62)  this 
quantity  is  found  to  be 


a = 4NX 


n • n 


SR 


3 ' 


Sn(0)  (n  + nSR) 


(3.67) 
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3.4.  Experimental  Results  and  Analysis 

The  experimental  setup  is  described  In  ref.  [28].  The  samples  are 
planar  n+-n-n+  silicon  resistors.  The  length  of  the  n-reglon  is  78  pm, 
whereas  the  surface  area  is  1 x 1 mm2.  The  silicon  is  intentionally 
doped  with  phosphorus  resulting  in  a resistivity  of  16  «o  at  room 
temperature  for  the  n-region.  Aluminum  is  evaporated  on  the  n+-contact 
surfaces  to  provide  ohmic  contacts.  In  Fig.  3.6  the  experimental  values 
of  the  resistance  are  presented  as  a function  of  the  temperature. 
Starting  from  the  high-temperature  side,  the  curve  follows  the  inverse 
temperature  dependence  of  the  mobility  for  nearly  intrinsic  Si  [31]. 
For  temperatures  below  250  K the  temperature  dependence  of  the  concen- 
tration of  mobile  carriers  starts  to  contribute  to  the  temperature 
dependence  of  the  resistance. 

In  Fig.  3.7  we  show  the  electron  density  as  a function  of  the 
inverse  temperature.  The  plotted  data  are  calculated  from  the  resis- 
tance measurements  by  using  the  mobility  measurements  of  Jacoboni  at  al. 
[31].  Starting  from  the  high-temperature  side,  we  find  an  initially 
constant  electron  density.  For  temperatures  around  250  K we  observe  the 
filling  of  a deep  trap  (0.30  eV)  resulting  in  the  decrease  of  the  free- 
electron  density.  When  the  temperature  is  decreased  even  further,  the 
filling  of  a shallower  (0.16  eV)  trap  is  observed.  In  the  latter  case 

the  electron  density  seems  to  decrease  by  more  than  the  density 
. 17  -3 

U.5  x io  m , as  we  will  later  find)  of  the  shallower  (0.16  eV) 
trap.  This  can  be  attributed  to  the  fact  that  the  temperature 
dependence  of  the  mobility  is  not  so  accurately  known  in  this  tempera- 
ture range,  which  results  in  an  error  in  the  calculated  values  of  the 
mobile  carrier  concentration. 
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Fig.  3.6. 


Measured  resistance  values  as  a function 
of  the  temperature. 


electron  density  (m“3) 
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Fig.  3.7.  The  density  of  mobile  electrons  in  the 
conduction  band  calculated  from  the 
resistance-temperature  measurement . 
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Noise  measurements  performed  on  a homogeneous  semiconductor  cannot 
distinguish  between  accept or- like  behavior  and  donor— like  behavior  of  a 
trap.  In  the  case  of  two  independent  trap  levels,  there  are  four 
possible  configurations.  Each  of  these  possibilities  results  in  a 
different  value  for  the  shallow-donor  concentration,  as  calculated  from 
the  resistance  value  at  high  temperatures.  In  the  case  that  both  traps 
are  donors,  the  high-temperature  value  of  the  electron  density  would 

equal  the  sum  of  the  shallow  donor  concentration  (N_)  and  the  concen- 

u 

trations  (N,^ , N,^)  of  the  traps.  On  the  other  hand,  if  both  traps  were 

acceptors,  then  the  electron  density  on  the  high-temperature  side  of  the 

curve  would  just  equal  the  shallow  donor  concentration.  The  former  case 

would  result  in  Np  * 2.0  x lO^m'3  while  the  latter  would  give 
18  ~3 

ND  = 2.3  X 10  rn  . The  description  of  the  g-r  noise  is  the  same  in 

both  cases  as  it  is  identical  to  the  other  two  cases  where  one  trap  is 

an  acceptor  and  the  other  a donor. 

The  measured  zero-frequency  plateau  values  of  the  spectral  density 

of  the  voltage  fluctuations  are  shown  in  Fig.  3.8  and  the  corresponding 

characteristic  times  are  shown  in  Fig.  3.9.  From  an  inspection  of 

Fig.  3.5  the  data  can  be  divided  into  two  groups,  each  corresponding  to 

S (0)  2 

one  peak.  A semilog  plot  of  — versus  1/T  for  each  of  the  two 

V2  T 

data  sets,  as  we  have  done  in  Fig.  3.10,  has  the  number  of  considered 
Impurities  as  a peak  value.  Note  that  for  an  ohmic  device 

S (<*»)  S (w) 
n V 


where  is  the  spectral  density  of  the  voltage  fluctuations  and  V 

the  DC  voltage  applied.  For  the  high-temperature  noise  peak  the  density 
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Fig.  3.8.  The  normalized  zero  frequency  noise  plateau  levels 
versus  the  Inverse  temperature.  The  measured  data 
are  Indicated  by  ■ and  the  simulated  data  by  the 
full  line. 
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Fig.  3.9.  The  characteristic  times  belonging  to  the  zero- 
frequency  noise  plateau  levels.  The  measured 
data  are  indicated  by  ■ and  the  simulated  data 
by  the  full  line. 
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Fig.  3.10.  The  Copeland  function  (see  text)  S (0)/x  as  a 
function  of  the  inverse  temperature.  The  peak 
values  determine  the  impurity  densities.  The  full 
lines  are  drawn  to  guide  the  eye.  The  dashed  line 
tentatively  describes  a third  trap.  This  is 
discussed  in  more  detail  in  the  text. 
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17  3 

is  found  to  be  6 x io  /m  , the  temperature  at  which  the  peak  occurs 

being  Tq  = 235  K.  Most  of  the  measured  data  for  this  peak  lie  on  the 

low-temperature  side.  Accordingly,  we  choose  to  plot  for  these 
. , Sv(0)  5 

data  (see  previous  section)  — versus  1/T.  The 

TV  (quR)JN  A 

results  are  shown  in  Fig.  3.11.  From  the  slope  we  find  that  Eq  - Up  = 
0.30  eV.  The  spin  degeneracy  is  computed  using  eq.  (3.66)  and  found  to 
be  equal  to  0.3. 

The  last  parameter  to  be  extracted  is  the  capture  coefficient 
(a)  or  the  capture  cross  section  (6).  The  latter  quancity  is  found  from 
plotting  eq.  (3.67)  but  taking  into  account  the  temperature  dependence 
of  the  thermal  velocity.  We  show  these  results  in  Fig.  3.12. 


Analogously  the  parameters  of  the  second  trap  level  are  found.  In  this 

case  most  of  the  data  are  on  the  high-temperature  side  of  the  peak.  To 

S (0)  l3n  A 

obtain  the  energy  difference  Ec  - , we  now  plotted  -L_ £_  versus 

1/T.  In  Table  3.1  the  values  found  "for  the  trap  parameters  for  both 
traps  are  listed. 

To  make  the  analysis  complete,  the  two  assumptions  made  in  the 

previous  section  and  the  conditions  derived  in  sec.  3.2  pertaining  to 

the  validity  of  the  decomposition  of  the  measured  noise  spectra  into  two 

Lorentzians  need  to  be  checked.  The  condition  N^/N^  « 4 is  clearly 

fulfilled  for  both  trap  densities  and  for  either  combination  of  donor- 

like  traps  or  acceptor-like  traps. 

From  the  resistance  measurements  the  position  of  the  Fermi  level 

is  obtained  as  a function  of  temperature.  The  minimum  value  of  E_  - E„ 

C F 

is  found  to  be  0.1  eV  at  T = 78  K.  Clearly  the  Fermi  level  is  more  than 
4kgT  below  the  conduction  band,  and  the  assumption  of  a Boltzmann 
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The  modified  Copeland  function  versus  the 
inverse  temperature  for  one  of  the  traps. 
From  the  slope  the  ionization  energy  of  the 
trap  level  can  be  calculated. 


Fig.  3.11. 
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Fig.  3.12.  The  capture  cross  section  as  a function 
of  the  temperature,  for  one  of  the  traps. 
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Table  3.1 

Values  of  the  Trap  Parameters 


Trap  1 

Trap  2 

Unit 

Density 

6*1017 

1. 5*1017 

nf3 

Activation  Energy 

0.30 

0.15 

eV 

Spin  Degeneracy 

0.3 

1.5 

— 

Capture  Cross 
Section 

3x10-14xT-2-3 

2xlo-21 
1 . 5x10-27xT2 

T < 115  K 
T > 115  K 

O 
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distribution  of  the  free  elect rondensity  is  valid.  The  last  condition 
to  be  fulfilled  is  eq.  (3.42).  This  condition  reads: 


where 


Pi  _ “ tU,  - 'L  ~ n , number  of  empty  trap  1 states, 

1 1 1 1 

and 

P9  = nt  _ n = N ~ n , number  of  empty  trap  2 states. 

2 l2  l2 

With  the  parameters  provided  in  Table  3.1  and  the  electron  density  in 
the  conduction  band,  it  can  be  shown  with  the  help  of  a simple  computer 
program  that  the  inequality  holds  for  the  range  of  temperatures 
considered  in  our  experiments.  We  would  like  to  make  a final  remark  on 
the  interpretation  of  Fig.  3.10.  It  is  possible  to  sea  a third  noise 
peak  with  a tentative  energy  level  of  0.2  eV  and  trap  density  of 

8 x 10  m . The  occurrance  of  a third  trap  might  explain  the  more 
gradual  decrease  of  the  electron  density  as  1/T  increases  without 
reaching  a plateau  in  the  temperature  range  below  200  K (Fig.  3.7). 
However,  if  we  look  at  Figs.  3.3  and  3.9,  we  see  that  these  data 

probably  have  no  physical  meaning  but  are  artifacts  of  decomposing  the 
spectra. 

The  ultimate  test  of  our  analysis  is  simulating  the  noise  data  by 
writing  a computer  program  that  describes  the  noise  exactly,  as  men- 
tioned in  ref.  [28],  using  the  trap  parameters  from  Table  3.1.  Figures 
3.8  and  3.9  show  for  comparison  both  the  measured  and  the  simulated 
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data.  v*/e  find  that  the  simulated  data  fit  the  measured  data  very  well, 
leaving  no  doubt  in  the  interpretation  concerning  these  two  traps. 

3.5.  Double  Donor  Model 

For  reasons  given  in  the  introduction  we  consider  in  this  section 
the  case  of  an  impurity  that  can  donate  two  electrons,  a so-called 
double  donor.  As  explained  in  the  preceding  sections,  noise  measurements 
do  not  discriminate  on  the  electrical  nature  of  trap  centers.  Con- 
sequently, a double  donor  and  a double  acceptor  give  the  same  noise 
spectra.  In  this  section  we  opted  to  describe  the  double-donor  model 
since  this  situation  is  more  likely  to  occur  in  n-type  Si.  Analogous  to 
what  is  described  in  sec.  3.2,  one  can  derive  expressions  for  the  char- 
acteristic times  and  zero-frequency  plateau  levels  that  completely 
describe  the  g-r  noise  spectra  as  a function  of  temperature  for  a 
double-donor  impurity,  using  modified  rate  equations  and  compounded 
Fermi-Dirac  statistics  (see,  for  instance,  van  Vliet  [32]).  To  dis- 
tinguish between  the  intended  dopants  (donors  here)  and  the  unintended 
dopants  we  call  impurities  (also  donors  here),  the  first  have  a 
subscript  D and  the  latter  have  superscripts  o,  +,  and  ++.  Ue  denote 
the  number  of  neutral,  singly  charged,  and  doubly  charged  impurities  by 
d®,  N+,  and  N++,  respectively.  The  total  number  of  impurities  is  N. 
Then  the  rate  equations  read: 


dn 

dt  “ S1  ~ rl  + s2  r2 

(3.68a) 

dN+ 

dt  gl  rl  “ g2  + r2 

(3.68b) 

dN**  _ 

dt  g2  r2 

(3.63c) 
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where  the  gj  and  rt  are  described  in  Table  3.2.  Neglecting  the  contri- 
bution of  the  holes  to  the  space  charge,  the  condition  of  space-charge 
neutrality  is  given  by 


where  nD  is  the  number  of  electron-occupied  intended  donors.  For 


Equation  (3.69)  is  used  to  eliminate  one  of  the  three  dependent 


spectral  density  of  current  or  voltage  fluctuations  is  directly  related 
to  fluctuations  in  the  number  of  mobile  electrons,  n is  kept  as  an 
independent  variable.  The  choice  of  the  other  independent  variable  is 
arbitrary. 

In  the  same  notation  as  in  sec.  3.3  we  find  the  elements  of  the 
relaxation  matrix  M 


4*  ++ 

n = N + 2N 


(3.69) 


shallow  donors  n^  « N . 


variables  (n,  N+,  and  N++)  from  the  rate  equations.  Since  the  measured 


(3.70a) 


(3.70c) 


(3.70b) 


(3.70d) 


The  subscript  t stands  for  taking  the  derivatives  in  thermal  equili- 
brium. The  measured  corner  frequencies,  from  the  spectra,  times  2tt  are 
the  eigenvalues  X^  of  the  matrix  M.  The  noise  spectra  are  given  by 
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Table  3.2 

Generation  and  Recombination  Rates  in  Terms  of 
(Ionized)  Impurity  Numbers  and  Free  Electron  Numbers 


Rate 

Bimolecular 
"Chemical"  Reaction 
Involved 

Description 
of  the  Rates 

§1 

1 

1 

o 

25 

> N+  + e 

a1(N-N+-N++)nsr,1 

rl 

N+  + e — 

1 

1 

1 

V 

% 

a^N+n 

g2 

N+  - 

> N++  + e 

^9 

r2 

N++  + e — 

> N+ 

a?N++n 

92 


S (w)  = — — 
n 2 


2 "r  2 2 
a + X^  u U, 


(3.71) 


where 


'1  ‘ [1<M12-M22>2  - ¥ <W  + U:'12+f,22)2  - ¥ <¥>] 

(3.72a) 


X2  - X1 


2 2 2 ^X2  ^M12_M22^  ^ <?1?1>  + (X2  _ 


X2  - X1 


12^22^  ^ <?2F2^ 
(3.72b) 


<F1?1>  = 4gl 


(3.72c) 


<?2F2>  = 4S2 


( 3 . 72d ) 


The  number  of  singly  and  doubly  Ionized  Lmpurities  can  be  calculated 
using  the  results  Champness  [33]  obtained. 


N+  = N 


„„  rE 1 ~ EF> 
28 1 e*p(-kgr-) 


Ei  ' ef. 


1 + ZBj  exp ( — T ) + BxB2  exp(- 


E1  + E2  - 2V 


V 


(3.73a) 


N = N 


E + E - 2E 
!1B2  «*»(- i^T ) 


E,  " Ev  E + E.  - 2E  • 

1 + 2B1  exp(_k^-'  + B1B2  axp( i£f £) 


(3.73b) 


Ei  and  E2  are  associated  with  the  singly  and  doubly  Ionized  Impurity 
levels,  respectively. 
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From  a theoretical  point  of  view  one  would  expect  a difference  in 
applying  the  results  derived  here  as  compared  to  the  two  single-level 
impurities.  In  the  double-donor  case  only  one  trap  density  has  to  be 
found  as  opposed  to  the  two  densities  that  need  to  be  found  in  the  two- 
single-level  case.  However,  when  we  simulate  the  noise  spectra,  speci- 
fically the  zero-frequency  plateau  values  and  the  characteristic  times, 
using  the  parameters  of  Table  3.1  and  making  N_  = NT  , we  find  the  same 
excellent  agreement  with  the  measurements  for  the  double-donor  model  as 
in  the  case  of  two  single-level  impurities.  The  reason  that  we  find  no 
significant  difference  between  the  simulation  results  of  the  two- 
independent  -trap  model  (model  I)  and  the  double-donor  model  (model  II) 
is  twofold.  Firstly,  in  the  case  of  the  double  donor  the  possibility 
for  a second  electron  to  leave  the  donor  only  exists  when  the  Impurity 
is  already  singly  ionized.  In  other  words,  transitions  between  the 
conduction  band  and  the  lower,  second  energy  level  only  occur  if  the 
higher,  first  one  is  emptied.  But  this  condition  is  also  fulfilled  In 
the  case  of  the  two  independently  interacting  traps  when  they  are 
separated  by  more  than  kgT.  Only  when  the  energy  levels  are  closely 
spaced  will  changes  in  the  occupancy  of  one  of  the  levels  be  strongly 
correlated  with  changes  in  the  occupancy  of  the  other  level,  and  now  the 
availibility  of  other  states  for  the  electron  to  go  to  should  make  it 
possible  to  distinguish  between  the  two  models. 

Secondly,  although  one  less  parameter  has  to  be  found  in  the 
analysis  applied  to  a double-donor  model  (only  one  trap  density  instead 
of  two)  this  subtle  difference  might  be  obscured  by  the  temperature 
dependence  of  the  capture  cross  section.  No  conclusive  evidence  is 
found  in  favor  of  either  one  of  the  models  to  describe  the  measured 
data,  both  agree  very  well. 


CHAPTER  IV 

MONTE  CARLO  CALCULATIONS 

4.1.  Introduction 

In  the  history  of  semiconductor  devices,  a major  goal  has  always 
been  to  achieve  short  response  times  to  enable  the  increase  of  the 
number  of  logic  operations  in  digital  circuits  and  the  increase  of  the 
bandwidth  in  analog  amplifiers.  Improved  processing  techniques,  along 
with  higher  degrees  of  integration,  have  greatly  contributed  to  achiev- 
ing the  higher  device  speeds.  The  latest,  almost  revolutionary  develop- 
ment has  been  the  reduction  of  the  dimensions  of  the  devices  to  the 
point  where  the  description  of  the  charge  transport  can  be  properly 
given  only  in  terms  of  quantum  mechanics,  rather  than  classical  physics. 
The  layers  of  semiconductor  material  have  shrunk  to  thicknesses  as  small 
as  10  atomic  layers.  The  dramatic  impact  of  the  size  quantization  has 
triggered  both  experimentalists  and  theorists  to  explore  the  possible 
applications.  The  device  that  has  received  the  greatest  attention  is 
the  modulation-doped  field-effect  transistor  (MODFET).  The  channel  in 
which  the  conduction  takes  place  has  a width  of  the  order  of  100  A in 
the  device.  The  typical  gate-to-channel  distance  is  of  the  order  of 
300  A in  these  transistors. 

Presently,  several  techniques,  such  as  metal  organic  chemical 
vapor  deposition  (MOCVD),  liquid  phase  epitaxy  (LPE),  molecular  beam 
epitaxy  (MBE),  etc.,  are  available  for  the  fabrication  of  these 
devices.  For  theorists  a new  territory  needed  to  be  covered.  The  size 
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quantization  gives  a new  dimension  to  the  description  of  the  charge 
transport.  i.he  Boltzmann  equation  needs  to  be  solved  taking  quantum 
mechanical  effects  into  account.  The  confinement  of  the  charge  carriers 
is  reflected  in  a quantization  of  the  wave  vector  of  the  carrier.  Two 
direct  implications  of  the  quantization  of  the  wave  vector  of  the  charge 
carr^er  are  the  different  densities  of  states  and  the  different  rates  at 
which  the  carriers  interact  with  phonons  and  other  carriers. 

Before  large  computer  systems  were  available  and  numerical  solu- 
tions to  various  complicated  mathematical  problems  could  be  found,  the 
only  approach  for  solving  the  Boltzmann  transport  equation  (3TE)  was  the 
analytical  way.  The  advantage  of  analytical  solutions  is  that  they  can 
be  described  by  explicit  mathematical  formulae.  However,  finding 
analytical  solutions  for  the  3TE  is  only  possible  for  some  particular 
cases  of  high  symmetry  and  simplicity. 

A few  solution  techniques  will  be  discussed.  One  of  the  tech- 
niques used  is  to  expand  the  distribution  function,  which  is  the  solu- 
tion of  the  3TE,  in  Legendre  polynomials  [9]: 


f(k)  = f0(E)P0(cos9)  + f1(E)P  (cos0)  + (4.1) 

It  can  then  be  shown  that 

fx(E)  = eFx( 2E/m*) 1/2  (4.2) 


where  E is  the  kinetic  energy  of  the  particle,  k its  wave  vector,  e the 
charge  of  an  electron,  t the  relaxation  time  of  the  scattering  pro- 
cess, F the  applied  electric  field  strength,  and  m"  the  particle's 
effective  mass.  In  sections  4.1,  4.2,  and  4.2.1,  the  small  k denotes 
the  three-dimensional  wave  vector.  In  the  other  sections  of  this  chap- 
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ter  small  k is  reserved  for  the  two-dimensional  wave  vector  parallel  to 
the  hetero  junction  interface.  Then  the  assumption  is  made  that  the 
applied  electric  field  is  small  and  the  zero-order  velocity  distribution 
function  (fg)  is  Gaussian.  Macroscopic  observables  like  the  drift 
velocity  and  diffusion  coefficient,  parameters  that  determine  the  charge 
transport,  can  be  found  from  the  first  and  second  moments  of  the 
distribution  function.  However,  for  some  scattering  mechanisms  the 
relaxation  time  cannot  be  defined.  An  example  is  the  polar  optical 
phonon-electron  interaction.  Accordingly,  this  technique  cannot  be 
employed  to  describe  charge-transport  properties  of  (III-V)  and  (II-VI) 
semiconductors . 

Another  approach  is  that  of  forcing  the  solution  of  the  BTE  to  fit 
a displaced  Maxwellian  [34] 


f(k)  = const  exp[- 


k-k„ 


knT 
B e 


2m 


-)] 


(4.3) 


where  is  Planck's  constant  divided  by  2m,  kQ  is  Boltzmann's  constant, 
and  kjj  is  the  center  of  the  displaced  distribution  function.  The  tem- 
perature Tg  is  now  the  electron  temperature  which  might  be  elevated  with 
respect  to  the  lattice  temperature.  The  two  degrees  of  freedom  T and 
are  found  from  the  first  and  second  moments.  Again,  only  approximate 
solutions  can  be  found.  The  extension  to  time—  and  space— dependent 
phenomena  is  not  easily  available.  Also,  electron-electron  interaction 
is  difficult  to  implement  since  this  type  of  interaction  makes  the  BTE 
nonlinear. 

The  next  technique  is  an  iterative  process  to  obtain  a solution  of 
the  BTE  [35].  The  method  needs  fast  computers  to  fully  exploit  the 
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technique,  but  the  result  will  be  an  exact  solution.  Also,  time- 
dependent  phenomena  can  be  easily  included.  The  BIT  reads 


_ r3f(k,t) ' , r3f(k,tK 

3t  1 at  J + l— — J 

coll  field 


(4.4) 


The  evolution  of  the  distribution  of  particles  in  wave  vector  space  and 
time  is  determined  by  the  rate  at  which  the  particles  suffer  collisions 
and  the  rate  at  which  energy  from  the  applied  electric  field  is  supplied 
to  the  carriers.  The  collision  term  can  be  split  up  into  two  terms: 

/ f(k'  ,t)S(k,ic'  )dk'  - A(k)f(k,t)  . (4.5) 


The  first  term  of  eq.  (4.5)  is  the  rate  at  which  the  distribution  func- 
tion changes  due  to  scattering  of  particles  with  an  initial  wave  vector 
k'  to  a state  of  wave  vector  k.  5(k,k')  is  the  transition  probability 
from  a state  k'  to  k.  The  second  term  of  eq.  (4.5)  represents  the  rate 
at  which  particles  with  initial  wave  vector  k are  scattered  to  any  other 
state, 

X(k)  = / S(k,k')dk'  (a 


In  the  field  term  the  partial  derivative  with  respect  to  time  can  be  re 
placed  by  — . Then,  a formal  integration  yields  [35] 


f (k) 


=/  exp { / A(k  - f-  x)dx}  / S(k  - f-  t",k')f(k')dk'df 


Now  write 


(4.7) 


and 


g(k)  = / S(k,k' )f (k1 )dk' 
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f(k)  = / g(k  - j-  t)exp{  - / X(k  - x )dx}dt  ; 

0 n 0 A 

then  gn+1  = Af^  and  fn+1  = Bgn+1  in  operator  notation.  In  this  scheme 
the  development  of  f is  monitored  in  time,  fn+1,  being  the  distribution 
function  at  a time  t^+^ . 

An  even  more  powerful  technique  is  the  Monte  Carlo  method.  Numer- 
ically the  technique  is  very  simple  and  no  integrations  are  required. 
Both  the  space  and  time  development  of  the  distribution  function  can  be 
simulated.  Again,  this  method  requires  a large  and  fast  computer.  In 
the  Monte  Carlo  method  used  in  the  transport  theory,  the  calculation  of 
the  distribution  function  is  replaced  by  a computer  simulation  of  the 
motion  of  individual  particles.  These  particles  then  are  treated  as 
sample  members  of  an  ensemble.  Expectation  values  of  physical  quanti- 
ties, such  as  the  drift  velocity  and  diffusion  coefficient,  are  given  by 
time  averages.  The  distribution  function  is  obtained  by  counting  the 
number  of  appearances  of  a variable  value  in  each  "histogram  interval." 
The  motion  of  a mobile  charge  carrier  in  a solid  is  a sequence  of 
ballistic  flights  governed  by  the  applied  external  fields  and  terminated 
by  scattering  processes  determined  by  phonon  interactions.  Both  the 
duration  of  the  intercollision  flight  and  the  choice  of  the  scattering 
process  terminating  the  ballistic  flight,  as  well  as  the  determination 
of  the  wave  vector  after  scattering,  are  treated  as  random  variables 
with  given  distributions.  These  random  variables  are  found  from  deter- 
minate functions  of  one  or  more  random  numbers. 

The  success  of  the  Monte  Carlo  method  is  very  much  related  to  the 
fact  that  particle  variables  can  be  explicitly  computed  without  having 
to  calculate  the  distribution  function  first.  In  this  method  com- 
plicated paths  and  physical  boundaries  can  be  included,  as  well  as 
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complicated  generation,  recombination,  and  scattering  mechanisms.  The 
generated  particle  states  tend  to  favor  the  ranges  of  particle  variables 
that  are  appropriate  to  the  situation  being  simulated.  A serious  dis- 
advantage of  the  Monte  Carlo  method  is  that  the  computed  results  are 
estimators.  They  are  random  variables  themselves.  Further  calculations 
do  not  steadily  decrease  the  difference  between  estimator  values  and 

true  values.  Consequently,  convergence  is  difficult  to  monitor  and 
assess. 

In  the  next  section  the  details  of  the  Monte  Carlo  method  are  ex- 
plained. In  section  4.3  we  focus  on  the  application  of  this  technique 
to  the  description  of  the  charge  transport  in  an  infinitly  high  quantum 
well.  In  the  final  section  we  present  results  for  the  drift  velocity 
and  diffusion  coefficient  as  a function  of  the  electric  field  and  well 
width.  The  outcome  of  our  simulations  is  then  discussed. 

^ • 2 . The  Monte  Carlo  Procedure 

To  describe  the  structure  of  a typical  Monte  Carlo  program,  the 
focus  will  be  on  the  motion  of  an  electron  in  a solid.  To  be  more  pre- 
cise, the  simulation  of  a stationary  transport  process  under  the  con- 
dition of  a homogeneous  electric  field  will  be  considered.  The  simula- 
tion starts  with  one  electron  in  given  initial  conditions  with  wave 
vector  kQ.  Then  the  duration  of  a free  flight  under  the  influence  of 
the  applied  fields  is  chosen  with  a probability  distribution  determined 
by  the  scattering  probabilities.  Each  electron  state  is  specified  by 
the  wave  vector  k of  the  electron  and  hence  a momentum 


P = Wc 


(4.3) 


can  be  defined.  The  energy  of  the  electron  will  be  a function  of  the 
wave  vector  E(k)  or  E(p).  During  the  flight  under  the  influence  of  the 
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external  forces  the  equation  of  motion  describes  the  change  of  momentum 


of  the  electron 


Force 


(4.9) 


where  Force  is  the  externally  applied  force,  which  may  be  taken  as  given 
by 


Force  = -q(F  + - v * B)  - 7E(p,r) 


(4.10) 


where  F and  B are  the  applied  electric  and  magnetic  field,  respectively. 
The  last  term  of  eq.  (4.10)  is  due  to  a strain  gradient  in  the  solid 
(such  as  due  to  phonons)  or  an  inhomogeneity  of  composition  and  hence  of 
the  band  structure  [36].  We  drop  this  term  because  we  only  consider 
homogeneous  semiconductors.  Furthermore,  we  choose  B = 0 and  take  the 
electric  field  constant  in  time  and  space.  Under  these  conditions,  the 
wave  vector  of  the  electron  changes  during  the  ballistic  flight  accord- 
ing to: 


where  tf  is  Planck's  constant  divided  by  2m  and  q,  the  numerical  value  of 
one  electron  charge.  The  change  of  the  wave  vector  results  in  an  evolu- 
tion of  the  velocity  (v  = V^E  where  the  subscript  p denotes  taking  the 
derivative  in  momentum  space),  the  energy,  E,  and  other  quantities  of 
interest.  These  developments  are  recorded  and  a scattering  mechanism  is 
chosen  ending  the  flight.  The  choice  of  a particular  scattering  mechan- 
ism is  governed  by  the  relative  probability  of  occurrence  of  this  pro- 
cess with  respect  to  the  other  scattering  processes  considered.  Next  a 
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new  k state  is  selected  depending  on  the  particular  scattering  mechan- 
ism, and  this  state  is  the  initial  wave  vector  for  the  following  drift, 
and  the  whole  process  is  repeated.  The  total  simulation  time  is  divided 
into  successive  subhistories  of  equal  time  duration,  and  the  quantities 
of  interest  are  calculated  for  each  subhistory.  This  method  allows 
calculating  the  value  for  the  quantity  by  averaging  over  the  sub- 
histories, but  more  importantly,  estimating  the  uncertainty  in  this 
value  by  computing  the  standard  deviation. 

4.2.1.  Determination  of  the  flight  time 

The  subject  of  this  section  is  to  describe  how  the  length  of  an 
intercollision  flight  time  is  found  from  generating  a random  number. 
Obviously,  there  is  not  a linear  relation  between  the  two  quantities 
since  the  flight  time  can  take  on  any  value  from  zero  to  infinity, 
whereas  a random  number  is  uniformly  defined  on  a finite  interval,  most 
commonly  on  [0,1],  In  general,  if  p ( 9 ) and  p(r)  are  the  probability 
densities  associated  with  the  physical  distribution  and  the  random 
number  distribution,  respectively,  then 

9 r 

/ p(9')d9'  = / p(r’  )dr ' . (4.12) 

0 0 

Specifically,  in  the  case  that  the  physical  quantity  of  interest  is  the 
time  t and  p(r)  is  uniform,  one  finds 

t 

r = / p(t'  )dt'  , (4.13) 

0 


which  is  an  implicit  expression  for  t. 
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Suppose  now  that  every  flight  time  t is  the  sum  of  incrementally 
small  time  intervals 

n 

1 = E 5ti  • (4.14) 

The  probability  that  the  electron  is  scattered  within  the  time  interval 
6t1  is  equal  to  the  product  of  the  total  scattering  rate  and  the  length 
of  that  time  interval: 

A(k)5t1 

where  A(k)  is  the  rate  at  which  electrons  with  wave  vector  k scatter. 

N N 

*(k)  = l A (k)  = l f S (k,k')dk*.  (4.15) 

j=l  3 j=l  J 

In  eq.  (4.15)  S^(k,k')  is  the  transition  rate  from  a state  with  wave 
■> 

vector  k to  a state  with  wave  vector  k'  caused  by  scattering  mechanism 
j.  We  consider  N different  scattering  processes.  Clearly,  A . (k)  is  the 
scattering  rate  from  an  initial  state  k to  any  other  possible  final 
state  due  to  process  j.  The  probability,  P(t),  that  an  electron  travels 
a time  t without  suffering  any  collisions  equals  the  product  of  the 
probabilities  that  an  electron  does  not  scatter  in  time  interval  5t 
which  is  1 - A(k)5ti,  hence 

n 

P(t)  = n ( l-A (k)5t  ) . (4. 16) 

1=1  1 

Taking  the  natural  logarithm  and  realizing  that  log(l-x)  = -x  for 
x « 1,  one  finds 
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t 

P(t)  = exp{-  / A(k)dt ' } . (4.17) 

0 


Note  that  k is  a function  of  time,  see  eq.  (4.11).  The  probability  that 
an  electron  travels  a time  t without  any  collisions,  and  then  at  time  t 
scatters,  is  given  by  the  probability  density 

t 

p(t)  = A(k)  exp{-  / A(k)dt'}  . (4.13) 

0 

Then,  by  using  eq.  (4.13)  we  find 


t 

r = / p(t' )dt ' 

0 

t t' 

= / A(k)  exp{-  / \(k)dt"}dt ' 

0 0 

t t' 

= / ~ TTT  [exp{-  / A(k)dt" } ]dt ' 

0 0 

t ’ + t ' =t 

= - exp{-  / A(k)dy"} 

0 t ' =0 

t 

= 1 - exp{-  / A(k)dt"}  . (4.19) 

0 


It  is  not  possible  to  find  an  explicit  form  that  expresses  t as  a 
function  of  the  random  number  r.  Recall  that  A(k)  is  the  sum  of  all 
scattering  rates  A^(k).  Suppose  that  we  add  another  scattering  rate, 
» a virtual  process  that  does  not  change  the  electron  wave  vector, 


(j  - 0,  1,  2,  ...  N)  is  a constant,  T 


N 

I A.(k)  = r . 

j=0  3 


(4.20) 
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We  can  calculate  Aj  (k)  for  j = 1,  2,  . . . ,N  since  they  describe  real 
physical  processes  for  which  the  rates  are  known.  Choose  r not  smaller 
than  the  maximum  of  the  sum  of  the  real  scattering  rates: 


N 

r > max  l A (k)  , 
{k}  j=l  3 


(4.21) 


then  X^(k) 


can  be  found  from 


x0(k)  = r 


N 

I X (k)  . 
j = l 3 


(4.22) 


The  probability  of  the  occurrence  of  a virtual  scattering  process  can  be 
determined  since  all  of  the  rates  are  known. 

Eq.  (4.19)  will  now  read 

t 

r = 1 - exp{ - / A(k)dt "} 

0 

t N 

= 1 - exp{ - / [ l X (k)  + X (k)]dt" 

0 j=l  3 U 

t 

= 1 - exp{-  / rdt”}  , (4.23) 

0 

which  can  be  readily  integrated  and  t is  found  to  be 

f = f • (4.24) 

After  the  electron  drifts  for  a time  t,  in  accordance  with  eq.  (4.24), 
the  electron  is  scattered  either  by  a real  process  with  rate  X^(k) 

(j  = 1,  2,  ...  N)  or  by  a virtual  process  with  rate  XQ(k).  To  deter- 
mine which  scattering  mechanism  is  chosen,  all  of  the  rates  are  normal- 
ized by  dividing  by  T.  The  sum  of  the  normalized  rates  (sum  over  j = 0, 
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•••  N)  equals  1,  and  each  normalized  rate  X.(ic)/r  e [0,1].  The 
interval  [0,1]  is  divided  into  subintervals  of  lengths  A (k)/T  (see 
Fig.  4.1).  The  scattering  process  with  the  highest  rate  occupies  the 
largest  subinterval,  and  when  a random  number  is  generated  with  a 
uniform  distribution  on  interval  [0,1],  ic  will  have  the  most  likelihood 
of  falling  in  the  interval  associated  with  that  rate.  The  scattering 
rate  A (k)  is 

X0(k)  = / So^fic'Jdk'  (4.25) 

where 

Sg(k,k')  = X^(k ' )5(k-k' ) . (4.26) 

The  last  factor  on  the  right-hand  side  is  the  Dirac  delta  function  which 
shows  that  the  wave  vector  does  not  change  in  the  process,  hence  the 
name  virtual  process. 

The  introduction  of  the  virtual  process,  or  self-scattering 
process,  has  the  advantage  that  eq.  (4.19)  can  be  solved  very  easily, 
dowever,  since  the  virtual  processes  do  not  contribute  to  the  estimation 
of  the  transport  quantities,  the  price  that  has  to  be  paid  is  the  com- 
puting time  spent  on  the  evaluation  of  these  processes.  Therefore,  it 
is  advisable  to  keep  T as  small  as  possible  but  to  avoid  letting 
Xg(k)  become  negative. 

After  a real  scattering  mechanism  has  been  chosen,  the  final  wave 
vector  state,  that  is  the  state  after  scattering,  has  to  be  determined. 
The  discussion  of  how  the  scattering  rates  and  the  final 

-states  are  generated  continues  in  the  next  two  sections,  4.2.2  and 
4.2.3. 
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Fig.  4.1. 


Selection  of  the  scattering  channel 
terminating  the  electron  flight. 
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4.2.2.  Scattering  rates  for  a two— dimensional  electron  gas 
Heterolayer  semiconductor  materials,  with  alternating  layers  of 
different  semiconductors,  tens  to  hundreds  of  angstroms  thick,  have 
become  the  subject  of  many  experimental  and  theoretical  research  pro- 
jects. The  electron  states  in  these  structures  differ  from  the  conven- 
tional Bloch  states  and  result  in  the  following,  interesting  electronic 
phenomena:  (i)  formation  of  superlattice  bands;  (ii)  the  introduction 

of  energy  subbands  due  to  quantization  of  the  wave  vector  perpendicular 
to  the  interfaces;  (iii)  formation  of  a two-dimensional  electron  gas 
(2DEG)  at  the  heterojunction  interfaces.  This  chapter  is  mainly 

concerned  with  the  electron  motion  in  the  two-dimensional  gas,  as  it 
shows  up  in  AlGaAs/GaAs  devices. 

In  Fig.  4.2  the  doping  profile  and  the  layer  configuration  are 
sketched  for  several  periods  of  a superlattice.  The  alternating  doping 
profile  is  called  modulation  doping.  In  these  structures  the  material 
with  the  larger  bandgap  (AlGaAs)  is  heavily  doped  (~  1024  m'3),  and  the 
mobile  electrons  donated  by  the  ionized  dopants  spill  over  to  the  lower 
bandgap  material  since  this  Is  energetically  more  advantageous.  The 
separation  of  the  electrons  from  their  donors  gives  rise  to  an  electric 
field  across  the  interface.  Due  to  the  electric  field,  the  conduction 
band  bends,  the  curvature  is  concave  where  the  space  charge  is  positive 
(in  the  doped  AlGaAs),  and  convex  where  the  space  charge  is  negative, 
i.e.,  in  the  GaAs.  The  resulting  band  profile  for  one  interface  is 

shown  in  Fig.  4.3.  The  separation  of  the  electrons  from  their  donors 
reduces  the  impurity  scattering  and  thus  Improves  the  electron  transport 
in  the  GaAs  channel  parallel  to  the  Interface.  A further  reduction  of 
this  Coulomb-like  interaction  is  obtained  by  growing  an  undoped  AlGaAs 
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doping 


Fig.  4.2.  Typical  doping  profile  and  conduction  band 
profile  for  a heterolayer  structure. 
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Fig.  4.3. 


Qualitative  band  profile  of  one 
heterointerface. 


no 


layer  between  the  doped  AlGaAs  and  the  undoped  GaAs,  in  this  way 
increasing  the  separation  between  the  donors  and  the  electron  gas.  This 
spacer  layer,  on  the  other  hand,  must  not  be  too  thick  because  it  ob- 
structs the  transfer  of  electrons  from  the  doped  AlGaAs  to  the  GaAs. 

For  a high  doping  level  the  band  bending  is  large,  resulting  in  a 
strong  confinement  of  the  electrons  in  the  vicinity  of  the  interface. 
If  in  addition  the  width  of  the  GaAs  layer  becomes  of  the  order  of 
200  A,  then  quantization  effects  of  the  wave  vector  need  to  be  taken 
into  account,  and  solutions  of  the  Schrodinger  equation  need  to  be  found 
in  addition  to  solving  Poisson's  equation  to  account  for  the  band  pro- 
file. The  simultaneous  solution  of  these  equations  is  possible  but  very 
complicated.  To  gain  physical  insight  into  the  electron  transport  in 
these  systems,  we  chose  a model  in  which  the  band  profile  of  the  GaAs  is 
approximated  by  an  infinitely  high  quantum  well.  For  this  system  the 
energy  levels  and  the  wave  functions  are  well  known,  and  this  simplifies 
the  analysis  greatly.  It  is  for  this  system  that  the  interactions 
between  electrons  and  phonons  are  considered.  Although  the  quantization 
of  the  electron  wave  vector  is  taken  into  account,  the  phonon  gas  is 
still  considered  to  be  three-dimensional  since  the  lattice  matching 
between  the  AlGaAs  and  GaAs  Is  very  good. 

Analytic  expressions  for  the  scattering  rates  for  acoustic  phonon 
interaction  via  the  deformation  potential,  optical  (or  intervalley) 
phonons  via  the  deformation  potential,  and  optical  phonons  via  the  polar 
Interaction,  are  presented.  For  a more  detailed  discussion,  see  Riddoch 
[37].  An  analogous  discussion  can  be  found  in  ref.  [36].  We  consider  a 
parabolic  band;  the  effective  mass  is  independent  of  the  wave  vector. 
The  direction  of  confinement  is  along  the  x-axis,  while  the  electron 


1 1 1 


transport  takes  place  along  the  z-axls.  The  electron  wave  function 
«(x,r)  vanishes  at  x = 0 and  x = L for  an  Infinitely  high  well.  Hence 
the  normalized  wave  function  is 


1/2 


>Kx,r)  - ( ) u+(x,r)  exp(ik*r)  sin(k  *x) 


(4.27) 


where 


and 


and 


K (kx’VV  " (kx’k) 


R = (x,y,z)  = (x,r) 


, niT 

k = — . 
x L 


V is  the  volume;  u^  is  the  Bloch  function.  The  energy  of  an  electron  is 
K 

given  by 


E = 


Jl_  v2  _ tf2  ^2 
K ; — * <Vk> 

2m  2m 


,,2  2 2 

- A [=4-  + *2] 

2m  L 


2 2 

0_  . 2 

r + n E 

2m"  1 


(4.28) 


The  density  of  states  is  equal  to  the  number  of  energy  or  wave-vector 


states  per  unit  volume. 
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N(E+)dE 

K ! 


g(E+)dE 

K 


| = g (K)dK 
V 


K 


V 


(4.29) 


* 


k 


(4.30) 


The  two-dimensional  density  of  staces  per  unit  energy  is  sketched  in 
Fig.  4.4.  The  scattering  rate  for  a particular  electron-phonon  scatter- 


where Hep  Is  the  Hamiltonian  describing  the  electron  phonon  interaction, 
and  the  matrix  element  describes  the  probability  of  an  electron  with 


ence  of  interaction  Hep.  The  6-function  provides  for  energy  conserva- 
tion during  the  process,  and  hu^  is  the  energy  of  the  phonon  involved. 
The  integration  is  performed  over  the  available  final  states.  The 
matrix  element  for  each  scattering  process  of  interest  can  be  expressed 
as  follows  [37] : 


ing  process  represented  by  Hep  is  [37] 


A(K)  = ~ f S (K , K ' )dK' 
h J ep 

2tt  r I | . ■+•  , 9 

- r I l<K’ !«  IK>I  5<E>  - 

* Lr  1 v 


K*  K 


(4.31) 


initial  wave  vector  K being  scattered  to  final  state  K'  under  the  influ- 


!<^*lHep|K>!2 


= cons 


t[n(coQ)  + 5^G2(qx)  . 


1 _ 1 


(4.32) 


q 
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Fig.  4.4.  Two-dimensional  density  of  states  compared 
with  three-dimensional  density  of  states. 
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The  constant  Is  comprised  of  a factor  , the  overlap  integral,  and 

the  coupling  strength.  M is  the  oscillator  mass,  N the  number  of  lat- 
tice modes,  and  the  angular  frequency  of  the  lattice  wave.  The 
Kronecker  symbol  provides  for  conservation  of  momentum  in  the  (y,z) 
plane,  k'  = k ± q . The  upper  sign  is  for  phonon  absorption  and  the 
lower  for  phonon  emission.  The  last  factor  in  eq.  (4.32)  is 


G<V  -f  / 


L iq  x 

e x sin(k'x)  sin(k  x)dx 

X X 


(4.33) 


•■/hen  calculating  the  matrix  element,  one  needs  to  compute 

/ K K K • (4.34) 

K'  Q K 

In  the  case  of  a three-dimensional  electron  gas,  this  results  in 
/ exp(+iIC'*R)  exp(-iQ«R)  exp(-iK*R)dR 
= <5(K'  - $ - k) 

= 6(k'  - q - k) 6(k^  - qx  - k*)  . (4.35) 

It  is  clear  that  this  expression  describes  the  momentum  conservation. 
In  the  case  of  a two-dimensional  electron  gas,  expression  (4.34)  will 
yield 

fK  ^ ^ dR  = / / exp(ik'*r)  sin(!y!-)  exp(-iq-r)  exp(-iq  x) 

K'  Q K 0 L x 


> ,nTTs  •> 

x exp(-ik*r)  sin(— Jdr  dx 


sinfx  -^1 

= 6(k'-q-k)  • l { — exp( i k y)} 

+/-  ‘x  2 


(4.36) 
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where  <x  - qx  ± (kx  ± k^)  and  the  summation  Is  performed  over  the  four 
possible  combinations  of  + and  - signs  in  k . sin(<  x)/(<  reaches 

a maximum  for  = 0,  which  means  that  conservation  of  the  momentum  in 
the  x direction  (quantization  direction)  is  still  preferred.  However, 
there  is  a finite  contribution  to  the  matrix  element  for  final  wave 
vectors  K'  that  do  not  conserve  the  momentum  in  the  x-diraction,  i.e. , 
for  * 0.  In  the  limit  that  the  width  of  the  well  goes  to  infinity 
(L  ^ °°),  sin(<x  j) / ( tends  to  the  5-function  which  is  consistent 
with  the  three-dimensional  description,  and  momentum  is  conserved  in  the 
x-direction.  However,  in  the  limit  that  L goes  to  zero,  the  width  of 
the  interval  of  kx  that  contributes  significantly  to  the  matrix  element 
is  large  and  momentum  conservation  does  not  hold  strictly,  but  is  rather 
fuzzy.  In  the  cases  that  k'  = k or  k'  = -k  , two  of  the  terms  in  aq. 
(4.36)  maximize  at  the  same  time,  that  is  for  qx  = 0.  These  two  cases 
are  associated  with  intrasubband  (n'  = n)  scattering;  in  the  former  case 
the  wave  vector  is  not  reflected,  while  in  the  latter  case  it  is.  In 
calculating  | G(q^) | , we  will  distinguish  between  these  two  cases.  To 
find  tractable  analytical  expressions  for  the  scattering  rates,  a few 
approximations  have  to  be  made.  The  interference  between  the  four 

processes  will  be  neglected,  except  when  k^  = ± 1^.  Then,  in  the  case 

kx  = * V 


lG(qx) 


"x  2 


<v2V  i 


Sln( (q+2k)  |) 

+ [ * *,  2 ]2} 


<V2kx>  i 


and  for  k'  * ± k 
x x 


(4.37) 
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G(qx} 


j 1 [si°t(V(k^x)>  ?)]2 
+/-  4 (,xi(k’ikx»  i 


(4.38) 


This  truncation  is  referred  to  as  the  momentum  conservation  (MCA)  [37]. 
The  second  approximation  is  to  let  |G2(qx)|  act  as  a delta  function  in 
such  a way  that  for  any  function  f(qp 


/ f(<i;>|G2(q  )|dq'  * / f(q’)|G2(q  )|dq' 

q * x _oo  x xx 

max 


2lT  F ( 

r f(qx 


>7  I 


+/- 


(5 


q ±(k'±k  )• 
x x x 


(4.39) 


More  specifically,  if  we  choose 


f (lx)  - -2A(<l)2  • (4.40) 

qx  + q 

the  deformation  potential  case  is  described  for  A(q)  = q2  when  q is  very 
large.  The  polar  coupling  is  described  by  choosing  A(q)  = 1 and  keeping 
q finite.  The  calculation  of 

/ f<qx)|G(qx)|2dqx 

where  G(qx)  is  defined  in  eq.  (4.36)  is  tedious  but  straightforward.  It 
can  be  shown  that  the  scattering  rate  for  acoustic  phonon  absorption  and 
emission  is  equal  to 


X 

ac 


(K) 


ra  (2+5  )ir  H2k„TN(E  ) 
max  m,n  B + 


m=l 


tfcT 


(4.41) 


where  H is  the  deformation  potential  and  cL  is  the  elastic  constant 
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associated  with  acoustic  vibration  which  is  equal  to  ps2,  where  p Is  the 

mass  density  and  s the  speed  of  sound.  X (K)  Is  a summation  over  all 

ac 

of  the  allowed  final  K planes.  This  result  (4.41)  is  exact  In  the  long 

k T 

wavelength  approximation.  We  used  u = s-Q,  and  n(0  = JL  rather  than 

'(  Q ps  Q 

the  Bose-Einstein  distribution.  Here,  ra^  is  given  by  the  largest 
integer  that  does  not  exceed  (S  /E  )^2, 


m = Ent(E  /E„) 
max  + 0 

K 


1/2 


(4.42) 


For  optical  phonon  scattering  one  finds: 


m o n(u  )N(E  +’4o)  ) absorption 

* max  (2+5  )mD2  ° K 0 

X (K)  = l 2.  x J 

opt  __  i 2pcj  1 

[n(a)g)+lJN(E+-!iu^)  emission 
K 


. (4.43) 


Here  Ent{[(E_^  ± ^J/Eq]  2/^2}  . These  results  are  exact  in  the 

K 

case  of  long  wavelength  phonons;  no  use  has  been  made  of  the  MCA. 
However,  in  the  case  of  the  polar  optical  interaction,  the  MCA  is  a 
necessity  to  obtain  an  analytical  approximation  [37]: 


m 2 
max  e to. 


x<ib  . r ["<v + i 5 4j 

polopt.  m=l  p 


{- 


1 + 5 


n,m 


[(m-n)4E0  + 2(m-n)2EQ(2E  ± Jiu*)  + (H<o*)2]1/2 


[(m+n)4E2  + 2(m+n)2EQ(2EK  ± flu*)  + (tfuj*)2]1/2 


} 


(4.44) 


2 2, 


where  ± (n^")EQ  and  ER  = tf2(k)2/2m*  . Note  that  Jiu,*  can  be 

negative.  Again, 
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\a*-E"tf[<\i*V/Eo]1/2> 

K 

. 1_  _ 1 1 

n e “e  " e : eoo  13  the  high-f requency  permitivity  and  e the 

P 00  g S 

static  permitivity. 

The  scattering  rates  are  calculated  for  the  following  electron 
phonon  Interactions  that  are  considered:  acoustic  phonon,  equivalent 

and  nonequivalent  intervalley  phonon  interaction  via  the  optical 
deformation  potential,  and  the  polar  optical  interaction.  Although  the 
MCA  is  exact  for  long  wavelengths  in  the  case  of  deformation  potential 
scattering,  it  is  not  quite  correct  for  the  polar  optical  interaction. 
The  error  increases  from  10%  for  qL  = 4 to  60%  for  qL  = 1.  In  GaAs  the 
wave  vector  is  estimated  by 


q = (2 mVj/yl)172  » 2 x 10  V1 


Accordingly,  the  MCA  gives  an  error  ~ 10%  in  the  polar  optical  scatter- 
ing rates  for  well  widths  ~ 200A. 

4-2.3.  Final  wave  vector  after  scattering  event 

After  a real  scattering,  it  is  assumed  that  the  electron  makes 
an  instantaneous  transition  in  k-space  from  initial  wave  vector 
K = ky>  kz)  to  the  final  wave  vector  K'  = (-2^,  k'  , k^)  , i.e., 

intracollision  field  effects  are  neglected  [38].  These  final  wave 
vectors  are  randomly  distributed  and  found  again  by  generating  random 
numbers  with  a uniform  distribution. 

In  a bulk  material,  without  quantization,  the  possible  wave 
vectors  are  distributed  uniformly  over  the  k-space.  However,  for  a 2DEG 
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the  wave  vector  Is  quantized  in  one  direction  (x-directlon) , and  the 
possible  wave  vectors  are  uniformly  distributed  across  parallel  planes, 
perpendicular  to  the  x-direction  (see  Fig.  4.5).  Once  a scattering 
channel  is  chosen,  the  final  energy  of  the  electron,  after  scattering, 
gives  the  length  of  the  final  wave  vector  via 

|K'|  = (2m*E’/tf2)1/2  . (4.45) 


The  body  of  constant  energy  E'  is  a sphere  of  radius  |K' | in  wave  vector 

space.  The  circles  which  are  the  intersections  of  the  sphere  and  the 

wave  vector  planes  form  the  collection  of  allowed  final  wave  vectors 

from  which  one  has  to  be  chosen  (see  Fig.  4.6).  For  the  deformation 

potential  scattering  processes,  acoustic  and  intervalley,  the  choice  is 

straightforward.  The  allowed  final  K'  states  are  uniformly  distributed 

on  the  circles  in  the  n'  planes  (n1  =1,  2,  ...,  n . where  n 

•nax  max 

1/2 

Ent{ (E'/Eq)  })  . Hence,  the  circles  with  largest  circumference,  or 

equivalently  the  largest  radius,  contain  the  largest  number  of  states. 
The  probability  of  choosing  the  plane  should  be  proportional  to  the 
radius  of  the  circle.  The  radius  of  the  circle  is  proportional 

to  sin(Y^()  where  is  the  half  angle  of  the  cone  described  by 
vectors  starting  in  the  origin  and  ending  on  the  circle  of  constant  E' 
in  plane  n'  (see  Fig.  4.6). 


cos (y  ,)  = 
n 


n'  ir/L 


(4.46) 


The  radius  of  the  circles  is  proportional  to  sin(Y  ) 

n' 


sin(Y  ,)  = /l  - cos^(y  ) . 
n n' 


(4.47) 
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kx 


Fig-  4.5.  Distribution  of  allowed  wave  vectors. 
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k'x 


Fig.  4.6.  Definition  of  the  geometry  Involved  In 
obtaining  the  final  wave  vector. 
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The  angle  Y^,  e [ 0 , tt/ 2 ] ; therefore,  sin^,)  > 0.  The  normalized  proba- 
bility for  choosing  plane  k'  = ~ is  then 

x L 

sin(y  ,) 

p(n.)=- 3 . (4.48) 

max 

I 3in(y  () 

n'=l  n 


The  relation  between  the  generated  random  number,  r,  and  the  probability 
for  choosing  plane  n'  is  (see  eq.  4.12) 


r 


I 


n'=l 


m 

max 

sin(y  )/  T sin(y  ,) 

n , , m ' 

m=l 


(4.49) 


which  has  to  be  solved  for  m.  There  is  no  analytical  solution  for  this 
problem;  therefore,  we  follow  the  same  method  that  was  used  to  find  the 
scattering  channels.  Divide  the  interval  [0,1]  into  subintervals  of 
length 

m 

max 

sin(y  )/  l sin( y ,) 
n m * =1 

and  generate  a random  number  0 < r < 1.  When  r falls  in  the  interval  of 
p(n'),  then  plane  n'  is  chosen  and  V is  chosen.  The  other  two  compo- 
nents  of  K'  are  resolved  by  generating  a random  number,  s,  and  finding  <p 
(see  Fig.  4.6)  according  to 


4>  = 2irs  . 
It  follows  that 


(4.50) 


% = /|K,|2  - k;2  cos(*> 

and 


(4.51) 
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V = / |K'|2  - k;2  sin(<»  . 


(4.52) 


For  the  polar  optical  interaction  the  resolution  process  for  the 
three  components  of  K'  is  a little  bit  more  complicated  due  to  the 
|K  - K*  | dependence  of  the  transition  rate  S(K,ic'  ) . Let  Q be  the 
angle  between  the  initial  wave  vector  K (just  before  scattering)  and  the 
final  wave  vector  K'  (just  after  scattering),  see  Fig.  4.7.  The  factor 

I"*"  * i ~2 

|K  K | maximizes  lot  8 = 0;  in  other  words  polar  optical  interaction 
tends  to  favor  scattering  over  small  angles.  Consequently,  polar  opti- 
cal scattering  is  not  randomizing.  The  probability  p(S)  for  an  angle 
8 is  [39] 


p(8) 


sin(B) 


E + S'  - 2/ ( EE ' ) cos  8 


(4.53) 


Accordingly,  cos(8)  can  be  found  with  the  help  of  eq.  (4.12), 


cos(8)  . i-LJ  - (!-*»• 


(4.54) 


where  r is  a generated  random  numsber  and  f equals  2/ EE ' / (/E  - /E5”)2 
The  angle  6 defines  a cone  of  half  angle  8 and  base-to-top  distance 
I K ' |cos8  around  K.  The  circle  at  the  base  is  the  set  of  possible  wave 
vectors  after  scattering.  The  collection  of  points  in  k-space  that  are 
intersections  of  this  circle  with  the  constant  energy  circles  (Fig.  4.6) 
forms  the  solution  for  the  problem  of  finding  the  wave  vector  after 
scattering  with  a polar  optical  phonon.  The  circle  at  the  base  of  the 
cone  (Fig.  4.7)  can  be  described  as  follows: 

U + XV  + yW 


(4.55) 


124 


K' 


Fig.  4.7.  Geometry  for  a polar  optical  scattering  process. 
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where 


-V 


U 


|K' |cos(3)  > 

|K| 


- IsiL  cos(3) 
1*1 


V and  W are  two  normalized  (length  = 1)  vectors  perpendicular  to  U and 
perpendicular  to  each  other.  X and  y are  equal  to  |K' |sin(3)cos(<f>) 
and  |K' |sin(3)sin($),  respectively,  where  0 < <+>  < 2m.  p is  defined 
with  respect  to  any  arbitrary  line  starting  in  the  center  and  inter- 
secting the  circle;  | K * | s in( 0 ) is  the  radius  of  the  circle.  The 
constant  energy  circles  in  the  n'  planes  are  parametrized  as  follows: 


-n'  tt 


(lx- 


2 2 

n'  tt  N 1/2  ...  ,2 

~— ) cosOji),  ( X* 

L 


4A1/2sin(*>) 

h 


(4.56) 


where  ip  is  defined  on  [0,2m]  with  respect  to  some  arbitrary  axis  through 
the  center  of  the  circle,  in  the  plane  of  the  circle. 

Equation  of  the  first  components  of  the  vectors  in  expressions 
(4.55)  and  (4.56)  yields 


U + XV  + yW 

X X X 


n'  tt 

IT  ; 


Ik*  I 

. cos(B)  ^ + | K’  |sin(3)cos(<j>)V  + | K*  | sin(3 )sia($ )W  = ^ . 

I K I X Li 

(4.57) 

From  this  equation  <p  can  be  found,  and  it  has  either  two  solutions  or  no 
solutions  at  all: 


<f>,  = arcsin  ( 


(V2  + U2)1/2 


•)  - a 


$2  = 17  “ a ~ arcsin  ( 


(V2  + W2) 1/2 


) 


(4.58a) 


(4.58b) 
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where 


A 


1 

K1 | sin( S ) 


-tl'TT  _ |K' 


(4.59a) 


and 


sin(a) 


(V2  +U2)in 
X X 


(4.59b) 


For  each  plane,  n'  - 1,  2,  ...,  nmax,  one  finds  two  intersections 

provided  that 


-1  < 


(v2  + y2)1/2 

X X 


< 1 . 


Since  all  of  these  points  are  equally  likely  to  be  final  wave  vectors, 

the  interval  [0,1]  is  divided  into  2*n  subintervals  of  equal  length’ 

when  the  generated  random  number  falls  into  one  of  the  intervals,  we 

choose  the  associated  final  wave  vector  as  the  solution.  This  choice 

automatically  resolves  k'  and  k'  by  analogy  with  eq.  (4.57).  3efore 

y z 

continuing  with  tne  next  section,  in  which  we  present  the  results  of  our 
simulations,  we  would  like  to  discuss  what  we  might  expect  from  our 
model.  The  two-dimensional  density  of  states  is  smaller  than  the  three- 
dimensional  density  of  states  for  a given  energy,  especially  for  smaller 
energies.  In  calculating  the  scattering  rates,  we  performed  an  inte- 
gration over  the  density  of  states.  Hence  we  expect  that  the  scattering 
rates  for  a 2DEG  are  smaller  than  in  a bulk  electron  gas.  The  reduced 
scattering  rates  should  enhance  the  mobility  of  the  charge  carriers. 
The  MCA  has  mainly  an  influence  on  the  results  for  the  polar  optical 

In  the  exact  expression  for  G(q^)  we  see  that 


phonon  interaction. 
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instead  of  delta  functions  5 ±(k,±k  ) os  in  the  bulk  case,  we  find 

sin(i<L/2)/ (kL/2)  which  has  a finite  width  around  = 0.  Accordingly, 
there  are  more  allowed  phonon  wave  vectors  to  cause  the  transition,  thus 
increasing  the  scattering  rate.  The  MCA  limits  the  availability  of 
-inal  states  resulting  in  an  underestimation  of  the  scattering  rates  and 
an  overestimation  of  the  mobility.  However,  for  a well  width  of  200  A 
the  error  should  not  be  larger  than  10%. 

4.3.  Discussion  of  Results 

The  computer  program  used  to  simulate  the  charge  transport  in  a 
2DEG  is  a modification  of  the  program  for  electron  transport  in  bulk 
GaAs  published  by  Boardman  [39].  Apart  from  the  different  scattering 
rates,  different  mechanisms  to  find  the  wave  vector  after  scattering,  as 
discussed  in  the  previous  sections,  had  to  be  implemented.  The  material 
for  which  we  simulated  the  charge  transport  is  GaAs,  a direct  gap  semi- 
conductor with  an  energy  difference  of  1.43  eV  between  the  top  of  the 
valence  band  and  the  bottom  of  the  central  valley  (T  point).  In  the 
band  model  that  has  been  employed  six  satellite  valleys  are  included, 
positioned  at  the  edge  of  the  Brillouin  zone  along  the  main  crystallo- 
graphic axes,  offset  by  0.36  eV. 

The  quantization  of  the  wave  vector  splits  up  the  conduction  band 
in  subbands.  We  only  consider  charge  transport  in  the  conduction  band. 
The  bottom  of  these  parabolic  valleys  is  located  at  n^Eg  (Fig.  4.3). 
The  scattering  mechanisms  considered  are: 

- acoustic  phonon  emission  and  absorption; 

- polar  optical  phonon  emission  and  absorption; 
equivalent  intervalley  phonon  emission  and  absorption; 

- nonequivalent  intervalley  phonon  emission  and  absorption. 
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Fig.  4.8.  Dispersion  relation;  not  on  scale. 
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The  acoustic  phonon  scattering  Is  an  elastic  process;  the  energy 
of  an  electron  does  not  change  in  the  process  (the  phonon  energy  Is 
negligible),  and  the  distribution  of  the  scattering  angle  is  uniform. 
The  acoustic  scattering  occurs  via  the  deformation  potential  in  the  long 
phonon  wavelength  approximation.  Hence,  only  small  phonon  wave  vectors 
are  considered,  and  the  acoustic  phonon  scatterings  involve  intravalley 
transitions.  The  approximation  excludes  the  short  wavelength  or  large 
wave  vector  phonons  necessary  for  an  Intervalley  transition. 

ihe  polar  optical  phonon  scattering  is  inelastic  and  nonrandom- 
izing.  In  the  case  of  emission  of  such  a phonon,  the  electron  loses  an 
energy  (~  30  meV) , while  for  an  absorption  process  the  electron 
energy  is  increased  by  Contrary  to  deformation  potential  scatter- 
ing, there  is  a preference  for  scattering  over  small  angles  in  the  event 
of  a polar  optical  phonon  scattering.  In  the  process  of  such  a scatter- 
ing the  electron  wave  vector  cannot  change  very  much.  Therefore,  only 
electrons  with  wave  vectors  very  close  to  the  point  where  the  central 
and  satellite  valleys  cross  over  (see  dispersion  diagram,  Fig.  4.8)  will 
invoke  an  Intervalley  transition.  The  probability  of  finding  an  elec- 
tron there  is  very  small  and  neglected  in  the  model.  Only  polar  optical 
intravalley  transitions  are  considered. 

The  equivalent  intervalley  phonon  emission  and  absorption  is 
inelastic  and  randomizing.  For  the  transition  from  one  satellite  valley 
to  another  satellite  valley,  a large  phonon  wave  vector  is  needed  to 
accommodate  the  large  change  in  the  electron  wave  vector  which  changes 
from  one  Brillouin  zone  edge  to  another  zone  edge.  The  phonon  branch 
that  also  includes  large  wave  vectors  is  the  optical  branch.  This 
branch  is  approximated  in  the  dispersion  relation  by  a constant  energy 
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*“e  for  all  phonon  wave  vectors.  The  optical  scattering  Is  via  the 
deformation  potential. 

The  nonequivalent  Intervalley  emission  and  absorption  processes 
are  also  optical  phonon  scatterings  via  the  deformation  potential. 
Again  the  processes  are  inelastic  and  randomizing.  The  energy  asso- 
ciated with  a nonequivalent  intervalley  (satellite  valley  <->•  central 

valley)  transition  is  Hu  . 

n 

Simulations  are  reported  here  for  different  values  of  the  non- 
equivalent intervalley  deformation  potential  E using  a well  width  of 
200  A and  compared  with  the  results  of  a computer  simulation  that 

describes  charge  transport  in  bulk  GaAs . In  Fig.  4.9  the  drift  velocity 
is  plotted  as  a function  of  the  electric  field.  As  the  deformation 

potential  is  decreased,  the  low-field  mobility  does  not  change  because 
the  dominating  scattering  processes  are  polar-optical  emission  and 
absorption.  A change  of  E hardly  influences  the  total  scattering  rate. 
Vhen  the  electron  kinetic  energy  becomes  comparable  with  the  energy 
offset  between  the  central  and  satellite  valleys,  intervalley  transfer 
occurs.  At  these  higher  energies,  the  polar-optical  phonon  scattering 
is  no  longer  the  dominant  scattering  mechanism  (see  Fig.  4.10)  because 
intervalley  scatterings  are  much  more  abundant.  When  the  deformation 
potential  Is  reduced,  the  transfer  of  an  electron  back  to  the  central 
valley  is  reduced.  In  short,  the  relative  probability  of  an  electron 
transferring  from  the  central  valley  to  the  satellite  valley  is  only 

slightly  affected  by  a change  in  the  deformation  potential.  But  once  it 
Is  in  the  satellite  valley,  it  stays  there  much  longer  for  smaller  E 
because  equivalent  intervalley  scattering  is  dominant.  Therefore,  the 

velocity  field  curves  are  identical  at  smaller  values  of  the  electric 
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Fig.  4.9.  Tlie  drift  velocity  vs.  electric  field  strength 

for  different  values  of  S in  109  eV/cm.  The  well 
width  is  200  A.  For  comparison  the  curve  for 
bulk  GaAs  is  included. 
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Fig.  4.10.  The  scattering  rates  in  the  T valley  for  the  dif- 
ferent processes  considered  as  a function  of  the 
energy.  1 = polar-optical  emission,  2 = polar- 
optical  absorption,  3 = acoustic  phonon*  4/5  = inter- 
valley emission/absorption  H = 1.0  x 10^  eV/cm, 

6/7  = intervalley  emission/absorption  S = 1.0 
x 10y  eV/cm. 
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rield  strength,  but  the  smaller  coupling  constant  reduces  the  peak 
velocity  and  the  field  strength  at  which  this  peak  is  located.  Finally, 
the  higher  relative  occupancy  of  the  satellite  valleys  for  smaller  H 
reduces  the  saturation  velocity.  The  relative  occupancy  is  plotted  in 
Fig.  4.11  for  the  different  values  of  E. 

The  straightforward  application  of  Ridley's  theory  has  resulted  in 
a higher  low-field  mobility  for  2DE3  as  compared  with  bulk  (E  = l.Q 
10  eV/cm) . This  effect  can  be  attributed  to  the  fact  that  fewer 
states  are  available  to  scatter  to  in  the  case  of  a 2DEG.  Results  for 
the  diffusion  coefficient  are  shown  in  Fig.  4.12.  The  low-field  values 
are  consistent  with  the  Einstein  relation  for  a two-dimensional  .as: 
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where  Ep  denotes  the  Fermi  level.  For  a nondegenerate 
E0  “ aF  >>  kBT>  e(l*  (4.60)  complies  with  the  well-known 


gas,  Ep  <<  Eq  or 
result: 
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(4.61) 


Although  the  shape  of  quantum  wells  in  practical  AlGaAs  devices 
differs  from  the  one  we  have  studied,  the  general  trend  of  obtaining  a 
higher  low-field  mobility  and  diffusion  coefficient  and  a lower  satura- 
tion velocity  for  decreasing  dimensionality  is  probably  generally  valid. 
This  phenomenon  can  be  used  to  enhance  the  device  performance  of,  for 
example,  a high-frequency  FET.  A quantum  well  will  lower  the  value  of 
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Fig.  4.11. 


Relative  occupancy  of  the  satellite  valleys 
for  different  values  of  3 In  10^  eV/cm. 
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Fig.  4.12.  Diffusion  coefficient  as  a function  of  the 

electric  field  strength  for  different  values 
of  H in  l(r  eV/cm. 
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the  parasitic  source  and  drain  resistance  in  which  only  small  electric 
fields  are  preserved.  In  addition,  a Gunn  diode  employing  a 2-0  channel 
will  operate  efficiently  at  lower  bias  voltages  than  a 3-D  device,  since 
the  2-D  negative  differential  mobility  regime  extends  to  lower  field 
strengths. 

In  conclusion,  a selective  use  of  quantum  wells  in  device  fabrica- 
tion can  significantly  enhance  device  performance.  The  Monte  Carlo 
program  records  the  history  of  an  electron  in  K— space.  The  momentum  of 

an  electron  is  tfK  and  the  velocity  of  the  electron  is  . It  was 

m 

suggested  by  Handel  [4]  that  the  1,/f  noise  observed  in  many  physical 
systems  [40-42]  can  oe  explained  i ti  terras  of  velocity  changes  of  the 
particles  considered.  The  particles  suffer  collisions  and  in  the 

process  their  velocities  change.  He  suggested  that  the  Brerasstrahlung 
generated  in  the  process  of  accelerating  and  delerating  the  particles 
gives  rise  to  a 1/f  dependence  of  the  spectral  density  of  the  current 
fluctuations.  The  magnitude  of  the  1/f  noise  is  proportional  to 
(Av)  , where  Av  is  the  vectorial  change  of  velocity  suffered  in  a 
collision. 

In  the  simulation  the  electron  is  followed  from  a scattering  event 
to  a scattering  event,  and  it  seems  straightforward  to  calculate  the 
velocity  change  that  the  electron  experiences  in  the  process.  For  an 
electron  that  stays  in  the  same  conduction  valley,  the  mass  of  the  elec- 
tron does  not  change  in  a scattering  event  (parabolic  bands)  and 

Av  = <*a  - Kb)  (4.62) 

m 

+ -f 

where  Kg , are  the  wave  vector  after  and  before  scattering,  respec- 
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tively.  However,  la  the  case  of  a transition  from  the  central  to  the 
satellite  valley  or  vice  versa  or  for  a scattering  within  a nonparabolic 
band,  the  mass  changes  and  it  is  not  clear  any  more  what  should  be  taken 
for  the  mass  in  eq.  (4.62).  Only  for  very  low  electric  fields,  when 
scattering  only  occurs  within  the  central  valley  and  the  electrons  only 
occupy  states  close  to  the  minimum  where  the  band  can  be  approximated  by 
a parabola,  would  eq.  (4.62)  be  valid.  For  these  field  strengths 
aH’  the  measure  for  the  magnitude  of  the  1/f  noise,  is  calculated 

e 

„ _ 4a  (Av) 

H " 3m  ~2  • (4.63) 

The  superscript  e denotes  ensemble  averaging.  In  the  limit  for  low 
field  strengths,  is  found  to  be  ~ 2 x 10~8,  which  agrees  very  well 
with  the  value  Kousik  [43]  calculates. 


CHAPTER  V 

HIGH-FREQUENCY  NOISE  MEASUREMENTS  OF  AN  ALGAAS/GAAS  GATELESS  MODFET 
5.1.  Introduction 

In  the  preceding  chapter  (IV)  the  computer  simulations  of  the  path 
of  an  electron  In  a quasi-two-dlmensional  wave-vector  space  were 
described.  This  investigation  resulted  in  the  presentation  of  the 
dependence  on  the  electric  field  strength  of  two  important  charge-trans- 
port parameters,  the  drift  velocity  and  the  diffusion  coefficient.  In 
the  present  chapter  experiments  are  described  to  ascertain  the  field  and 
frequency  dependence  of  one  of  these  parameters,  the  diffusion  coeffi- 
cient. Specifically,  high-frequency  (0.5  - 12  GHz)  noise-temperature 
measurements  were  performed  on  a heterostructure  with  a conduction-band 
profile  that  approximates  as  well  as  possible  the  ideal,  infinitely 
high,  square  quantum  well  used  in  the  simulations.  The  diffusion  co- 
efficient will  be  found  from  a combination  of  noise  measurements  and 
impedance  measurements  of  the  structure  under  investigation.  To  perform 
the  high-frequency  noise  measurements,  use  is  made  of  a technique  devel- 
oped by  Gasquet  [44].  By  employing  circulators  with  a sufficiently 
large  bandwidth,  short  pulse-biased  noise  measurements  allow  the  deter- 
mination of  the  noise  temperature  up  into  the  hot-electron  regime. 
After  a description  of  the  device  under  test  in  section  5.3,  we  will 
proceed  in  section  5.4  by  discussing  the  setup  to  measure  the  impedance. 
In  that  section  the  measurements  that  were  performed  are  presented  and 
analyzed.  The  technique  to  measure  high-frequency  noise  of  a pulse- 
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biased  device  is  the  topic  of  section  5.5,  where  measurements  are  pre- 
sented also.  The  results  will  be  discussed  in  section  5.6,  where  the 
dir  fusion  coefficient  is  presented  and  compared  with  the  outcome  of  the 
computer  simulations  of  Chapter  IV. 

5.2.  Theory 

Of  the  three  types  of  noise  encountered  in  spectra  of  the  power 
density  of  current  fluctuations,  1/f  noise,  generation-recombination 
noise,  and  velocity-fluctuation  noise,  the  last  one  is  related  to  the 
velocity  changes  charge  carriers  suffer  due  to  collisions.  The  spectral 
density  of  the  velocity  fluctuations  is  a direct  measure  for  the  diffu- 
sion coefficient: 

Sy(w,F)  = 4D(co,F)  (5.1) 

where  Sy  is  the  spectral  density  of  the  velocity  fluctuations,  D the 
diffusion  coefficient,  u the  angular  frequency,  and  F the  electric  field 
strength.  To  find  the  diffusion  coefficient  from  noise  temperature 
measurements,  the  expression  of  the  spectral  density  of  current  fluctua- 
tions in  terms  of  velocity  fluctuations  is  compared  with  the  well-known 
thermal  noise  expression  in  terms  of  the  noise  temperature  and  the 
impedance  of  the  device. 

The  current  i through  a channel  of  cross-sectional  area  A perpen- 
dicular to  the  current  flow  is 

1 ' ! • (5.2) 

where  q is  the  electron  charge,  L the  length  of  the  conducting  channel, 

N the  total  number  of  carriers,  and  v.  the  velocity  of  the  jth  electron. 


140 


Neglecting  fluctuations  in  the  number  of  charge  carriers,  a fluctuation 
of  the  current,  Ai,  is 


Ai 


(5.3) 


The  spectral  density  of  the  current  fluctuations  S,  is  then 


= «Ll 

•7 

L“ 


N S 


(5.4) 


which  is  expressed  in  terms  of  the  diffusion  coefficient  by  using  eq 
(5.1): 


In  the  derivation  of  eq.  (5.5),  it  is  assumed  that  there  is  no  correla- 
tion between  the  velocity  fluctuations  and  that  the  electric  field 
inside  the  device  is  independent  of  position. 

Another  name  for  velocity-fluctuation  noise  is  thermal  noise.  The 
spectral  density  of  the  current  fluctuations  in  this  context  is  usually 

expressed  in  terms  of  the  noise  temperature  Tn  and  the  device  Impedance 
Z: 


where  kB  is  Boltzmann's  constant.  It  should  be  noted  here  that  Tn  is 
not  a physical  quantity;  it  does  not,  for  instance,  signify  the  electron 
temperature  or  lattice  temperature.  It  merely  is  a measure  for  the 
total  noise  produced  and  should  be  thought  of  as  an  equivalent  of  this 
noise  power.  A comparison  of  eqs.  (5.5)  and  (5.6)  will  yield  an  equa- 
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tion  that  relates  tne  diffusion  coefficient  to  the  noise  temperature: 
L^k_T 

D = -T  ^e[Z_1]  . (5.7) 

N 

It  follows  from  eq.  (5.7)  that  in  order  to  determine  the  diffusion 
coefficient,  measurements  of  the  noise  temperature  and  the  device  impe— 
dance  need  to  be  performed. 

5.3.  Device  Description 

To  simplify  the  analysis,  measurements  were  performed  on  gateless 
AlGaAs/GaAs  MODFETs.  In  order  to  allow  for  a comparison  between  the 
experimental  results  and  the  outcome  of  the  Monte  Carlo  simulations,  the 
structure  would  have  to  exhibit  a square  conduction-band  well  and  an 

electron  density  small  enough  for  the  electron  gas  to  be  nondegenerate. 
The  shape  of  the  quantum  well  strongly  affects  the  density  of  states  as 
a function  of  energy  and  therefore  the  scattering  rates.  The  reason 

that  a nondegenerate  electron  gas  is  preferred  is  that  the  computer 

simulations  do  not  account  for  carrier-carrier  interactions  which  become 
more  important  when  the  carrier  density  increases. 

In  Fig.  5.1  the  composition  of  the  different  layers  that  make  up 
the  heterostructure  is  drawn.  On  an  undoped  GaAs  substrate  a superlat- 
tice of  five  periods  of  a 50  A layer  of  intrinsic  GaAs  alternated  with  a 
50  A layer  of  intrinsic  AlAs  is  grown.  These  superlattices  have  been 
shown  to  enhance  the  quality  of  the  structure  on  top  of  the  superlattice 
since  possible  lattice  errors  in  the  substrate  will  not  propagate.  The 
channels  in  which  the  conduction  of  current  takes  place  are  the  250  A- 
wide  intrinsic  GaAs  layers.  On  each  side  of  these  layers  100  A-thick 
intrinsic  AlQ.3Ga0.7As  buffer  layers  have  been  grown,  followed  by 
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Fig.  5.1.  Layer  configuration  of  the  gateless  MODFET. 
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175  A -thick  intentionally  doped  n-type  A1Q  ^Gag^As  layers.  The  doping 
concentration  results  in  1018  electrons  per  cm3  in  these  layers. 

The  rather  thick  buffer  layers  prevent  an  effective  transfer  of 


electrons  from  the  doped  AlGaAs  region  to  the  GaAs  well.  This  helps  to 

tteep  the  electron  concentration  in  the  well  low.  From  a measurement  of 

the  dc  source  to  drain  resistance  the  estimated  electron  concentration 
16  3 

is  3 x io  cm  , resulting  in  a sheet  carrier  concentration  of  7.5  x 

10  9 

10  cm  . As  a result  of  the  low  carrier  concentration  in  the  well, 
the  charge  on  either  side  of  the  AlGaAs/GaAs  interface  will  be  small, 
and  hence  the  electric  field  normal  to  the  interface  caused  by  these 
charges  will  be  relatively  minor.  The  bending  of  the  bends  will  there- 
fore be  insignificant,  guaranteeing  an  almost  perfectly  shaped  square 
well.  The  conduction-band  structure  resulting  from  a first-order 

analysis  is  shown  in  Fig.  5.2.  The  large  41  mole  fraction  (0.3)  results 
in  a conduction-band  discontinuity  of  about  200  raeV  [7]. 

The  drain-to-source  spacing  is  6 pm  and  the  width  is  100  um. 

5.4.  Impedance  Measurement s 

As  a first  step  to  characterize  the  impedance  of  the  device, 
pulsed  current-voltage  (I-V)  characteristics  were  measured.  The  setup 
is  depicted  in  Fig.  5.3.  The  applied  pulses  were  typically  60  ps  long 
and  the  repetition  rate  was  900  Hz.  The  desktop  computer,  a Tektronix 
4054A,  drives  the  Tektronix  7603  oscilloscope  with  the  programmable 
digitizer  7D20  which  measures  the  two  voltages.  From  the  difference  of 

the  two  voltages  the  current  through  the  device  is  calculated.  These 


computations  are  performed  by  the  computer,  which  also  calculates  the 
derivative  of  the  I-V  curve.  The  output  of  the  printer  is  a list  of 
current  and  voltage  points  and  dV/dl.  Figure  5.4  shows  the  results  of 
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Fig.  5.2.  Schematic  band  diagram  (not  to  scale). 
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Fig.  5.3. 


Pulsed  I-V 


measurement  setup. 
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these  experiments  for  pulses  of  both  polarities.  Both  curves  show  a 
resistance  of  220  £2  in  the  ohmic  regime.  For  voltages  larger  than  1 V 
(F  - 1.7  kV/cm),  the  current  increases  sublinearly  with  Increasing  bias. 
This  can  be  understood  as  follows.  For  large  enough  electric  field 
strengths  the  electrons  gain  so  much  kinetic  energy  that  it  becomes 
comparable  to  the  energy  offset  between  the  central  valley,  where  the 
electrons  originally  reside,  and  the  satellite  valley.  The  scattering 
of  nonequivalent  intervalley  phonons  with  these  hot  electrons  will 
result  in  the  transfer  of  these  electrons  to  the  satellite  valleys. 
Since  the  effective  mass  is  larger  in  the  satellite  valley,  the  mobility 
is  reduced,  which  in  some  cases  might  result  in  a maximum  in  the 
velocity-field  curve.  In  bulk  GaAs  the  maximum  velocity  occurs  around 
field  strengths  of  3 kV/cm.  This  electric  field  value  is  also  indicated 
in  Fig.  5.4.  The  derivative  dl/dV  for  both  I-V  characteristics  is 
plotted  in  Fig.  5.5.  For  voltages  larger  than  0.7  V,  or  electric  field 
strengths  larger  than  1.2  kV/cm,  the  differential  conductance  decreases 
very  rapidly.  Note  that  in  this  analysis  the  electric  field  strength  is 
assumed  to  be  independent  of  position  and  equal  to  V/L,  i.e.,  space- 
charge  injection  from  the  source  and  drain  contacts  is  neglected.  That 
this  assumption  seems  to  be  warranted  follows  from  the  next  argument. 
The  onset  of  space-charge  injection  occurs  when  the  transition  time  t 
of  the  electron  becomes  of  the  order  of  the  ohmic  relaxation  time  t^ : 

r = k 

tr  v (5.8) 

and 


t 


pe  . 


(5.9) 
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Fig.  5.4.  Pulsed  current -voltage  characteristics  for 

both  polarities.  The  pulse  duration  Is  60  us 
and  the  repetition  rate  Is  900  Hz. 


(tXJ)  AP/IP 
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Fig.  5.5.  The  derivative  of  the  I-V 


curves  In  Fig.  5.4. 
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It  follows  from  these  two  equations  that  the  Injection  sets  on  for  vol- 
tages 

? 

V = 

e * (5.10) 

For  the  device  under  consideration,  the  value  for  V is  of  the  order  of 
2 kV. 

tub  impedance  of  che  hecerostruccure  device  u»s  also  measured  as  a 
function  of  frequency  for  aero  bias.  For  this  purpose  Che  reflection 
coefficient  of  the  device,  mounted  on  a stripline,  was  measured  for  fre- 
quencies between  100  MHz  and  2 GHz.  The  setup  for  these  measurements  Is 
shown  in  Fig.  5.6.  The  network  analyzer  measures  the  magnitude  and  the 
phase  angle  of  the  reflection  coefficient  as  a function  of  frequency. 

The  reflection  coefficient  is  a function  of  the  Impedance  that  cermln- 
ates  the  stripline: 


In  eq.  (5.11)  ^ is  the  reflection  coefficient,  7^  the  impedance  that 

terminates  the  stripline,  and  Zq  the  characteristic  impedance  of  the 

stripline  (50  ft).  It  is  assumed  here,  and  also  experimentally  verified, 

that  the  stripline  is  virtually  lossless  « 0.1  dB)  for  all  frequencies 

of  interest.  Therefore,  the  measured  reflection  coefficient  T has 

exp 

the  same  magnitude  as  T The  phase  angle  of  T (/  r ),  however 

exp  — exp  C| 

differs  by  2 8£  from  the  phase  angle  of  where  8 is  the  propagation 
constant  (2rf/c),  f is  the  frequency,  c the  velocity  of  light,  and  £ the 
electrical  length  the  signal  travels  between  the  reference  plane  and  the 
Since  the  reference  plane  was  chosen  at  the  SMA 


device. 
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Fig.  5.6.  Experimental  setup  to  measure  the  reflection 
coefficient  of  the  device  under  test  (DUT). 


corresponds  co  the  length  of  the 


the  sample  mount,  l approximately 
stripline. 


2 34 


(5.12a) 


(5.12b) 


The  real  aad  Imaginary  part  oE  the  teralaatlng  Impedance  ^ can  be  found 

from  T : 

T 


(5.13a) 


(5.13b) 


But  before  the  real  and  imaginary  part  of  the  terminating  impedance  can 
be  found,  the  electrical  length  has  to  be  determined.  This  is  done  by 
plotting  the  phase  angle  of  the  reflection  coefficient  versus  the  fre- 
quency. For  frequencies  of  the  order  of  1 GHz  and  electrical  lengths  of 
the  order  of  5 cm,  one  calculates  2 84  - 2 rad  = 120°.  As  the  frequency 
increases  one  sees  from  eq.  (5.12a)  that  the  phase  angle  will  be  totally 
determined  by  the  path  length  1 the  signal  has  to  travel.  In  the  high- 
frequency  limit  one  expects  a linear  relationship  between  the  measured 
phase  angle  / and  the  frequency  f.  From  the  slope  the  electrical 

length  is  computed.  This  procedure  has  been  performed  on  the  measured 
reflection  coefficient  data.  In  Fig.  5.7  the  phase  angle  is  plotted  as 

a function  of  the  frequency  and  the  data  are  fitted  to  a linear-least- 
squares  curve  through  the  origin  of  the 
is  found  to  be  7.500  ± 0.017  cm. 


axes.  The  electrical  length  l 
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Fig.  5.7.  Measured  phase  angle  of  the  reflection  coeffi- 
cient as  a function  of  the  frequency.  The  solid 
line  Is  the  linear-least-squares  fit  to  the  data 
that  Intersects  with  the  origins  of  the  axes. 
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lJow  chac  the  value  of  the  electrical  length  has  been  established, 
It  Is  used  to  calculate  j_  1^,  and  with  the  help  of  eqs.  (5.13a)  and 
(5.13b)  the  real  and  imaginary  part  of  the  impedance  %T  can  be  evalu- 
ated. The  results  of  these  computations  are  shown  in  Figs.  5.3  and  5.9. 

The  impedance  found  in  this  process  can  be  approximated  by  an 
equivalent  circuit.  This  equivalent  circuit  consists  of  a 100  pF  capac- 
itor, mounted  to  provide  an  ac  path  to  ground  in  series  with  the  device 
which  is  basically  a resistor  with  a possible  geometrical  capacitance  in 
parallel.  The  bonding  wire  that  connects  the  stripline  and  the  source 

of  the  FET  is  modeled  as  an  inductor.  This  equivalent  circuit  is  drawn 
in  Fig.  5.10. 

The  impedance  of  the  equivalent  circuit  is  split  in  its  real  and 
imaginary  parts.  The  magnitude  of  these  two  components  is  simulated  as 
a function  of  the  frequency  and  plotted  in  Figs.  5.3  and  5.9,  respec- 
tively (solid  line).  The  values  for  the  parasitic  capacitance  and 
Inductance  are  reasonable.  The  capacitance  can  be  associated  with  the 
geometry  of  the  device.  The  conductance  is  associated  with  the  bonding 
wire,  which  in  our  case  was  ~ 0.5  cm. 

It  should  be  noted  that  although  the  error  in  the  calculated  elec- 
trical length  is  of  the  order  of  0.2%,  the  subtraction  of  two  large 
numbers,  [_  and  2 S£,  gives  rise  to  a large  relative  error  in  the 
answer,  [_  iy  Especially  for  higher  frequencies,  this  error  will  become 
large.  Hence,  the  evaluation  of  the  impedance  will  be  less  accurate  as 
the  frequency  increases.  For  this  reason  efforts  to  evaluate  the  impe- 
dance for  frequencies  larger  than  2 GHz  were  not  undertaken. 

In  the  previous  section  it  was  shown  that  in  order  to  calculate 
the  diffusion  coefficient,  the  real  part  of  the  admittance  had  to  be 
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Fig.  5.8.  The  real  part  of  the  Impedance  ZT  as  a 

function  of  the  frequency.  The  solid  line 
represents  the  equivalent  circuit  that 
describes  the  impedance  (see  text). 
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Fig.  5.9.  The  imaginary  part  of  the  Impedance  as  a 

function  of  the  frequency.  The  solid  line 
represents  the  equivalent  circuit  that 
describes  the  impedance  (see  text). 
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Fig.  5.10.  Equivalent  circuit  representing  the  impedance 

that  terminates  the  stripline.  The  two  elements 
R and  C2,  in  the  dashed  box  model  the  device. 

L - 15  nH,  R = 220  12,  = 100  pF,  and  C2  = 0.3  pF. 
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found.  In  Fig.  5.11  the  measured  real  part  of  1/ZT  Is  plotted  versus 
frequency.  The  solid  line  represents  the  real  part  of  1/ZT,  as  It  Is 
calculated  from  the  equivalent  circuit  of  Fig.  5.10. 

5-5*  Noise  Temperature  Measurements 

Measurements  of  the  noise  temperature  were  performed  as  a function 
of  the  electric  field  strength  (F  < 3 kV/cm)  and  as  a function  of  fre- 
quency (0.5  GHz  < f < 12.0  GHz).  In  order  to  be  able  to  measure  the 
noise  at  these  field  strengths,  the  bias  needs  to  be  applied  in  short 
pulses  to  avoid  any  effects  due  to  Joule  heating  of  the  device.  In  this 
range  of  field  strengths  the  device  is  operated  in  the  nonohmic  regime, 
l.c.,  its  impedance  changes  when  the  bias  is  switched  on  and  off.  Con- 
sequently, the  amplifier  will  see  different  impedances  when  the  bias  is 
varied.  As  a result  the  noise  signal  contribution  of  the  amplifier  to 
the  noise  signal  at  the  output  of  the  measuring  setup  is  different  for 
the  different  biases.  To  circumvent  this  problem,  Gasquet  et  al.  [44] 
designed  the  so-called  circulator  method.  In  Fig.  5.12  the  setup,  as  it 
was  employed  in  the  measurements  presented  here,  is  shown. 

The  first  pulse  generator  (PG1),  an  HP  214A,  supplies  100  ys- 
long  pulses  at  a rate  of  50-100  Hz  to  the  device.  The  capacitor 
(100  PF)  provides  an  ac  path  to  ground,  making  it  possible  to  measure 
the  noise.  The  circulator,  a three-port  device,  conducts  the  signal 
only  in  one  direction.  A noise  signal  generated  by  the  noise  source 
(NS),  HP  346B , will  travel  through  the  stripline  on  which  the  device  is 
mounted,  be  reflected  at  the  stripline-device  interface,  and  then  travel 
the  circulator  where  it  is  passed  on  to  the  bias  tee,  HP11590A, 
and  into  the  low-noise  amplifier  (LNA) . The  big  advantage  of  using  the 
circulator  is  that  independent  of  the  settings  of  switches  SI  and  S2, 
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Fig.  5.11.  Calculated  values  of  RefZ^1]  from  reflection 
coefficient  measurements.  The  solid  line 
represents  RefZ^1]  according  to  the  equivalent 
circuit  in  Fig.  5.10. 
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Fig.  5.12.  Noise  measurement  setup. 
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and  independent  of  whether  the  bias  is  on  or  off,  the  amplifier  will 
always  see  the  Impedance  of  the  noise  source  (NS),  which  is  50  n. 

The  bias  tee  provides  a dc  path  to  ground  for  the  bias  pulses. 
The  signal  is  fed  into  a low-noise  amplifier,  then  passed  on  to  the 
spectrum  analyzer  (HP  S53A).  The  spectrum  analyzer  mixes  the  signal 

down  to  21  MHz,  at  which  frequency  it  can  be  detected  by  the  power  meter 
(HP  8484A) . 

Since  the  bias  is  supplied  in  pulses,  the  detection  of  the  noise 
power  has  to  be  synchronized  with  these  pulses.  For  this  purpose  a 
second  pulse  generator,  PG2,  an  HP  222A,  triggered  by  the  first  pulse 
generator,  drives  a PIN  diode  switch  which  is  closed  only  when  the  bias 
-s  applied.  The  pulses  of  the  second  pulse  generator  were  usually  40  us 
long  and  delayed  with  respect  to  the  rising  edge  of  the  biasing  pulse. 

Therefore,  risetime  effects  at  the  beginning  of  the  biasing  pulse  could 
be  avoided  during  detection. 

In  the  setup  discussed  here,  there  are  four  unknowns:  the  noise 
temperature  of  the  device  (Tn>,  the  gain-bandwidth  product  of  the  con- 
figuration (GB),  the  noise  temperature  of  the  low-noise  amplifier  (T  ), 
and  the  impedance  of  the  device.  The  impedance  shows  uP  in  the  reflec- 
tion coefficient  T.  To  solve  for  this  system  of  unknowns  four  measure- 
ments have  to  be  performed.  These  are  just  the  four  different  combina- 
tions of  settings  of  the  switches  SI  and  S2.  Switch  SI  shorts  out  the 
device  and  switch  S2  switches  between  noise  source  ON  and  noise  source 
OFF,  resulting  in  two  different  noise  signals,  Th  or  T.,  that  are  fed 
into  the  circulator.  The  four  detected 


power  signals  are 
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m = GB(T  + T ) 
a Ac' 

M2  = GB(T  + T) 
a h' 

M3  = GB(Ta  + I F | 2 T + (1  - | r | 2 )T  ) 

n 

M4  = GB(Ta  + |r|2Th  + (1  - !r|2)Tn)  . 


(5.14b) 
(5.14c) 
( 5 . 1 4d  ) 


For  measurements  Ml  and  M2,  the  position  of  switch  SI  Is  such  that  the 
device  Is  shorted  out.  For  measurements  Ml  and  M3,  switch  S2  Is  open 
and  the  noise  source  is  off  (T  = 300  K) . 

From  eqs . (5.14a  - 5.14d),  the  square  of  the  absolute  value  of  the 
reflection  coefficient  and  the  noise  temperature  can  be  calculated: 


M4  - M3 
M2  - Ml 


T,  - T 


T - T = 
n c 


1 - Ir 


c M3  - Ml 
2 M2  - Ml 


(5.15a) 


(5.15b) 


The  reflection  coefficient  Is  measured  as  the  absolute  value  squared, 
and  therefore  an  evaluation  of  the  real  part  of  the  Impedance  Is  only 
possible  when  the  Imaginary  part  Is  negligible  with  respect  to  the  real 
part.  In  that  case  the  real  part  of  the  admittance  Is  equal  to  the 
Inverse  of  the  real  part  of  the  Impedance. 

The  frequency  rang,  of  operation  of  the  bias  tee  Is  from  1 - 12.4 
GHZ.  For  frequencies  below  1 GHz,  the  bias  tee  was  replaced  by  a tuner 
between  the  device  and  switch  SI  to  provide  a dc  path  to  ground,  with 
the  tuner  the  Impedance  of  the  device  was  matched  as  well  as  possible  to 
the  50  0 of  the  system.  This  Improved  the  accuracy  of  the  measured 
noise  temperature.  However,  all  Information  concerning  the  Impedance  of 
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the  device  Is  lost  in  the  matching  process.  Therefore,  the  reflection 
coefficient  cannot  be  used  to  find  the  impedance  of  the  device  for 
frequencies  smaller  than  1 GHz. 


In  the  frequency  range  of  0.5  - 1.0  GHz,  a P & H Laboratories  cir- 
culator model  31-436310  was  used.  The  insertion  loss  between  the  ports 
Is  0.5  dB  and  the  isolation  between  the  ports  ranges  from  18  to  24  d3 . 
The  low-noise  amplifier  used  is  a Trontech  W1G211. 

In  the  frequency  range  of  2 - 4 GHz,  a P & H Laboratories  circu- 


lator model  31-S26322  was  used.  The  insertion  loss  varies  between  0.25 
and  0.40  dB , and  the  isolation  ranges  from  24  to  30  d3.  The  insertion 
loss  is  lowest  at  the  center  of  the  frequency  band  where  the  isolation 
Is  largest,  also.  Because  of  the  wide  band  (2-4  GHz)  of  the  amplifier 
(Amplica,  Inc,  ASM625301),  filters  were  used  to  avoid  overloading  of 
the  amplifier.  The  filters,  K & L Microwave,  Inc.,  model  53120,  have  a 
bandwidth  of  0.2  GHz  and  center  frequencies  2.2,  2.6,  and  3.4  GHz. 

For  the  highest  frequency  range  (8  - 12.4  GHz),  a P i H Labora- 
tories model  31-X16311  circulator  was  used  in  combination  with  a tunable 
filter,  K & L 3FCT  8000/124000,  of  bandwidth  0.1  GHz.  The  amplifier 
used  is  an  Amplica,  Inc.,  AXM545304. 

The  results  of  the  measurements  of  the  noise  temperature  are  dis- 
played in  Fig.  5.13.  The  noise  is  white  up  to  2.2  GHz,  and  for  frequen- 
cies larger  than  2.2  GHz  the  noise  temperature  decays  slowly.  The 
measurements  seem  to  indicate  a more  rapid  decrease  for 
frequencies  exceeding  8 GHz. 

The  data  do  not  seem  to  fit  a 1/f  spectrum,  but  might  be  fitted  to 
one  or  two  Lorentzian-shaped  spectra.  As  opposed  to  generation-recom- 
bination noise  processes,  it  is  not  expected  at  these  high  frequencies 
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Fig.  5.13.  Noise  temperature  as  a function  of  the 

frequency.  The  error  In  each  measurement 
Is  10%.  The  solid  lines  are  Lorentzians 
3dB  = 5 GHz. 
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that  trapping  and  detrapping  of  carriers  will  give  rise  to  the 
Lorentzians.  There  are,  however,  three  other  effects  that  would  give 
rise  to  Lorentzian-shaped  spectra:  transit-time  effects,  real-space- 

charge  transfer,  and  wave-vector-space-charge  transfer. 

For  a device  length  of  6 pm  and  a peak  velocity  of  2 x 105  m/s  at 
the  largest  field  strength,  the  transit  time  is  30  ps.  For  JFETs , 
MESFETs  and  MOSFETs  a maximum  operating  frequency  is  defined  as  (2ttt)_1 
[45],  where  t is  the  transit  time.  This  would  result  in  a maximum  fre- 
quency of  j GHz,  which  is  right  in  the  frequency  range  where  the  noise 
temperature  roll-off  is  observed.  It  should  be  noted  that  transit-time 
effects  can  only  be  observed  when  the  injection  or  extraction  of  a 
single  carrier  gives  rise  to  a velocity  fluctuation  that  is  significant 
with  respect  to  fluctuations  caused  by  scatterings.  For  a carrier  that 
suffers  only  a few  collisions,  the  influence  of  the  extraction  or  injec- 
tion will  be  larger  than  for  a carrier  that  suffers  many  collisions. 
For  a bulk  scattering  rate  of  1013  s"1  [46],  the  electron  would  suffer 
300  scattering  events  per  transit  from  source  to  drain.  From  the  simu- 
lations discussed  in  Chapter  IV,  it  follows  that  in  the  case  of  a 2DEG 
the  scattering  rates  are  smaller,  enhancing  the  transit-time  effect. 
For  electric  field  strengths  smaller  than  2.5  kV/cm,  the  electron 
travels  with  a velocity  roughly  proportional  to  the  electric  field 
strength.  Therefore,  the  maximum  frequency  would  have  to  decrease  with 

the  field  strength.  This  effect  has  not  been  observed  in  the  experi- 
mental data. 

The  two  other  effects  that  would  give  rise  to  a Lorentzian  in  the 
spectrum  are  real-space-charge  transfer  and  wave-vector-space-charge 
transfer.  In  the  former  case  the  electrons  are  transferred  out  of  the 
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GaAs  well  Into  the  AlGaAs . In  the  latter  case  the  electron  stays  In  the 
GaAs  well,  but  transfers  from  the  central  conduction-hand  valley  into  a 
satellite  valley.  For  both  of  these  processes  the  electron  needs  to 

acquire  a minimum  kinetic  energy  before  it  can  transfer.  Hence,  it  is 
expected  that  a critical  field  strength  for  which  these  transfer 

processes  start  to  occur  would  be  observed.  This  is  not  the  case. 

The  roll-off  of  the  noise  temperature  cannot  be  attributed  to  the 
parasitic  elements  of  the  equivalent  circuit.  For  lossless  parasitics 
it  is  shown  in  Appendix  B and  ref.  [47]  that  the  measured  noise  power  is 
not  affected.  It  was  experimentally  verified  that  both  the  100  pF 
capacitor  and  the  inductance  introduced  by  the  bonding  wire  behaved 

ideally  for  frequencies  up  to  8 GHz. 

In  Fig.  5.14  the  measured  noise  temperatures  are  plotted  as  a 
function  of  the  electric  field  strength  for  the  different  frequency 
bands.  The  data  for  the  frequency  band  0.5  - 1.0  GHz  are  also  plotted 
in  Fig.  5.15.  In  this  figure  the  experimental  results  are  compared  with 
experimental  results  [48]  for  bulk  GaAs  and  results  from  computer  simu- 
lations of  the  2DEG  for  different  deformation  potential  values  and  of 
bulk  charge  transport  in  GaAs.  The  noise  temperature  calculated  from 
the  simulations  in  a 2DEG  agree  well  with  the  experimentally  obtained 
data.  It  is  clear  that  the  onset  for  the  increase  of  the  noise  tempera- 
ture occurs  at  smaller  field  strengths  in  the  case  of  the  2D  EG  as  com- 
pared with  bulk  GaAs.  Both  the  measurements  and  the  computer  simula- 
tions indicate  this.  The  occurrence  of  this  hot  electron  effect  can  be 
explained  by  the  higher  mobility  found  from  the  simulations.  Due  to  the 
higher  mobility  the  electron  gas  heats  uP  faster,  and  therefore  hot 
electron  effects  will  occur  for  smaller  field  strengths. 
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Fig.  5.14.  Noise  temperature  as  a function  of  the 

electric  field  strength.  The  numbers  along 
the  curves  refer  to  the  frequency  bands. 
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Fig.  5.15.  Noise  temperature  as  a function  of  the 

electric  field.  The  solid  lines  are  simula- 
tions for  different  values  of  S In  lCr  eV/cm 
The  curves  are  explained  In  the  text. 
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5-6.  Diffusion  Coefficient 

The  ultimate  goal  of  the  experiments  described  In  the  previous 
sections  is  to  ascertain  the  electric  field-strength  dependence  of  the 
diffusion  coefficient  and  if  possible  its  frequency  dependence.  Under 
the  assumption  that  the  measured  noise  is  due  to  velocity  fluctuations, 
the  diffusion  coefficient  is  calculated.  The  low-frequency  value  of  the 
noise  temperature  (Fig.  5.13)  is  multiplied  by  dl/dV  (Fig.  5.5)  for  each 
field  strength.  The  result,  a normalized  plot  of  the  diffusion  coeffi- 
cient as  a function  of  the  electric  field  strength,  is  presented  in  Fig. 
5.16.  For  larger  values  of  the  electric  field  strength,  the  error  in 
dl/dV  increases  and  hence  the  uncertainty  in  the  diffusion  coefficient. 
This  is  indicated  by  the  error  bars  in  Fig.  5.16.  The  most  significant 
feature  of  this  graph  is  the  increase  of  the  diffusion  coefficient  as 
Che  field  strength  increases.  The  increase  of  the  diffusion  coefficient 
has  been  observed  in  bulk  GaAs  [49]  and  is  attributed  to  the  effect  of 
polar  runaway  [50].  Also,  in  Monte  Carlo  simulations  of  charge  trans- 
port in  both  bulk  [51,52]  and  2DEG  GaAs,  an  increase  of  D is  found  when 
the  deformation  potential,  E , for  nonequivalent  intervalley  scattering 
is  reduced  to  0.18  x 10*  eV/cm.  These  results  are  shown  in  Fig.  5.16 
for  both  S - 1.0  x 109  eV/cm,  0.3  x 109  eV/cm,  and  = = 0.13  x i09 
eV/cm.  The  latter  values  of  the  deformation  potential  shows  the  same 
trend  for  the  diffusion  coefficient  as  the  experiments,  except  for  the 

onset  of  the  increase  of  D.  The  optimum  value  for  S will  be  between 
0.13  and  0.30  x iq9  eV/cm. 

In  Fig.  5.16  diffusion-coefficient  data  are  included  as  they  were 
reported  by  Gasquet  [48].  These  experimental  results  for 


bulk  GaAs  show 


NORMALIZED  DIFFUSION  COEFFICIENT 
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ELECTRIC  FIELD  (kV/cm) 


Fig.  5.16.  The  diffusion  coefficient  normalized  to  Its  low- 

flj1jTTalUe*  The  l°w-f re<iuency  noise  temperatures 
and  dl/dV  are  used  to  calculate  D(F)/D(0).  The 
solid  lines  are  Monte  Carlo  simulations  for  dif- 
ferent values  of  S. 
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an  Increase  In  the  diffusion  coefficient.  Hie  occurrence  of  a peak  in 
the  diffusion  coefficient  of  bulk  GaAs  is  associated  with  polar  runaway 
followed  by  wave-vector-space-charge  transfer.  The  onset  for  the  in- 
crease of  the  diffusion  coefficient  of  bulk  GaAs  occurs  for  F = 2 kV/cm 
(Fig.  2.16),  which  is  larger  than  the  field  strength  (0.5  kV/cm)  for 

which  the  increase  of  the  experimental  2 DEG  diffusion  coefficient  hap- 
pens . 

Hess  [54]  explains  that  in  heterostructure  devices  real-space- 
charge  transfer  shows  up,  in  general,  at  electric  field  strengths 
smaller  than  the  field  for  which  wave-vector-space-charge  transfer  would 
occur.  Since  the  presented  experimental  2DEG  diffusion  coefficient 
shows  an  increase  at  field  strengths  smaller  than  for  the  bulk  case, 
real-space-charge-transfer  effects  might  be  present  and  responsible  for 
the  increase  in  D.  The  effect  of  real-space-charge  transfer  is  not 

included  in  the  Monte  Carlo  simulations  since  an  infinitely  high  well 
was  considered. 

The  critical  parameter  in  this  entire  discussion  is  the  low-field 
mobility.  The  simulations  presented  in  Chapter  IV  indicate  a mobility 
which  is  higher  than  the  mobility  that  is  measured  in  heterostructure 
devices  exhibiting  a 2DEG.  Usually,  the  2 DEG  mobility  is  measured  to  be 
8000  cm2 /Vs  [55],  as  compared  with  11000  cm2/Vs  found  from  the  simula- 
tions. As  a result  of  the  high  mobility,  the  hot-electron  effects  occur 
at  smaller  field  strengths,  and  the  increase  of  both  the  noise  tempera- 
ture and  the  diffusion  coefficient  sets  on  at  lower  fields.  The 
measurements  of  the  noise  temperature  as  a function  of  the  electric 
field  strength  are  well  explained  by  the  simulations,  and  it  is  shown 
that  the  noise  temperature  of  a 2 DEG  increases 


more  rapidly  than  the 


noise  temperature  of  bulk  GaAs . The  accuracy  of  the  noise  measurements 
is  within  10%.  The  experimental  diffusion  coefficient  data  are  derived 
from  the  noise  measurements  and  from  dl/dV  because  it  is  assumed  that 
dl/dV  scales  with  the  differential  mobility  u ' = dvd/dF  [53]  in 


This  assumption  follows  from  neglecting  any  type  of  space  charge  in  the 
devices.  Only  a careful  deconvolution  of  the  I-V  curves  can  determine 
the  contribution  of  space-charge  injection  to  the  current.  This  is 
another  reason  why  the  determination  of  the  diffusion  coefficient  has  a 
relatively  large  uncertainty.  An  accurate  measurement  of  the  mobility 
will  be  crucial  to  verify  this  analysis. 

The  higher  mobility  found  in  the  simulations  might  be  explained  by 
the  fact  that  several  scattering  mechanisms  were  not  included.  Possible 
candidates  to  reduce  the  mobility  are  alloy  scattering,  interface  scat- 
tering, or  significant  real-space-charge  transfer.  Less  likely  are 
electron-electron  interactions  and  scattering  with  a two-dimensional 
phonon  gas  instead  of  a three-dimensional  phonon  gas.  Since  the  elec- 
tron concentration  is  small,  electron-electron  interactions  are  insigni- 
ficant. Riddoch  and  Ridley  [55]  showed  that  the  influence  of  a two- 
dimensional  phonon  gas  on  the  scattering  rates  is  maximal  10%.  When  all 
of  these  effects  are  included  and  the  mobility  is  reduced,  it  is  to  be 
expected  that  the  simulations  of  the  noise  temperature  and  diffusion 
coefficient  will  explain  the  presented  measured  data  better. 


CHAPTER  VI 

CONCLUSIONS  AND  SUGGESTIONS  FOR  FURTHER  RESEARCH 

6.1.  1/f  Noise  In  Pseudomorphlc  UODFETs 

Measurements  of  the  power  spectral  density  of  the  current  fluctua- 
tions in  AlGaAs/ InGaAs/GaAs  MODFETs  at  300  K showed  a strong  dependence 
of  the  shape  of  the  spectra  on  the  applied  gate  bias.  For  a gate  vol- 
tage of  -94  mV,  perfect  1/f  spectra  were  measured  for  over  four  decades 
of  frequency.  When  the  gate  voltage  was  either  increased  or  decreased, 
additional  spectral  components  were  observed.  These  components  were 
associated  with  generation-recombination  (g-r)  processes.  Two  interpre- 
tations are  presented,  one  in  which  a resonance  in  the  1/f  noise  causes 
the  anomaly  in  the  measured  spectra  and  one  in  which  the  g-r  noise  in- 
creases when  one  of  the  quantization  levels  in  the  InGaAs  well  is 
reached  by  the  Fermi  level. 

It  seems  worthwhile  to  extend  the  measurements  to  lower  tempera- 
tures. The  Fermi  level  will  go  up  even  higher  in  the  conduction  band, 
and  higher  quantization  levels  can  be  probed.  As  a result,  evidence  can 
be  collected  either  to  support  or  to  refute  the  Interpretation  presented 
in  Chapter  II. 

6.2.  Generation-Recombination  Noise  in  Bulk  Silicon 

It  has  been  shown  that  the  separation  of  a measured  g-r  noise 
spectrum  in  its  constituting  Lorentzians  is  almost  always  warranted. 
This  knowledge  then  was  applied  to  g-r  noise  measurements  performed  as  a 
function  of  temperature  on  planar  silicon  resistors.  The  separated 
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notse  components  were  used  to  find  the  relevant  parameters  that  describe 

the  8lec”lcal  °£  «>e  localised  energy  levels  that  act  as  trap- 

ping centers.  The  acquired  trap  data  were  used  as  Input  of  an  erect 
multitrap -level  computer  program  to  simulate  the  noise  spectra.  The 
zero-frequency  plateau  levels  and  characteristic  times  resulting  from 
the  simulation  are  In  excellent  agreement  with  the  measured  data.  The 
extraction  technique  described  here  Is  shown  to  be  a new  way.  in  addi- 
tion to  DLTS,  to  obtain  Information  pertaining  to  Impurity  traps.  It  Is 
known  that  the  detection  limit  of  DLTS  Is  10’*  Impurities  per  Intended 

dopant.  The  detection  limit  of  the  g-r  noise  method  Is  unknown.  Perhaps 

future  research  in  this  area  will  ascen-^n  m, 

ascertain  the  sensitivity  of  this 

method. 

6*3*  Monte  Carlo  Simulations 

The  simulations  of  ao  electron  traveling  through  the  wave-vectot 

apace  of  a two-dimensional  electron  gas  (2DEG)  resulted  lu  drift 

velocity-field  curves  and  the  electrical  r . 

electrical  field  strength  dependence  of 

the  diffusion  coefficient.  The  presented  mobility  is  higher  than  Is 
found  experimentally  for  devices  that  exhibit  a 2DEG.  It  was  shown  l„ 
Chapter  V that  the  mobility  Is  a critical  parameter  In  the  determination 
of  the  diffusion  coefficient.  Therefore,  It  is  suggested  that  scatter- 
ing mechanisms  that  are  not  included  In  the  simulation  program  will  be 
added  to  It.  Especially,  alloy-scattering  and  Interface-scattering 
proceases  might  Influence  the  drift  velocity-field  curve  strongly. 
Also,  the  Inclusion  of  real-space-charge-transfer  effects  In  the 
simulation  will  strongly  modify  the  mobility  profile,  because  carriers 
move  Into  the  larger  bandgap  material  as  soon  as  their  kinetic  energy  is 
of  the  order  of  the  conduction-band  discontinuity. 
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The  Monte  Carlo  simulation  technique  is  a powerful  tool  to  inves- 
tigate charge  transport  in  all  kinds  of  devices.  In  this  dissertation 
simulations  are  presented  for  charge  transport  that  takes  place  parallel 
to  the  heterointerfaces.  Most  of  the  published  literature  in  this  area 
so  far  is  concerned  wich  parallel  transport.  Recent  technological 
improvements  in  MBE  fabrication  made  it  possible  to  produce  very  high- 
quality  interfaces  which  allow  for  device  applications  in  which  the 
electrons  propagate  perpendicular  to  the  interface.  Well-known  examples 
are  superlattice  avalanche  photodetectors , effective-mass  filters,  and 
superlattices.  Also  in  this  fast-developing  area,  simulations  of  this 
electron  travel  will  give  a deeper  understanding  of  the  processes 
involved  in  the  transport.  Since  carrier-carrier  interactions  are  very 
important  here,  the  simulations  should  involve  some  ensemble  averaging. 
More  than  one  electron  needs  to  be  followed  in  Its  travel  from  cathode 
to  anode. 

6.4.  Diffusion  Coefficient  of  AlOaAs/OaAs/AlGaAs 

Measurements  of  the  noise  temperature  in  the  frequency  range  of 
0.5  - 12  GHz,  together  with  impedance  measurements,  resulted  in  a dif- 
fusion electric-field  curve.  The  noise  temperature,  measured  for 
electric-field  strengths  up  to  3 kV/cm,  is  white  up  to  2.2  GHz,  and  then 
a roll-off  occurs  for  all  field  strengths.  Efforts  to  measure  the  impe- 
dance of  the  device  resulted  in  the  establishment  of  the  values  of  the 
parasitic  elements  that  are  present  in  the  termination  of  the  stripline. 
The  low-frequency  differential  conductance  was  obtained  from  the  pulse- 
biased  current-voltage  characteristics. 

The  resulting  diffusion-coefficient  electric-field  curve  differs 
from  measurements  on,  and  simulations  of,  bulk  GaAs.  The  observed  onset 
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of  Che  Increase  In  the  data  occurs  at  a smaller  eleccrlc-f leld  strength 
than  is  the  case  ior  bulk  GaA.s.  Simulations  of  the  charge  transport  in 
a 2DEG  GaAs  well  show  the  peak  to  occur  at  approximately  the  same  field 
strength. 

The  same  comparison  was  made  for  the  field  dependence  of  the  noise 
temperature.  Again  the  onset  of  the  increase  of  the  experimental  2DEG 
noise  temperature  occurs  at  lower  fields  when  compared  with  the  onset  of 
bulk  measurements  and  simulations.  However,  the  agreement  with  the 
simulated  noise  temperature  is  much  better. 

The  main  reason  for  this  discrepancy  is  the  high  mobility  that  was 
found  from  the  simulations.  This  causes  hot-electron  effects  to  occur 
at  lower  fields.  Therefore,  it  is  suggested  that  accurate  mobility 
measurements  be  performed  and  the  Monte  Carlo  program  be  adjusted  in 
such  a way  that,  by  including  necessary  scattering  and  transfer  mechan- 
isms, the  velocity  field  curves  agree.  It  is  expected  that  the  new 
program  will  yield  results  that  match  the  measured  diffusion-coefficient 


data  much  better. 


APPENDIX  A 
SUBROUTINE  POLOPT 


Subroutine  POLOPT 

The  subroutine  is  initiated  in  the  program  by  the  statement  CALL 

POLOPT  (RKXI,  RKYI,  RKZI,  REI,  REFIN),  where  RKXI  is  the  x component  of 

the  wave  vector  just  before  scattering  on  input.  On  output  RKXI  is  the 

x component  of  the  wave  vector  just  after  scattering;  likewise  for  RKYI 

and  RKZI.  REI  and  REFIN  are  the  energy  of  the  electron  just  before  and 

just  after  the  polar  optical  phonon  scattering,  respectively.  Vector  K 

(RKXI,  RKYI,  RKZI)  is  renamed  and  normalized:  A = if  • K .Also,  the 

energies  are  normalized,  but  they  keep  the  same  name: 

T 2 

REI  + REI  • (in 

T 2 

REFIN  > REFIN  • (— ) 

The  vector  9 ulll  describe  the  wave  vector  state  after  scattering.  See 
Fig.  A.l  for  the  definition  of  A and  3.  At  the  end  of  the  routine  3 

win  be  deno realized.  The  lengths  of  the  vectors  B and  A are  BdORM  and 
ANORM  respectively.  The  angle  3 generates  a cone  with  half  angle 
3 and  central  awls  A.  The  base  of  the  cone  la  a circle  with  center 

|5|  cos ( 3 ) -||-and  radius  W (w  - |||  sln(B)).  This  circle  Is  parame- 
trized  as: 


-T 

|3|  cos(3)  + lx  + 

1*1 


4 

ny 


(Al) 


176 


177 


Fig.  A. 1 . Definition  of  wave  vectors  before  and  after 
scattering:  (a)  actual  wave  vectors; 

(b)  the  normalized  wave  vectors. 
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where  x and  y are  vectors  In  the  plane  perpendicular  to  A.  Hence 
A y - 0 and  A»y  = 0.  We  will  also  require  x»y  = 0 
Choose: 


x = (1,  0,  - A1/A3); 


then 


|x|  = /l  + Al2/A32  = El 


(A2) 


and 


y = f Ai  _ Al2  + A32  1 . 

^ A3  ’ A2*A3  ’ LJ  ’ 


then 


01  + 

A3  (A2*A3)2 


E2 


(A3) 


El  and  E2  are  defined  in  the  routine  as  stated  above. 

Also,  A = and  u = W cos(^) 

El  ana  u £2  ‘ The  larSest  vector  that  can  be 

allowed  according  to  the  final  energy  REFIN  is: 


h2(K' 1 2 h2 (k2+k2 ) 2 „ 

E = — >, . _ z y , h ,2 

FIN  _ * 1T~  + — 


2m 


2m 


2m 


2 2 2 


> ^ ^ 2 _ h NMAX  *tt 


FIN  ■>  * X- 

2m*  max  2m 

L2  • E 2 

FIN  . h 2 

2 > ~ NMAX  > 

tt  2m 


(A4) 


2 2 

where  L is  the  normalized  final  energy 


2m  2 

— r • REFIN  > NMAX  . 

h 
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Remember  that  j 3 j = BNORM  = /2*m""*REFIN/h2  , and  we  find  NMAX  = INT 

(3 NORM) . The  largest  allowed  k;<  component  Is  NMAX-ir/L,  and  the  largest 

Bx  component  is  NMAX  * 1.0.  The  final  energy  REFIN  generates  a shell  of 

radius  /2‘m  *REFIN/h  in  It  space.  For  quantized  k components  this 

\ 

results  in  circles  of  constant  energy  in  each  plane  for  n = 1,  2,  ..., 
NMAX.  To  find  a final  wave  vector  state,  after  scattering,  the  inter- 
sections of  these  circles  wich  the  base  of  the  cone  have  to  be  found. 
The  constant  energy  circles  are  described  by 

(I,  0,  0)  +/jB|2  - I2  {(0,1,0)  cos  4>  + (0,0,1)  sin  (A5) 

where  1=1,  2,  ...,  NMAX  and  (0,1,0)  and  (0,0,1)  are  unit  vectors  along 
the  y and  z axes,  respectively.  Equating  the  first  components  of 
eq.  (Al)  and  eq.  (\5)  yields 


j B | cos(B) 


777 + w sin^>  It  + 

{ A | 


W cos(j») 
E2 


which  is  rewritten  in  the  form 


(A6) 


Pl*sin($)  + P2*cos('^)  = P3 


where 


(A7) 


PI  a W/El  (AS) 

P2  = W»A1/(E2*A3) 


(A9) 


and 


P3  = I - j B j cos(S) 


(A10) 


We  further  define 
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P4  = /?12  + P22  , 

and  it  can  be  seen  from  eq.  (A7)  that  there  are 
1 < = TEST  < 1 . Rewrite  eq.  (A7) 

cos(v)sin(iJ>)  + sin( v)cos(ii)  = P3/P4 


or 


sin(ij)  + v)  = P3/P4  . 


Then 

i|>  = arcsin(P3/?4)  - v 
or 

^ = it  - arcsin(P3/?4)  - v . 

These  two  solutions  for  'j)  are  stored  in  PSI(2I)  and 
tively;  for  each  plane  I we  find  two  solutions, 
with  eq.  (A12),  one  readily  sees 

cos(v)  = P1/P4 


and 


sin(v)  = P2/P4  . 
Solutions  for  eq.  (A14)  are 

v = ± arccos(Pl/P4)  . 


(All) 

only  solutions  when 
(A12) 

( A1 3) 


PSI(2I  + 1) , respec- 
Corapariag  eq.  (A7) 

(A14) 
(A15) 
( A16) 


Solutions  for  eq.  (A15)  are 
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v = arcsin(P2/P4) 


(A17a) 


or 


v = it  - arcsin(P2/P4)  . 


( A1 7 b ) 


All  of  the  solutions  (A16)  are  called  RNU1  and  RNU2,  respectively,  In 
the  routine.  RNU3  is  defined  as  (A17a)  and  RNU4  as  (A17b).  The  solu- 
tion sought  for  v Is  the  one  that  gives  the  identical  result  in  (A16) 
when  compared  with  (A17). 

All  of  the  lively  final  K'  states  are  characterized  by  the 
angle  ifi  for  each  plane.  We  assume  that  all  of  these  angles  * ace 

equally  likely  in  finding  the  final  state.  With  the  help  of  a random 
number  the  angle 

^0  ( = PSIO)  is  chosen,  and  all  components  of  B can  be  calculated.  The 

vector  B is  multiplied  by  £ to  find  K’  . .Also,  the  energies  RE  I and 
REFIN  are  denormalized  before  returning  to  the  main  program. 


Variables 

For 

simulation 

Variable 

AA 

BB 

C 

C1,C2, . . .C9 
CHECK 


in  Main  Program 

completeness  all  of  the  variables  used  in  the  Monte  Carlo 
program  are  listed  below: 

Physical  Quantity 

local  variable,  used  for  the  calculation  of  the 

emission  and  absorption  rates 

local  variable,  see  AA 

unit  of  charge  in  e.s.e.  units 

constants  in  the  scattering  rates,  as  defined  in  [39] 
~ 3 2—  — 3 

v - 3v  v + 2(v)  , this  quantity  should  be  zero  for  a 


Gaussian  distribution  function 
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COSGAM 

cos(y),  where  y is  the  angle  between  the  possible 
final  wave  vector  and  the  quantization  direction 

D 

energy  difference  between  the  bottoms  of  the  central 
and  satellite  valleys 

DIF 

diffusion  coefficient,  calculated  from  the  second 
moment  of  the  distribution  function 

DIF4 

diffusion  coefficient,  calculated  from  the  fourth 
moment  of  the  distribution  function 

DK 

mesh  size  in  wave  vector  space 

E 

charge  of  one  electron  in  SI  units 

EF(i) 

energy  of  the  electron  after  scattering  if  channel 
^ 1>2, . . . ,10)  were  chosen 

EFIN 

energy  of  the  electron  after  scattering  when  the 
scattering  mechanism  has  been  chosen 

El 

energy  of  the  electron  at  the  end  of  the  ballistic 
flight  just  before  scattering 

EM 

effective  mass  of  the  electron 

EMI , EM 2 

effective  mass  of  the  electron  in  the  central  and 
satellite  valleys,  respectively 

EMAX 

the  maximum  energy  the  electron  is  allowed  to  have; 
it  is  mainly  used  to  calculate  T (=  GAMMA) 

EMAX1 

the  maximum  energy  an  electron  reaches  in  the  course 
of  an  actual  simulation 

EMIN1 

the  minimum  energy  an  electron  has  during  an  actual 
simulation 

ETOT 

the  kinetic  energy  an  electron  gains  during  one 
sampling  interval  FT 
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EZERO(i) 

the  quantum  mechanically  lowest  energy  an  electron 
can  assume  in  the  central  (i=l)  or  satellite  (i=2) 
valley 

F(i) 

electric  field  strength 

FT 

length  of  one  sampling  time  interval 

GAMMA ( i ) 

sum  of  all  of  the  scattering  rates,  including  self- 
scattering, for  the  central  (i=l)  or  satellite  (t=2) 
valley 

CK(i) 

array  used  to  check  the  proper  operation  of  the  sub- 
routine POLOPT 

GMAX 

counter  to  keep  track  of  how  many  times  the  electron 
energy  exceeds  the  set  value  EMAX 

H 

Planck's  constant  divided  by  2tt 

HWE 

energy  of  an  equivalent  intervalley  phonon  in  eV 

HWI 

energy  of  a nonequivalent  intervalley  phonon  in  eV 

HWO 

energy  of  a polar  optical  phonon  in  eV 

HWSTAR 

effective  phonon  energy  used  in  the  calculation  of 
the  polar  optical  scattering  rates,  hoi*  in  ref.  [37] 

I 

counter,  index  of  the  scattering  rate  array  L 

IDIF 

counter,  number  of  time  samples 

II 

counter 

J 

counter,  number  of  electric  field  strengths 

JM 

counter,  used  in  the  calculation  of  the  scattering 
rates  to  sum  over  the  allowed  wave  vector  planes 

K 

counter 

KB 

Boltzmann's  constant 

KRHO 

(KXI*KXI+KYI*KYI)1/2 
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kxf,kyf,kzf 

x,  y,  and  z components  of  the  wave  vector  after 

KXI,KYI,KZI 

scattering,  at  the  beginning  of  next  ballistic  flight 
x,  y,  and  z components  of  the  wave  vector  before 

KZKY2 

scattering,  at  the  end  of  the  ballistic  flight 
(KZ*KZ+KY*KY)li/2 

KZTOT 

the  total  amount  the  z component  of  the  wave  vector 
changed  during  one  sampling  time 

L(i) 

scattering  rate  of  process  i (i=l,2, . . . ,10) 

LENGTH 

width  of  the  quantum  well 

LL 

counter 

M 

mass  of  a free  electron 

MAX(i) 

integer  that  determines  the  largest  x component  of 
the  wave  vector  under  the  condition  that  the  electron 
energy  is  smaller  or  equal  to  EF(i) 

N1 

integer  identifying  the  kx  plane  the  electron  is  in 

NAC 

number  of  deformation  potential  scatterings 

NDIF 

total  number  of  sampling  times 

NE 

equivalent  intervalley  phonon  occupancy 

NF 

total  number  of  electric  fields 

NI 

nonequivalent  intervalley  phonon  occupancy 

NMAX 

see  MAX(i) 

NO 

polar  optical  phonon  occupancy 

NOPT 

number  of  polar  optical  phonon  scatterings 

PHI 

angle  used  to  determine  the  wave-vector  state  after 
scattering 

PI 

7T 

PP 

local  variable 
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QQ 

local  variable 

R 

random  number 

R1 

high-frequency  relative  dielectric  constant 

R2 

low-frequency  relative  dielectric  constant 

RHE 

mass  density  of  the  material  (g/cm3) 

S 

velocity  of  sound  in  the  material 

SCATT  (I) 

scattering  rate  of  process  i 

SINGAM 

sin  (y),  see  COSGAM 

SR 

number  of  real  scatterings 

SS 

number  of  self-scatterings 

SUM 

local  variable,  used  to  calculate  the  scattering 
rates 

T 

temperature 

TF 

sampling  time,  time  interval  length  of  one  average 

THA 

deformation  potential  for  acoustic  phonon  scattering 
(eV) 

THE 

deformation  potential  for  equivalent  intervalley 
phonon  scattering  (eV/cm) 

THI 

deformation  potential  for  nonequivalent  intervalley 
phonon  scattering  (eV/cm) 

TLME 

ballistic  flight  time 

TLEV  (i,j) 

time  spent  in  the  kx  plane  of  the  i1-*1  valley 

T0CC1 

time  spent  in  the  central  valley,  i.e.,  summed  over 
the  planes 

TOCC2 

time  spent  in  the  satellite  valley,  i.e.,  summed  over 
the  k^  planes 
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TT 

flag;  if  TT-0  the  simulation  continues  and  if  TT=1 
the  initial  total  scattering  rates  GAMMA  (1)  and 
GAMMA  (2)  are  calculated. 

TTQT 

sum  of  the  ballistic  flight  times 

V 

indicates  in  which  valley  the  electron  is;  V=1 

V2AVG 

central  valley,  V=2  satellite  valley 

~2 

V 

V3AVG 

~3 

V 

V4AVG 

T 

V 

VAVG 

V 

VI 

ETOT/KZTOT 

VR 

distance  from  the  k^  axis 

VV 

see  V 

WE 

equivalent  intervalley  phonon  frequency 

WI 

nonequivalenc  intervalley  phonon  frequency 

WO 

polar  optical  phonon  fraciuency 

Variables 

in  Subroutine  POLOPT 

A1,A2,A3 

normalized  x,  y,  and  z component  of  the  wave  vector 
just  before  the  polar  optical  scattering  occurs 

ANORM 

length  of  the  vector  A = (A1,A2,A3) 

B1,B2,B3 

normalized  x,  y,  and  z components  of  the  wave  vector 
just  after  the  polar  optical  scattering  event 

BETA 

cos (3),  where  3 is  the  angle  between  the  wave  vectors 
before  and  after  the  scattering 

BN 

length  of  the  vector  B = (B1,B2,B3) 

BNORM 

see  BN 

El,  E2 

vector  norms  defined  in  this  appendix,  p.  177 
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EX 

defined  in  ref.  [39] 

F 

local  variable 

I 

counter 

11 

counter 

II 

counter 

J 

counter 

K 

counter 

KO 

counter 

N 

counter 

NMAX 

integer  that  determines  the  largest  x-component  of 
the  wave  vector  for  which  the  electron  energy  is  not 
larger  than  REFIN 

NPLANE 

number  of  the  k^  plane  that  has  been  chosen 

PI , P2 ,P3 

variables  in  the  equation  Pl*sin(^)+?2*cos(^)  = ?3 
used  to  solve  for  i(>. 

P4 

(P1*?1+P2*P2)1/2 

PSI  (i) 

solution  i for 

PS  10 

ip  chosen  after  scattering 

R 

random  number 

RE 

electron  charge 

REFIN 

energy  of  the  electron  after  scattering;  It  is  the 
same  as  EFIN  In  the  main  program 

REI 

energy  of  the  electron  before  scattering;  it  is  the 
same  as  El  in  the  main  program 

REM 

mass  of  the  electron;  it  is  the  same  as  EM  in  the 
main  program 

RH 

Planck's  constant 
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RKXI,RKYI,RKZI 

x,  y,  and  z components  of  the  wave  vector 

RLENGTH 

width  of  the  quantum  well 

RM 

free  electron  mass 

RNU 

chosen  angle  v 

RNU1 ,RNU2, 
RNU3.RNU4 

the  four  possible  solutions  for  v 

TEST 

criterion  to  decide  whether  solutions  for  'Jj  exist 

IJ 

random  number 

w 

radius  of  the  circle  generated  by  R,  K' , and  3. 
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c computer  simulation  of  the  drift  velocity'mean" 
c energy  and  distribution  function  of  hot 

c ELECTRON 


* 

* 


c= 

c 

c= 


C== 

c 

c== 

41 

42 

43 

44 
1 


3 

4 


j BB 

DOUBLE  PRECISION«6  F( 10) , GAMMA (2) ,EF(10)  SCATTfim 

™ecisio».6  Ei.mi.acsjfkS! 

1 ,DK(2) ,L( 10 ) ,EZER0(2) , TLEV( 1 :2, 0 : 20 ) 

NF’  VR» TT,  V,  SR , SS , K,  GMAX,  NDIF 

1 , NACC , NOPT , MAX ( 10) 

COMMON/POLOP/M, EM, E, H, LENGTH 

FUNDAMENTAL  CONSTANTS  =================== 

H=1 .05459  ============================= 

E=1 .6021 9 
C=4. 803213 
KB=1. 38062 
M=9. 109136 
PI=4.0*ATAN( 1.0) 

DATA  ON  MATERIAL  GaAs  ==========================: 

WRITE (6,41 ) ====================================: 

FORMAT ( 10X, 'MATERAL  DENSITY(GM. CM-3) » ) 

READ(8,689)RHE 

WRITE(6,42)RHE 

S™M89)s,0',’/',ra’'VELOCm  0F  s°™><io».5cm,s-i)') 

WRITE (6, 43 )S 

mDl(M89)Bj0',,’/’'OX’,HIGB  FKE0U£:,iCI  oielectbic  CONSTANT') 
WRITE(6,44)R1 

m“8!689)B2°',|’/’10X’’LO"  FBEWEIiCI  ^LECTBIC  CONSTANT') 

WRITE ( 6 , 1 ) R2 

31SS' 'V.10X, -OPTICAL  PHONON  FBEQUENCY<  10»13BAD.S-1 ) ') 
WRITE (6 ,2) WO 

ENTEBVALLET  PHONON  PBEO.  ' , 

WRITE ( 6 ,3 ) WE 

EEA“M89)W0'',’/’'0X’'ImmSIXEr  PH°"°“  FREQUENCY* ) 

WRITE ( 6 , 4 ) WI 

REA^8i689)TH^4’/,10X,,AC°USTIC  DEFORMAT*ON  POTENTAL(EV)  ' ) 
WRITE(6,5)THA 

^ 6 1 9) Jhe * ^ ' 1 °X ' ' EQ U* INTER VALL^Y  DEF. POTENTIAL(EV) • ) 


o o o o 
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6 

7 

8 

9 

689 

C = = : 
C 

C=== 

46 

10 

21 

11 

18 

533 

12 

534 

13 

14 
30 

23 

536 
22 

537 


WRITE (6, 6) THE 

READ( 8 (6 8 9) THI ’ ^ ’ ' ’ 1 °X ' ' INTERVALLEY  DEF- POTENTIAL ' ) 
WRITE(6,7)THI 

nuSI(JoX^E1°*4’/,10X’  'VALLEY  SEPARATION') 

READ(8,689)D 

WRITE (6, 8 )D 

READ^S^SgjEMl’^^^^'  'CEN™L  VALLEY  EFFECTIVE  MASS') 
WRITE(6 , 9) EMI 

vallei  eff-mass,) 

FORMAT (El  0.4) 


FINAL  DATA  INPUT 


WRITE (6 , 10 )EM2 

or?!!/o^1nX!E10'il  »10X>  'TEMPERATURE (KELVIN)  ') 

READ(8,689)  T 

WRITE (6, 21)  T 

FORMAT ( 15X,E10. 4 ,/ ,10X, 'MAXIMUM  ENERGY(EV) ' ) 

READ (8, 6 89)  EMAX 
WRITE (6,11)  EMAX 

FORMAT ( 1 5X , El  0 . 4 , / ,10X, 'NUMBER  OFAVERAGES') 

READ (8, 533)  NDIF 
FORMAT (15) 

WRITE(6, 12)  NDIF 

READ(8(534jIN^,/’10X’,NUIffiER  °F  ELECTRIC  FIELDS') 

FORMAT (12) 

DO  30  1=1 ,NF 
WRITE(6,13)  I 
FORMAT (1 OX, 'FIELD  # »,I2) 

READ(8,689)  F(I) 

WRITE(6, 1 4)  F(I) 

FORMAT(15X, 'FIELD  = 'E10.4,'  kV/cm') 

CONTINUE  ‘ 

WRITE (6, 23) 

READ(8!536)'wLLEY  F°R  DISTRIBUTI0N  UNCTION  :1  OR  2 • ) 
FORMAT (12) 

WRITE (6 ,22) 

REA^8!53^),?fTANCE  FR°M  KZ  AXIS  °F  DISYRIBUTION  :1  TO  21') 
FORMAT (12) 


SL™elLEi"e(2) vaTle?!.  GBOl,NI,LEVELS  “ 


LENGTH=2.0E-6 
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EZER0(1)=1.0E-14*H«H*PI*PI/(LENGTH*LENGTH*2.0«E*M»EM1) 

EZERO ( 2 ) =EZERO ( 1 ) «EM 1 / EM2 


C CALCULATE  PHONON  FREQUENCIES  AND  OCCUPATION  RATIOS 


HWO=H»WO/(E* 100.0) 

HWI=H»WI/(E*100.0) 

HWE=H*WE/(E» 100.0) 

IF(WO.  NE.0.0)  GO  TO  90  9 
N0=0 . 0 
GO  TO  908 

C=======  H*1E-27*W0*1E13/(E*1E-16)=W0»76.385 

909  NO=1/(EXP((WO»76.385)/T)-1) 

908  NI=1/(EXP((WI«76.385)/T)-1) 

NE=1/(EXP((WE*76.385)/T)-1) 

C CONSTANTS  FOR  SCATTERING  RATES 


,(1/R1-1/R2)/(4.0 *L  ENGTH  *E ) 


C1=1.0E+5*C*C*WO*(NO+1)i 
C2=C1*N0/( NO+1 ) 

C3  = 1.0E+3*THA»THA»E»E*KB*T»M/(i1.0*H»H*H*LENGTH»RHE»S«S 

C5-C4*^/*NEE*)HE*E*E*M*(f,E+1  ^/(i,'0*RHE*WE*H,H*LENGTH) 

C6  = 1.0e+7*thi*thi*e*e*m#(ni+i)/(24#°#rHe#Wi«h*h«lengTH) 

rrp72MT!^I*!HI#E#E*M*(NI+1)/(i<'0*RHE‘WI*H‘H!,LENGTH) 
oo=Of*NI/( NI+1 ) 

C9=C6*NI/( NI+1 ) 


t 


£ CALCULATE  K SPACE  MESH  ELEMENT  FOR  BOTH  VALLEYS 


DK ( 1 ) = 1 . 0 E+7  *SQ  RT ( 2 *EM 1 »M*EMAX*E ) / ( H*  20 . 0 ) 

DK ( 2) = 1 . OE+7 »SQRT ( 2»EM 2 *M*EMAX*E ) / ( H*20 . 0 ) 


C 

C 

C 

C: 


CALCULATE  VALUES  OF  KZ  AT  CENTRE  OF 
IN  CHOSEN  VALLEY  AND  WRITE  TO  FILE 


EACH  MESH  ELEMENT 


* 


KZI=-DK(W)«9.5 
DO  805  LL=1 ,20 
WRITE (5, 807)  KZI 
807  F0RMAT(E12.4) 

KZI=KZI+DK(W) 
805  CONTINUE 
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F0R  CENTRAL  valle^then  calculate  the 

OF  ENERCT^  rm  F°R  REAL  PR0CESSES(R)  FOR  A NUMBER 

TN  ™ T0  THE  SI2E- store  maximum  value  of  r 

IN  GAMMA ( TO  CALCULATE  PSEUDO  (SELF)  SCATTERING  RATE 

TT=1  ENABLES  PROGRAM  TO  RETURN  TO  LABEL  40 


31 

TT=1 

EM=EM1 

V=1 

GAMMA(V)=0.0 

35 

EI=EZER0( V) 

J=1 

EI=EI+(EMAX-EZERO( V) ) /20 . 0 

40 

GO  TO  100 
R=0.0 

50 

DO  50  1=1,10 

R=R+L(I) 

CONTINUE 

IF(R. GT. GAMMA(V)) 

GAMMA( V) =R 

J=J+1 

IF(J. NE.21 ) GO  TO 

35 

C=  = = : 

C 

C 

C==  = = 


S'S™  S4TEUJTE  A»Dlim7;R0CEsrra! 


71 
75 

72 


IF( V. EQ. 2)  GO  TO  71 

EM=EM2 

V=2 

GO  TO  31 

WRITE (5, 75)  GAMMA ( 1 ) , GAMMA (2) 
WRITE (7 ,75)  GAMMA (1), GAMMA (2) 
FORMAT ( 1 0X , ' GAMMA ( 1 ) = 1 f El  0 . 4 ' 
WRITE(5,72)  T, NDIF, EMAX 
WRITE (7 ,72)  T, NDIF, EMAX 

1 ) 1 °X  ’ ' TEMP=  ' ’ E1  ° ^ ’ 5X  ’ ' n° 


GAMMA ( 2) = ' ,E10 . 4 ) 


of  flights  = ' , 14 , 5X , ' MAXENERGY= • 


C = = = = : 

C 

C 

C 

C 

C===== 


After  outputing  Gamma(l)  and  Gamma(2)  both"are"divided=bv=in= 

30?tterin6  - in  tie6 

r.miifii  k negative  scattering  rates  and  then  increases 
O^aUJ  cr^Oa-MCs)  by  10  %,  depending  on  the  Valle) 


GAMMA ( 1 ) =GAMMA ( 1 ) / 1 0 . 0 

GAMMA ( 2) =GAMMA ( 2) / 1 0 . 0 
WRITE(5,75)  GAMMA (1), GAMMA (2) 
WRITE(7,75)  GAMMA( 1 ) , GAMMA(2) 


C SET  MESH  REGISTERS  TO  ZERO  AND  PLACE  ELECTRON  AT 


o o 
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STARTING  POINT  IN  MESH. TT=0  FOR  ITERATIVE  PROCESS 


TT=0 

J=1 

TTOT=0.0 
VAVG=0.0 
V2AVG=0 . 0 
DIF=0.0 
V3AVG=0.0 
V4AVG=0.0 
DIF4=0.0 
CHECK =0.0 
TF=4.0E-10 
FT=TF 
IDIFrO 
ETOT=0.0 
KZT0T=0 . 0 
NPRINT=0 
80  V=1 

EM=EM1 
N1  =1 

KXF=PI/ LENGTH 
KYF=1.0E6 
KZF=1 . 0E6 

KRHO=SQRT(KXF*KXF+XYF*KYF) 
KT=SQRT (KRHO*KRHO+KZF*KZF ) 

SR=0 
SS=0 
NOPT=0 
NACC=0 
NPOLOPT=0 
EMAX 1=0.0 
EMIN1=1000.0 
MNE=0.0 
DO  85  11=0,20 
TLEV( 1 ,11) =0.0 
TLEV(2,II)=0.0 
85  CONTINUE 
GMAX=0 

EFIN=H*H*KT*KT* 1 . 0E-1 4/ (E»2»EM»M) 


C IF  NO.  OF  REAL  PROCESSES  EQUALS  CHOSEN  VALUES*”" 

END  ITERATION  AND  GO  TO  FINAL  CALCULATION. 

c =================== 

CALL  RANDOM  NUMBER ( N0T=0)  AND  CALCULATE  TIME  OF 
FLIGHT  UNDER  ELECTRIC  FIELD  AND  NEW  POSITION  OF 
ELECTRON  K SPACE 


90  IF(TTOT.GT.FT)  GO  TO  470 
R=RANU(0.0 ,1.0) 


o o o o o 
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IF(R.EQ.O.O)  R=1 . 0E-20 
TIME=-ALOG( R ) / GAMMA ( V) 

KXI=KXF 

KYI=KYF 

KZI=KZF+(TIME*E*F (J)*1.0E+18)/H 
KRHO=SQRT(KXI »KXI+KYI «K  YI ) 
KT=SQRT (KRHO*KRHO+KZI*KZI ) 
EI=H*H*KT*KT*1 .OE-1 4/( E*2*EM*M) 
Em ax 1 =Dmax 1 ( Emax 1 , El ) 

Eminl =Dmin1 (Eminl ,EI) 

TLEV( V, N1 ) =TLEV( V, N1 )+TIME 


IF  ELECTRON  LEAVES  MESH  PLACE  IT  ON  EDGE  OF  MESH  AND 
REGISTER  OCCURRENCE  IN  COUNTER  GMAX 


IF(EI.LE.EMAX)  GO  TO  95 

GMAX=GMAX+1 

EI=EMAX 

KT= 1 . OE+7  *SQRT ( 2*EM*M*EMAX*E ) / H 
IF(KZI. GT.0.0)  GO  TO  94 
KZI=-SQRT(ABS(KT*KT-KRHO*KRHO) ) 
GO  TO  95 

94  KZI=SQRT(ABS(KT*KT-KRHO*KRHO)) 


STORE  FLIGHT  TIME  IN  TOTAL  TIME  REGISTER  FOR  APPROPRIATE 
VALLEY, THEN  REGISTER  PASSAGE  OF  ELECTRON  THROUGH  ELEMENTS 
OF  K SPACE  MESH 


95  TTOT=TTOT+TIME 

KZTOT=KZTOT+ABS ( KZI-KZF ) 

ETOT=ETOT+(KZI*KZI-KZF*KZF) * 10 . 0*H/ ( 2*EM*M) 
MNE=MNE+(KZI«KZI*KZI-KZF*KZF»KZF)«1.0E-14»H*H/(2.0*EM*M»E) 


SUM  TOTAL  CHANGES  IN  K SPACE  POSITION  AND  ENERGY  SPACE  POSITON  * 
AND  STORE  IN  KZTOT  AND  ETOT.SUM  MEAN  ENERGY  CHANGE  IN  MNE  • 
SUM  INDIVIDUAL  VALUES  FOR  EACH  VALLEY  IN  KZ(V)  AND  ET(V)  * * 


CHECK  FOR  ROUNDING  ERRORS  LEADING  TO  NEGATIVE  ENERGY  VALUES  * 
IF  THIS  OCCURS, PLACE  ELECTRON  AT  STARTING  POSITION.  * « 


IF(EI. GT. EZERO( V) ) GO  TO  100 
N1  =1 

KXF=PI/ LENGTH 
KYF=1.0E6 
KZF=1 . OE+6 

KRHO=SQRT(KXF*KXF+KYF*KYF) 

EM=EM1 
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V=1 

GO  TO  90 


c=======--------I!.!^AL  ENERGY  VALUE  F0R  eachIcattering'process 

100  EF(D=EI-HWO 
EF(2) =EI+HWO 
EF(3)=EI 
EF(i|)=EI 
EF(5) =EI-HWE 
EF(6) =EI+HWE 
EF(7) =EI-HWI+D 
EF(8)=EI-HWI-D 
EF(  9)  =EI+HWI-D 
EF(10)=EI+HWI+D 


C _ SCATTERING  RATES  FOR  REAL  PROCESSES 


EMISSION  OF  A POLAR  OPTICAL  PHONON 


110 


111 

c 

c 

c 

120 


125 


IF(EF( 1 ) . GT. EZERO( V) ) GO  TO  110 
L( 1 )=0.0 
GO  TO  120 

MAX( 1 ) =SQRT(EF( 1 ) /EZERO( V) ) 
SUM=0.0 

DO  111  JM=1 ,MAX( 1 ) 

HWSTAR=HWO-( N1 *N1 -JM*JM)*EZERO( V) 
IF(JM.EQ.NI)  THEN 
AA=2.0 
ELSE 


AA=1 . 0 


END  IF 

SUM=SUM+AA 
CONTINUE 
L(  1 )=C1*SUM 


ABSORPTION  OF  A POLAR  OPTICAL  PHONON 

IF(EF(2).GT.EZER0(V))  GO  TO  125 

L(2)=0.0 

GO  TO  130 

MAX(2) =SQRT(EF(2) /EZERO( V) ) 

SUM=0.0 

DO  121  JM=1,MAX(2) 

HWSTAR=HW0+(N1 *N1-JM*JM)*EZER0(V) 
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121 

C 

C 

C 

130 


135 


131 

C 

C 

C 

no 


145 

C 

C 

c 

150 


155 


160 


IF(JM.EQ.NI)  THEN 
AA=2. 0 
ELSE 


AAsl.o 
END  IF 

AA=AA/SQRT( ( JM-N1 ) **4*EZER0( V) *EZERO( V)+2.0*( JM-N1 )**2#EZER0f V) 
@*(2.0*EI+HWSTAR)+HWSTAR*HWSTAR)  ; EZERO(V) 

fl^;^!’?^QRT((JM+N1)”4*EZER0(V)*EZER0(V)+2.0»(JM+N1)««2 
§*EZERO ( V) « ( 2. 0 *EI+HWSTAR ) +HWSTAR*HWSTAR ) 

SUM=SUM+AA 

CONTINUE 

L(2)=C2*SUM 


EMISSION  OF  AN  ACCOUSTIC  PHONON 

IF(EF(3) . GT. EZERO( V) ) GO  TO  135 

L(3)=0.0 

GO  TO  140 

MAX(3)=SQRT(EF(3)/EZERO(V)) 

SUM=0 . 0 

DO  131  JM=1 ,MAX(3) 

IF(JM.EQ.NI)  THEN 
AA=3.0 
ELSE 

AA=2. 0 
END  IF 
SUM=SUM+AA 
CONTINUE 
L(3)=C3»EM*SUM 


ABSORPTION  OF  AN  ACOUSTIC  PHONON 

IF(EF(4) . GT.  EZERO( V) ) GO  TO  145 
L(4)=0.0 
GO  TO  150 
MAX(4) =MAX(3) 

L(4)=L(3) 


EMISSION  OF  AN  EQUIVALENT  INTERVALLEY  PHONON 

IF(EF(5) . GT. EZER0(2) ) GO  TO  155 

L(5)=0.0 

GO  TO  170 

IF(V. EQ.2)  GO  TO  1 6 0 

L(5)=0.0 

GO  TO  170 

MAX (5) =SQRT(EF(5)/EZER0(2)) 

SUM=0.0 

DO  151  JM=1,MAX(5) 

IF(JM. EQ.  N1)  THEN 
AA=3 . 0 
ELSE 

AA=2.0 
END  IF 
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SUM=SUM+AA 
151  CONTINUE 

L(5)=C4*EM2*SUM 

C 

c absorption  of  an  equivalent  intervalley  phonon 

170  IF(EF(6) . GT. EZER0(2) ) GO  TO  175 
L(6)=0.0 
GO  TO  1 90 

175  IF(  V.  EQ.  2)  GO  TO  1 80 
L(6)=0.0 
GO  TO  1 50 

180  MAX ( 6 ) =SQRT ( EF ( 6 ) / EZERO ( 2 ) ) 

SUM=0.0 

DO  171  JM=1 ,MAX(6) 

IF(JM.EQ.NI)  then 
AA=3 .0 
ELSE 
AA=2.0 
END  IF 
SUM=SUM+AA 
171  CONTINUE 

L(6)=C5*EM2*SUM 

C 

c EhIISSI0N  0F  NON-EQUIVALENT  INTERVALLEY  PHONON  (SAT.  TO  CENTRAL) 

190  IF(EF(7).GT.EZERO(1))  GO  TO  1 95 
L(7) =0.0 
GO  TO  210 

195  IF(  V.  EQ.  2)  GO  TO  200 
L(7)=0.0 
GO  TO  210 

200  MAX (7) =SQRT(EF(7) /EZERO ( 1 ) ) 

SUM=0.0 

DO  191  JM=1  ,MAX(7) 

IF(JM.EQ.NI)  THEN 
AA=3 .0 
ELSE 

AA=2.0 
END  IF 
SUM=SUM+AA 
1 91  CONTINUE 

L( 7) =C6*EM1 *SUM 
C 

c EMISSION  OF  NON-EQUIVALENT  INTERVALLEY  PHONON( CENTRAL  TO  SAT.) 

210  IF(EF( 8) . GT. EZERO ( 2) ) GO  TO  215 

L(8)=0.0 
GO  TO  230 

215  IF(V.EQ.I)  GO  TO  220 
L( 8) =0. 0 
GO  TO  230 

220  MAX (8) =SQRT(EF(8) /EZERO (2) ) 

SUM=0.0 


198 


211 

C 

C 

C 

230 


235 


240 


DO  211  JM=1 ,MAX(8) 
IF(JM.  EQ.N1)  THEN 
AA=3 . 0 
ELSE 
AA=2. 0 
END  IF 
SUM=SUM+AA 
CONTINUE 
L(8)=C7*EM2*SUM 


L(9)-09o'GT'EZEi!0(2))  G°  T°  235 
GO  TO  250 

IF(V. EQ. 1 ) GO  TO  240 
U 9)  =0.0 
GO  TO  250 

MAX ( 9) =SQRT ( EF ( 9)/EZERO(2)) 
SUM=0.0 


wuxvalujt  INTERVALLEY  PH0N0N(  CENTRAI 


DO  231  JM=1  ,MAX(9) 
IF(JM. EQ.  N1)  THEN 
AA=3 . 0 


TO  SAT. 


231 

C 

C 

C 

250 


255 


260 


ELSE 
AA=2. 0 
END  IF 
SUM=SUM+AA 
CONTINUE 
L(  9)  =C8*EM2*SUM 


ABSORPTION  OF  BOB-EQUIVALENT  INTERV4LLET  PHOKON(SAT. 

1 0 T-0 °Q  * GT*  EZER0 c 1 ) ) GO  TO  255 
GO  TO  270 

IF(V.  EQ.2)  GO  TO  260 
L( 10)=0.0 
GO  TO  270 

MAX(10)=sqRt(ef(i°)/eZEro(i)) 
oUM=0 • 0 

DO  251  JM=1,MAX(10) 

IF(JM. EQ.N1)  THEN 
AA=3 . 0 
ELSE 


TO  CEATRAJ 


AA=2.0 
END  IF 
SUM=SUM+AA 
251  CONTINUE 

L( 10)=C9*EM1«SUM 


270  IF(TT.EQ.I)  GO  TO  40 


CALCULATE  SUH  OF  REAL~PROCESS=SCATTERING=RATES====== 


« 
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C== 

280 

290 

C=  = : 
C 

C==: 

295 

296 

C=== 

c 

C = = = : 


c 

300 

310 

320 

330 


SCATT( 1 )=L( 1 )/ GAMMA (V) 

DO  290  K=2 , 1 0 

SCATT (K)=SCATT(K-1 )+L(K)/GAMMA( V) 
CONTINUE 


FLAG  FOR  NEGATIVE  SELFSCATTERING  RATE 


IF( 1 .0-SCATT (101 .LT. 0.0)  THEN 
WRITE ( 5 , 2 95 ) 

WRITE (7, 2 95) 

FORMAT ('  NEGATIVE  SELFSCATTERING  RATE') 
GAMMA ( V) =1 . 1 «GAMMA ( V) 

WRITE (5, 2 96 )V,  GAMMA ( 1 ) , GAMMA (2) 

WRITE (7 ,296) V, GAMMA ( 1 ) , GAMMA (2) 

FORMAT ( " Gamma  is  changed  in  valley  "II 
i"  Gamma(2)  = "E10.4) 

GO  TO  280 
END  IF 


Gamma ( 1 ) = "El  0.4 


CALL  RANDOM  NUMBER. SELECT  SCATTERING  CHANNEL 


R=RANU(0. 0,1.0) 

IF(R.LT. SCATT(I))  GO  TO  300 
IF(R. LT. SCATT(2) ) GO  TO  310 
IFCR.LT. SCATT(3))  GO  TO  320 
IF( R. LT. SCATT ( 4) ) GO  TO  330 
IF(R. LT. SCATT(5) ) GO  TO  340 
IF(R. LT. SCATTC6) ) GO  TO  350 
IF(R. LT. SCATT(7) ) GO  TO  360 
IF(R.LT. SCATT(8) ) GO  TO  370 
IFCR.LT.  SCATT  (9))  GO  TO  380 
IF(R. LT. SCATT ( 10))  GO  TO  3 90 
GO  TO  400 


SET  ENERGY  AFTER  SCATTERING  PROCESS 


EFIN=EF( 1 ) 
GO  TO  420 
EFIN=EF(2) 
GO  TO  420 
EFIN=EF(3) 
GO  TO  410 
EFIN=EF(4) 
GO  TO  410 
EFIN=EF(5) 
GO  TO  410 


340 


200 


350 

360 

370 

380 

390 

400 

410 

C 

C 

C=== 


411 


412 

413 

C4 1 3 

C 

C 

c 

c 


EFIN=EF(6) 
GO  TO  410 
EFIN=EF(7) 
GO  TO  430 
EFIN=EF( 8) 
GO  TO  430 
EFIN=EF(  9) 
GO  TO  430 
EFIN=EF( 10) 
GO  TO  430 
EFIN=EI 
GO  TO  450 
SR=SR+1 


REGISTER  REAL  COLLISION. CALCULATE  K SPACE  POSITION  AFTER  * 

ACOUSTIC, INTERVALLEY  OR  EQUIVALENT  INTERVALLEY  PHONON  SCATTERING 


NACC=NACC+1 

NMAX=SQRT(EFIN/EZERO( V) ) 

PP=0.0 

DO  411  11=1 , NMAX 

COSGAM=H»II *PI / ( LENGTH* 1 . 0 E7  *SQRT ( 2 . 0 *EM*M*EFIN*E ) ) 
SINGAM=SQRT(ABS( 1 .0-C0SGAM*C0SGAM) ) 

PP=PP+SINGAM 
CONTINUE 
R=RANU( 0. 0 ,1.0) 

QQ=0.0 

DO  412  11=1 , NMAX 

QQ=QQ+SQRT(ABS( 1 .0-(II*PI*H)**2/(LENGTH**2*1 .0E1 4*2.0 
it  *EM*M*E*EFIN)))/PP 
IF(R.LT.QQ)  GO  TO  413 
CONTINUE 


N1=II 


UU=  RANU (0.0,1 . 0 ) 

IF(UU. GT.0.5)  THEN 
N1=II 
ELSE 

N1 =-II 
END  IF 

KXF=PI*N1 /LENGTH 
S=RANU(0.0,1 .0) 

PHI=2.0*PI*S 

KZKY2=E*2.0*EM*1 . 0E1 4*M*(EFIN-N1 *N1 *EZERO( V) ) / ( H*H) 
KZF=SQRT(KZKY2)*C0S(PHI) 

KYF=SQRT(KZKY2) *SIN(PHI) 

KRHO=SQRT (KXF*KXF+KYF*KYF ) 

KT  = 1 . 0 E+7  *SQRT ( 2»EM*M*EFIN*E ) / H 
GO  TO  460 
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REGISTER  REAL  COLLISION. CALCULATE  K SPACE  POSITION  AFTER 

DPT  TP  AT  PWOMOXT  CPRT'PrnTvo 


420  AA=EI 

GK( 1 ) =KXI 
GK(2)=KYI 
GK(3)=KZI 
GK(4) =EI 
GK(5) =EFIN 

CALL  P0L0PT (KXI , KYI , KZI , El , EFIN) 

NPOL  OPT  = NPOL  0PT+ 1 
IF(ABS(EI-AA) . GT.  1 .OE-8)  THEN 
WRITE(6 ,421 )EI, AA 
WRITE (7, 421 )EI, AA 

421  FORMAT ( ' JL/EI= 'E20 . 1 4 , ' AA='E20.14) 

DO  424  11=1 ,5 

WRITE  ( 6 , 423 ) GK ( II ) 

WRITE ( 7 , 423 ) GK ( II ) 

423  FORMAT! ' GK='E10.4) 

424  CONTINUE 
STOP 

END  IF 
KXF  =KXI 
KYF=KYI 
KZF=KZI 

KRH0=SQRT (KXF  *KXF+KYF*KYF ) 

KT=SQRT(KRHO‘KRHO+KZF*KZF) 

N1  = NINT (KXF ‘LENGTH/ PI ) 

C 

IF(NI.EQ.O)  THEN 
WRITE (5, 426) 

WRITE (7, 426) 

WRITE ( 5 , 427 ) KXF , LENGTH , PI , KXF ‘LENGTH/ PI 
WRITE ( 7 , 427 ) KXF , LENGTH , PI , KXF ‘LENGTH/ PI 
END  IF 

426  FORMAT ("  ROUNDING  ERROR  !!!!!!") 

427  FORMAT (4(2X,E10.4)) 

C 

IF (EFIN- El.  EQ.  0 . 0)  THEN 
GO  TO  90 
ELSE 

N0PT=N0PT+1 
SR=SR+1 
END  IF 
GO  TO  460 

C CHANGE  VALLEY  PARAMETERS  FOR  INTERVALLEY  PROCESSES. 

430  IF(V.EQ.I)  GO  TO  440 
V=1 

EM=EM1 


o o o 
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440 

GO  TO  410 
V=2 

450 

EM=EM2 
GO  TO  410 
SS=SS+1 

P register  self  scattering  process. k space  position  unchanged  » 


KZF=KZI 
KXF  =KXI 
KYF=KYI 


KRHO=SQRT (KXF*KXF+KYF*KYF ) 
KT=SQRT(KZF*KZF+KRHO*KRHO) 
EFIN=H*H*KT*KT*1 . OE-1 4/ ( E*2*EM»M) 


C 

C 

C 

C 

C 

C= 


CHECK  IF  ELECTRON  IS  SCATTERED  OUT  OF  MESH. IF  SO, REGISTER 
PROCESS  ON  COUNTER  GMAX, AND  PLACE  ELECTRON  ON  EDGE  OF  MESH 
LABEL  90  REPEATS  ITERATIVE  PROCESS  STARTING  WITH  FREE 
ELECTRON  FLIGHT  UNDER  ELECTRIC  FIELD. 


9 

* 

» 

9 


460  Emax1=DMAX1 (Emaxl ,EFIN) 

Eminl =DMIN1 (Eminl , EFIN) 
IF(EFIN.LE.EMAX)  GO  TO  90 
GMAX=GMAX+1 

KT= 1 . 0 E+7  «SQRT ( 2*EM«M*EMAX*E ) / K 
IF(KXF. GT.KT)  THEN 

N1  = NINT ( LENGTH  *KT / PI ) 

KXF =N1* PI /LENGTH 
END  IF 

IF(KXF.LT.-KT)  THEN 
N1 = NINT ( LENGTH *KT /PI ) 

KXF =N1«PI/ LENGTH 
END  IF 

AA = KT  *K  T -KXF  *KXF 

BB=KYF*KYF+KZF*KZF 

AA=SQRT(AA/BB) 

KZF=KZF*AA 

KYF=KYF*AA 

KRHO= SQRT ( KXF  *KXF  +K YF#K YF ) 

GO  TO  90 


C====r====: 


FINAL  CALCULATION  OF  DRIFT  VELOCITY-VEL , TIME  SPENT  IN  EACH 
VALLEY- VTIM( V) ,MEAN  ENERGY-MNE 


C 

470  IF(IDIF.NE.O)  GO  TO  471 
VAVG=0 . 0 
V2AVG=0 . 0 
V3AVG=0.0 


9 

9 
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V4AVG=0.0 
471  IDIF=IDIF+1 

NPRINT=NPRINT+1 

ET0T=ET0T- ( KZI*KZI-KZF*KZF )*10.0*H/(  2*EM*M ) 
KZT0T=KZT0T-ABS(KZI-KZF ) 

KZI=KZF+(E*F ( J) * 1 . OE+1 8/  H) « ( (FT-TT0T)+TIME ) 
ET0T=ET0T+(KZI*KZI-KZF*KZF)*10.0*H/(2*EM*M) 
KZT0T=KZT0T+ABS(KZI-KZF ) 

VI=ET0T/KZT0T 

MNE=MNE/KZT0T 

VAVG=( VAVG* (IDIF-1 )+VI)/IDIF 
V2AVG= ( V2 AVG* ( IDIF-1 ) +VI*VI ) /IDIF 
V3 AVG= ( V3 AVG* ( IDIF-1 ) +VI*VI*VI ) /IDIF 
V4AVG=(V4AVG*(IDIF-1 )+VI*VI*VI*VI)/IDIF 
DIF=0.5  *FT*( V2AVG-VAVG*VAVG) 

DIF4=SQRT(ABS(  V4  AVG-4*V3  AVG*VAVG+6  *V2 AVG*VAVG**2-3*VAVG»*  4)  / 1 2 
1.0)*FT 

CHECK=V3AVG-3‘V2AVG*VAVG+2*VAVG**3 

ET0T=0.0 

KZTOT=0.0 

TT0T=0.0 

480  FORMAT(2X,I3,2(2X,E10.4) ,3(2X,I8)  ,3X,  12 ,3X,  18 ,3X, E10 .4) 

490  FORMAT (/,/, " Electric  Field  Strength  = "E8.2,"  kV/cn."/) 

520  FORMAT( IX, nFlt#n2X, "Ave  Vel"5X,"Dif  Coef"7X,"#  opt"5X,"#  acc" 

# 4X , se  scatt") 

530  FORMAT ("  Min.  Energy  = "E10.4,"eV  Max.  Energy  = nE10 . 4 , "eV") 

IFCIDIF.EQ. 1 ) THEM 
WRITE(5,490)F(J) 

WRITE (7 , 4 90 )F( J) 

WRITE (5, 520) 

WRITE (7, 520) 

END  IF 

WRITE  ( 5 , 4 80 ) IDIF , VAVG , DIF , NOPT , NACC , SS , J , NPOLOPT , MNE 
WRITE  ( 7 , 4 80 ) IDIF , VAVG , DIF , NOPT , NACC , SS , J , NPOLOPT , MNE 
IFCIDIF.EQ. NDIF)  THEN 

WRITE ( 5 i 530 ) Eminl ,Emax1 
WRITE(7,530)Emin1 ,Emax1 
END  IF 


NPRINT=0 

IFCIDIF.LT. NDIF)  GO  TO  90 
IDIF=0 

DO  500  II=-20 ,20 

WRITEC7, 510)11, TLEVCl, II), TLEVC2, II) 
510  FORMAT ( 2X , 14 , 2 ( 2X , El  0 . 4 ) ) 

500  CONTINUE 
J=J+1 

IF(J. NE.NF+1)  GO  TO  80 
999  STOP 
END 


o o o o 
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SUBROUTINE  POLOPT ( RKXI , RKYI , RKZI , REI , REFIN) 

SUBR.  TO  CALCULATE  THE  WAVE  VECTOR  AFTER  A POLAR  OPTICAL  » 
SCATTERING.  a 


DOUBLE  PRECISION*6  RKXI, RKYI, RKZI, REI, REFIN 

REAL*6  RM, REM, RE, RH, RLENGTH 

COMMON/POLOP/RM, REM, RE, RH, RLENGTH 

REAL *6  PI,  A1  , A2 , A3  , U, BETA, R,  BNORM, ANORM,  W,E1  ,E2 

REAL *6  P1,P2,P3,P4,TEST,RNU,F,PSI0,B1  ,B2,B3,BN 

REAL*6  RNU1 , RNU2 , RNU3 , RNU4 ,PSI( 100 ) 

INTEGER«6  II , NMAX, J, I, K, KO, II , N, NPLANE 
PI=ATAN(1.0)»4.0 


C NORMALIZATION 


REI=RLENGTH»RLENGTH*REI/(PI»PI) 
REFIN=RLENGTH'RLENGTH*REFIN/  ( PI'PI ) 

A1= RKXI 'RLENGTH/ PI 
A2= RKYI 'RLENGTH/ PI 
A3  = RKZI'RLENGTH/ PI 
C 

EX=2 . 0 *SQRT ( REFIN'REI ) / ( ( SQRT ( REI ) -SQRT ( REFIN) ) " 2) 
U=RANU( 0.0 ,1.0) 

BETA  = ( 1 . 0+EX-( 1 . 0+2 . 0 »EX) »»U) /EX 
BNORM=SQRT (2.0 'RM'REM'REFIN'RE ) » 1 . 0 E7 / RH 
ANORM=SQRT ( A 1 » A 1 +A2 'A2+A3  »A3 ) 

W=BNORM*SQRT( 1 . O-BETA'BETA) 

DO  70  11=1,100 
PSI(II)=1000.0 

70  CONTINUE 
E1=1.0+A1*A1/(A3'A3) 

El =SQRT(E1 ) 

E2=1 . 0+A1  »A1  / ( A3  'A3 ) +( A 1 »A  1 +A3  'A3 ) " 2/  ( A2  »A2  »A3  «A3 ) 
E2=SQRT(E2) 

NMAX=BNORM 
IFCNMAX.GT.20)  THEN 
WRITE (5,71) 

WRITE (7, 71 ) 

71  FORMAT (•  NMAX  EXCEEDS  20') 

STOP 

END  IF 
J=0 

DO  10  1=1 , NMAX 
PI =W/E1 

P2=W*A1/(E2*A3) 

P3=1 .0*I-BN0RM»BETA»A1 /ANORM 
P4=SQRT(P1'P1+P2»P2) 

TEST=P3/P4 

IF(ABS(TEST ) . LE. 1.0)  GOTO  20 
GOTO  10 


O O o 
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C 

C 

C 

C 

C 


20  RNU1 =AC0S(P1/P4) 

RNU2=-AC0S(P1/P4)+2.0“PI 

RNU3=ASIN(P2/P4) 

IF(RNU3  .LT.  0.0)  RNU3=RNU3+2.0«PI 
RNU4=PI-ASIN(P2/P4) 

IF(RNUi<  .LT.  0.0)  RNU4=2.0«PI+RNU4 
IF(ABS(RNU1 -RNU3 ) .LT.  1.0E-8)  THEN 
RNU=RNU1 
ELSE 

IF(ABS(RNU1-RNU4)  .LT.  1.0E-8)  THEN 
RNU=RNU1 
ELSE 


RNU=RNU2 
END  IF 
END  IF 

PSI(2*I)=ASIN(P3/P4)-RNU 
PSI(2*I+1 )=PI-RNU-ASIN(P3/P4) 
J= J+2 
10  CONTINUE 

R=RANU( 0.0 ,1.0) 

F=0.0 

IF(J.EQ.O)  GOTO  998 
DO  30  K=1,J 
F=F+1 .0/ J 
IF(R.LE.F)  GOTO  40 
30  CONTINUE 
40  KO=K 


11=0 


DO  50  N=1 , 1 00 

IF(PSI(N) . NE.  1000.0)  GOTO  60 
GOTO  50 
60  11=11+1 


IF(II.EQ.KO)  GOTO  80 
50  CONTINUE 
80  PSIO=PSI(N) 

IF(  N.  LT.  0 ) GOTO  11 
GOTO  12 

11  NN=N-1 
NPLANE=NN/2 
GOTO  13 

12  NPLANE=N/2 
IF(NPLANE.EQ.O)  STOP 

13  B1=1.0*NPLANE 


B2=BETA»BNORM»A2/ANORM-W»COS(PSIO)*(A3*A3+A1»A1)/(E2»A2*fl^1 

psi°  >,k'nE"k3>  *u*cos ' 1 > '« 

DENORMALIZATION  ======================== 


RKXI=B1 “PI/RLENGTH 
RKYI=B2*PI/RLENGTH 
RKZI=B3  *PI / RLENGTH 

refin=pi»pi«refin/(rlength»rlength) 

REI = PI »PI »REI / ( RL  ENGTH “RLENGTH ) 
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GOTO  999 

998  REI= REI *PI »PI / ( RLENGTK*RLENGTH ) 
REFIN= REI 

RKXI=A1 *PI/RLENGTH 
RKYI=A2*PI/RLENGTH 
RKZI=A3*PI/RLENGTH 

999  RETURN 
END 


APPENDIX  B 

REFLECTION  COEFFICIENT 


It  is  proven  here  that  a parasitic  capacitance  parallel  to  the 
device  under  test  or  a parasitic  inductance  in  series  with  the  device  do 
not  affect  the  noise  temperature  measurements.  Consider  the  two  cir- 
cuits drawn  in  Fig.  B.l.  In  this  figure  the  device  is  represented  by 
Zj),  the  inductor  by  Z^,  the  capacitance  by  Z^.,  and  the  characteristic 
impedance  of  the  stripline  by  Zq . The  power  per  unit  bandwidth  radiated 
into  the  stripline  will  be  calculated  in  both  cases:  v,v”  / Z^Af  . 

Case  a 

The  voltage  across  Zq  is 


v 


1 


Z0  + ZC  + ZD 


e 

n 


(3.1) 


1/2 

where  e^  = {4  k^T^  Re[ZQ]Af}  ' , kg  is  Boltzmann's  constant,  and  Tn  the 
noise  temperature.  The  reflection  coefficient  r equals 


ZT  ~ Z0 
ZT  +Z0 


ZL  + ZD  - Z0 
ZL  + ZD  + Z0 


(3.2) 


Combining  eqs.  (B.l)  and  (B.2),  one  finds 


v 


1 


( B.3) 


The  radiated  power  P into  the  stripline  is 
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Fig.  B.l.  Equivalent  circuits  for  (a),  an  impedance  in 
series  with  the  device,  and  (b),  an  impedance 
parallel  to  the  device. 
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P 


* 


V1V1 

Vf 


U - rr 

4 Z0Af 


* 

e e 
a n 


4 k T 
3 n 


Re  [ Z ] Af  . 
D 


(3.4) 


Under  the  assumption  that  Zg  is  real  and  by  using  efj.  (3.2)  to  express 
Z3  in  terras  of  T,  ZQ,  and  Z^,  eq.  (3.4)  will  read: 


?-kETJ‘  - r|2(^Re[-ZtJ  +R.ifi-f]} 

9 9 R*[Zr] 

= STn(1  - lI’l')  “ VJ1  - Fl  -zT"  • (3-5) 


When  ZL  is  purely  imaginary,  the  radiated  power  per  unit  frequency  is 
equal  to 


p - W1  - 


(3.6) 


Case  b 

In  this  example  the  noise  source  is  a current  source  and  the  for- 
malism will  be  in  admittances  rather  than  impedances.  The  current 
source  i equals 

‘n  ' !4  Vn  RelYDl4£>1/2  (B.7) 

where  YQ  = 1/Zg. 

It  follows  that 


= Vyd 


+ Y, 


Y0) 


-1 


(3.3) 


The  reflection  coefficient  T is 
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r = 


ZT  - ZQ 
ZT  +Z0 


(3.9) 


where 


yt  zt  zd//zg  (yd  + V 


(3.10) 


Equations  (B.9)  and  (3.10)  yield 


i - r 

Y + Y = Y - 

d c o i + r ’ 


(3.11) 


which  substituted  in  eq.  (3.8)  gives 


v,  = 


n 1 + r 


! Y0  2 • 


(3.12) 


The  radiated  power  per  unit  bandwidth  equals 


P = 


knT 
B n 


|1  + r|Z{Re[Y0  f-=-£]  - Re[Y  ]} 


(3.13) 


Under  the  assumption  that  YQ  is  real  and  Yc  is  imaginary,  the  radiated 
power  is 


P = kBTn(1  “ In*  • 


(3.14) 
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